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CHAPTER  I 
INTRODUCTION 


A.  Research  Goals 

Subsurface  radar  detection  and  identification  of  geological  and 
man-made  structures  Is  an  area  of  current  importance.  Examples  are: 
location  of  utility  pipes  such  as  plastic  and  metallic  gas  pipes  and 
water  p1pes[l ,2,3] , location  of  voids  and  tunnel s[4],  anthropology 
mapp1ng[5],  and  possible  exploration  of  energy  sources  such  as  oil . 
gas  and  coal.  Yet,  almost  all  of  the  work  done  In  this  area  was 
directed  toward  target  detection.  Pew  attempted  the  problem  of  tar- 
get identification.  Subsurface  target  identification  Is  a problem 
far  more  severe  than  the  Identification  of  aerospace  targets  by  con- 
ventional radars  where  the  target  can  literally  be  seen  and  the  class 
of  false  targets  is  limited  In  scope.  Underground  there  are  varltles 
of  unknown  false  or  undesired  targets  to  complicate  the  task.  Further- 
more, the  medium  involved,  i.e,,  the  ground,  is  usually  lossy,  inhomo- 
geneous and,  most  of  all,  electrically  weather-dependent.  These  pro- 
blems, together  with  the  presence  of  the  air-ground  Interface  makes 
the  task  of  subsurface  target  identification  truly  formidable.  It  is 
for  these  reasons  that,  to  date,  there  Is  no  single  technique  or 
system  capable  of  identifying  subsurface  targets  in  real  time. 

In  this  study,  a technique  for  subsurface  target  Identification 
Is  developed  and  extensively  tested  with  real  radar  measurements  col- 
lected using  a video  pulse  radar[l-4]  under  different  conditions  (I.e., 
different  ground  conditions,  different  antennas,  etc.).  This  techni- 
que is  Implemented  with  a "first- generation"  microcomputer  system  to 
demonstrate  the  feasibility  of  real-time  subsurface  target  Identifica- 
tion. 

The  technique  used  in  this  study  characterizes  subsurface  targets 
by  their  complex  natural  resonances[7-ll],  which  are  extracted  directly 
from  the  processed  time  domain  waveforms  via  Prony's  Mathod [12-1 5] . 

A predictor-correlator[40]  use5:  the  difference  equation  coefficients 
associated  with  these  complex  resonances  as  discriminants  to  generate 
a correlation  coefficient  for  target  identification.  This  characteri- 
zation and  identification  method  Is  attractive  for  It  characterizes  the 
response  of  a target  by  a set  of  complex  numbers  which  Is  Independent 
of  the  pulse  radar  location.  Furthermore,  the  complex  resonances  and 
the  difference  equation  coefficients  are  pre-determlned , thus,  only  sim- 
ple algebraic  operations  are  Involved  In  calculating  the  correlation 
coefficient  for  a real-time  Identification  decision. 
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[3 . Related  Research  In  Subsurface 

Target  Characterization  and 

Identification 

Electromagnetic  techniques  have  been  used  successfully  for 
many  years  for  probing  the  earth.  Keller  and  Freschnecht[25]  give 
an  excellent  summary  of  these  procedures.  A recent  summary  of  sub- 
surface probing  techniques  Is  given  in  a report  by  D.C.  Gates,  et 
a l.[  26  ] . 


Perhaps  the  earliest  documentation  of  a subsurface  electro- 
nagnetic  radar  system  is  contained  in  a patent  issued  in  1937  in 
which  an  electrical  analog  of  seismic  systems  is  descr i bed [ 16] . 

There  is,  however,  no  mention  of  successful  Implementation  of  such 
a radar.  There  have  been  attempts  to  use  bistatic  radar  configura- 
tions of  this  type  but  the  results  have  generally  not  been  highly 
successful [17,18] . The  reason,  recognized  by  Horton[19],  is  that 
the  tail  of  the  pulse  coupled  directly  from  the  transmit  to  the  re- 
ceive antenna  occurs  "just  at  the  time  when  the  maximum  of  the  re- 
flected pulse  (from  a burled  target)  must  be  accurately  timed. 

This  coincidence  tends  to  ruin  the  measurement". 

A significant  result  In  video  pulse  technology  is  described  In 
a patent  by  Lerner  in  which  a video  pulse  system  is  used  for  more 
moderate  depths[20].  Lerner's  scheme  differed  from  the  earlier  patent 
in  that  the  same  antenna  was  used  for  both  transmitting  and  receiving. 
Lerner  introduced  a combination  of  TR,  ATR  and  Hybrids  to  separate 
the  transmitted  pulse  and  the  received  target  signal  . 

In  this  study,  a crossed-dipole  antenna  system  Is  used  to  In- 
corporate this  function  into  the  antenna  Itself,  l.e.,  transmit- 
receive  isolation  is  achieved  by  isolating  the  antennas  themselves. 

The  crossed  dipole  is  an  orthogonal  dipole  pair,  one  horizontal  dipole 
for  transmission  and  another  orthogonal  horizontal  dipole  for  recep- 
tion, which  provides  substantial  reduction  of  the  primary  (directly 
coupled)  signal  on  the  receiving  antenna.  Many  measurements  have 
been  made  on  a variety  of  shallow  targets  (less  than  15  m)  using  these 
com epts[1 -4,21 ,22] . Targets  include  geological  structures  such  as 
faults,  joints,  sink  holes  and  man-made  structures  such  as  pipes.  A 
commercial  unit  for  pipe  detection  based  upon  the  research  and  design 
work  at  the  ElectroScience  Laboratory,  and  designated  as  Terrascan, 
is  being  produced  by  Microwave  Associates  Inc. [6].  The  pulse  radar 
used  in  this  study  for  subsurface  target  identification  'Is  a Terra- 
scan-like  radar  system. 

A major  reason  for  the  success  of  the  characterization  and  identi- 
fication procedures  discussed  in  this  dissertation  lies  in  the  Improve- 
ments In  the  antenna  system.  The  original  cross-type  antenna  struc- 
ture developed  at  the  ElectroScience  Laboratory  for  subsurface  radar 
applications  was  basically  a crossed  bowtle  geometry  wrapped  aroung 
a sphere[21],  but  this  was  subject  to  noise.  The  crossed  bowtie 
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evolved  Into  various  planar  crossed- dipole  arrangements  as  used  by 
Moffatt[4].  These  arrangements  were,  however,  too  awkward  for  use 
in  real-time  on-location  target  identification.  Later,  Young[l,3] 
introduced  the  loaded  folded  dipole  geometry  of  the  Terrascan  system. 
Tribuzi[66]  Improved  this  by  Introducing  the  loaded  folded  bowtle 
configuration.  Wald[59]  further  improved  the  electrical  character- 
istics of  this  structure  by  eliminating  part  of  the  supporting  struc- 
ture. He  also  constructed  the  small  antenna  used  in  a later  part  of 
this  dissertation  for  identification  of  mines.  The  design  of  the 
smaller  antenna  was  dictated  by  the  results  obtained  in  this  dis- 
sertation since  its  purpose  was  to  shift  the  antenna  resonance  to 
more  nearly  coincide  with  those  of  the  mine-like  target. 

One  of  the  first  studies  initiated  in  the  development  of  the 
pulse  radar  system  with  the  crossed-dipole  antenna  for  subsurface 
target  detection  and  identification  was  the  detection  and  identifi- 
cation of  TDMB  m1nes[28-35].  This  study  included  efforts  In  the 
development  of  antenna  systems  and  techniques  to  extract  the  char- 
acteristic spectra  of  the  electromagnetic  fields  scattered  by  the 
TDMB  mine.  In  1970,  Su 1 1 1 van[21 ] investigated  the  feasibility  of 
using  the  characteristic  real -frequency  resonances  of  subsurface 
targets  in  the  Identification  of  simple  burled  objects.  The  system 
used  in  Sullivan's  study  was  a video  pulse  radar  with  a crossed- 
polarized  antenna  system.  A similar  system  using  crossed  dipoles 
was  later  used  by  Moffatt,  et  a 1 . [4]  in  the  probing  of  man-made  and 
geological  subsurface  targets.  The  subsurface  video  pulse  radar 
system  was  then  modified  and  developed  to  be  the  existing  Terrascan 
system  in  a study  to  detect  gas  pi pes [1 -3] . To  automate  the  Terra- 
scan system  for  automatic  target  identification,  Chan[22]  Investi- 
gated a matched  filter  technique  for  automatic  identification  of 
plastic  pipes  using  a Terrascan-1 1 ke  radar  system. 

Other  methods  have  been  employed  in  the  detection  and  Identi- 
fication of  subsurface  targets.’  In  particular,  various  techniques 
of  Pattern  recogni t 1 on[37 -39]  were  used  by  Echard,  et  al.  for  the 
detection  and  Identification  of  buried  mi nes[36] . 

This  study  investigated  the  possibility  of  using  in  situ  target 
complex  natural  resonances  to  characterize  and  identify  subsurface 
targets.  The  basic  method  was  first  introduced  by  HI  11 [40]  In  the 
detection  of  targets  near  the  surface  of  the  earth,  and  later  used 
by  Chan,  et  al . [23,48-50]  in  the  characterization  and  Identification 
of  subsurface  targets. 
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C . Structure  of  This  Report 

The  structure  of  this  report  is  as  follows: 

Jn  chapter  II,  the  basic  radar  measurement  procedure  is  per- 
sented.  In  addition,  preliminary  signal  processing  is  discussed 
since  certain  preprocessing  does  improve  the  target  identification 
results. 

In  Chapter  III,  we  present  a method  for  extracting  from  response 
data  records  (l.e.,  the  time-domain  waveforms)  the  complex  natural 
resonances  associated  with  the  targets.  This  method,  known  as  Prony's 
method,  is  outlined  and  applied  to  extract  target  resonances  from  ™e 
backscattered  waveforms  in  the  time  domain. 

In  Chapter  IV,  the  predictor-correlator  method  for  target 
identification  is  discussed  and  applied  to  the  waveforms  collected 
in  this  study.  Detail  identification  statistics  are  given. 

In  Chapter  V,  the  effects  of  radar  bandwidth  on  the  character- 
ization and  identification  method  are  studied. 

In  Chapter  VI,  the  effects  of  target  size  and  depth  on  the 
characteri zation  and  identification  method  are  discussed. 

In  Chapter  VII,  we  direct  attention  to  the  detecion  and  identi- 
fication of  mine-like  targets  in  practical  situations.  Improvements 
on  the  pulse  radar  system  are  made  for  the  implementation  of  a port- 
able, real-time  on-location  subsurface  target  identification  radar. 


In  Chapter  VIII,  we  discuss  the  implementation  of  the  sub- 
surface identification  radar  as  a microcomputer  system.  Detailed 
prodecures  of  the  implementation  are  presented.  Real-time  target 
identification  results  using  the  microcomputer  system  are  given. 

In  Chapter  IX,  a method  for  automatic  tuning  of  the  identi- 
fication radar  to  the  ground  condition  in  real  time  is  discussed. 
This  method  is  simple  and  can  be  easily  incorporated  into  the  micro- 
computer system  for  real-time  subsurface  target  identification. 

In  Chapter  X,  major  achievements  accomplished  in  this  work  are 
summarized.  Conclusions  and  recommendations  are  made. 
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CHAPTER  II 

MEASURED  AND  PROCESSED  WAVEFORMS 
FROM  THE  SUBSURFACE  TARGETS 


A. 


The  objectives  of  this  chapter  are  the  following: 

1.  To  give  a description  of  the  subsurface  pulse  radar 
and  the  subsurface  targets  selected  for  this  study  and 
to  summarize  the  procedures  taken  to  measure  the  back- 
scattered  waveforms  from  these  targets.  Raw  (unpro- 
cessed) waveforms  from  the  targets  are  shown. 

2.  To  summarize  the  signal -procecsing  techniques  used 
to  partially  suppress  noise  and  clutter  in  the  raw 
waveforms.  Processed  waveforms  are  given. 


R.  Subsurface  Electromagnetic  Video 

Pulse  Radar  System 

The  video  pulse  radar  system  used  to  collect  measurements  for 
this  study  basically  consists  of  three  components:  the  energy  source, 
the  antenna  system  for  signal  transmitting  and  receiving  and  the 
receiver  for  signal  processing.  The  design  of  these  components  is 
dictated  by  the  electrical  properties  of  the  ground,  the  depth  of 
the  target  of  interest  as  well  as  the  target,  clutter  and  noise 
characteristics.  A picture  of  the  Terrascan-l i ke  subsurface  pulse 
radar  used  in  this  study  together  with  a basic  block  diagram  is 
shown  in  Figure  1.  The  basic  components  are:  the  impulse  generator, 
the  crossed-dipole  antenna  system  and  the  receiver.  Basic  operation 
is  as  follows.  The  impulse  generator  transmits  short  pulses  of  energy 
through  the  transmit  antenna  into  the  ground.  The  presence  of  a 
target  scatters  the  incident  energy  toward  the  receive  antenna.  This 
scattered  energy  is  received  as  a sampled  time-domain  waveform  for 
target  characterization  and  identification. 

In  the  current  study  a short  video  pulse  of  approximately  150 
ps  duration  (at  3 dB  points)  and  a nominal  1Q00V  peak  amplitude 
wo 5 used  (see  Figure  2).  This  pulse  duration  is  much  shorter  than 
those  used  in  conventional  radar  practice.  Furthermore,  a con- 
ventional radar  has  a few  percent  bandwidth  about  its  carrier 
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Figure  1.  The  subsurface  pulse  radar  and  its  block  diagram. 
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frequency  whereas  this  video  pulse  output  spectrum  spreads  from 
essentially  dc  (repetition  rate  = 256  Hz)  to  beyond  3 GHz.  It  is 
this  broad  band  of  frequencies  inherent  in  this  narrow  pulse  that 
makes  tarqet  identification  a possibility,  i.e.,  the  scattered  fields 
from  the  targets  can  be  sampled  over  a very  broad  frequency  band  and 
each  sample  contains  information  about  the  targets.  The  use  of  such 
a narrow  pulse  also  has  a second  substantial  advantage  for  the  de- 
tection and  identification  of  shallow  targets,  in  that  the  transmit- 
ted pulse  magnitude  has  fallen  to  a low  value  before  the  pulse  re- 
flected from  the  target  returns  to  the  antenna.  This  "time  isola- 
tion" effectively  minimizes  the  width  of  the  radar  "dead  zone"  and 
enables  the  scattered  pulse  to  be  observed  over  a wide  time  span. 

A second  form  of  isolation  exists  in  the  choice  of  the  antenna  system. 
The  pulse  radar  uses  a pair  of  crossed,  loaded,  folded  dipoles  with 
0.6m  (2  feet)  long  arms  lying  flush  with  the  ground  surface  (see 

Figure  1).  The  crossed-dipole  antenna  system  achieves  substantial 
isolation  between  transmit  and  receive  antennas.  For  a perfectly 
orthogonal  pair  and  no  target  perturbation,  the  transmitted  pulse 
would  not  be  observed  on  the  receive  antenna.  In  practice,  antenna 
isolation  on  the  order  of  60  dB  below  the  nominal  pulser  voltage  is 
achieved  routinely.  Such  isolation  further  minimizes  the  width  of 
the  "dead  zone"  and  is  essential  for  shallow-depth  target  identifica- 
tion. The  dipoles  are  heavily  loaded,  with  both  resistors  and 
incorporated  absorber  in  the  antenna  to  reduce  multiple  reflections 
and  consequently  reduce  pulse  distortion  caused  by  the  antenna.  The 
crossed-dipole  antenna  system  has  two  additional  advantages  for  the 
identification  of  subsurface  targets.  First,  being  a cross-polarized 
system  it  is  insensitive  to  reflections  from  layers  which  are  parallel 
to  the  antenna  arms.  An  important  example  is  the  ground  surface 
whose  reflection  of  the  incident  pulse  energy  would  produce  extraneous 
signals  in  other  non-orthogonal  systems.  Second,  received  waveforms 
obtained  from  objects  which  have  no  symmetry  with  respect  to  the 
antenna  arms  go  through  a polarity  reversal  as  the  crossed-dipole 
antenna  system  is  rotated  (about  its  vertical  axis)  by  90°.  This 
feature  represents  a valid  method  by  which  a target  can  be  separated 
from  an  extended  no-target  echo  which  is  introduced  by  multiple  re- 
flections on  the  antenna  structure^]. 

Typical  raw  time-domain  waveforms  received  with  the  antenna 
oriented  at  0"  and  at  90"  over  the  center  of  a plastic  mine-like 
target  are  shown  in  Figure  3.  Several  waveform  features  can  be 
described.  The  first  sharp  impulsive-type  portion  of  the  waveform 
is  due  to  direct  coupling  between  the  transmit  and  receive  antennas. 
From  the  amplitude  of  the  coupling  signal  and  the  pulser  output,  it 


*A  smaller  antenna  (0.15  m)  is  later  used  for  a more  specific  purpose 
of  mine  idonti fi cation. 


is  seen  that  about  60  dB  isolation  is  achieved.  Note  that  this  sig- 
nal does  not  change  significantly  with  antenna  orientation  and  thus 
can  be  removed  by  forming  the  difference  of  the  two  waveforms.  The 
coupling  signal  is  a source  of  clutter  but  it  is  also  useful  as  a 
time  reference  for  target  depth  and  range,  since  it  occurs  at  essenti- 
ally the  time  the  source  impulse  is  radiated  from  the  feed  terminals 
of  the  transmit  antenna.  The  next  feature  of  the  waveform  beyond 
the  impulse  is  the  random  clutter  due  to  the  ground  surface  irregu- 
larities directly  beneath  the  antenna.  Because  of  the  short  duration 
of  the  impulse  and  good  antenna  design,  this  clutter  feature  dies  out 
in  a few  nanoseconds.  Thus,  only  a small  portion  of  the  clutter  over- 
laps the  return  from  the  mine-like  target.  The  signal  from  the  mine- 
like target  appears  to  be  rather  strong.  Furthermore,  it  reverses 
polarity  when  the  antenna  is  rotated  by  90°,  thus,  its  amplitude  will 
double  in  the  difference  waveform.  The  mine-like  target  signal  ex- 
tends through  a time  window  of  approximately  30  ns  and  falls  to 
a negligible  level  at  the  time  the  balun  reflection  occurs.  The  balun 
is  necessary  for  connecting  the  unbalanced  impulse  generator  to  the 
balanced  dipole  antenna.  The  impedance  mismatch  at  this  connection 
is  the  source  of  the  balun  reflection.  The  balun  reflection  limits 
the  width  of  the  reflectionloss  time  window  of  the  system.  In  the 
present  system,  the  width  of  the  reflectionless  time  window  is  36.5 
ns,  which  turned  out  to  be  wide  enough  for  the  identification  of  the 
shallow  subsurface  targets  considered  in  this  study.  For  a wider 
window,  one  can  lengthen  the  delay  cable  at  the  balun-antenna  con- 
nection. One  can  also  greatly  suppress  the  effects  of  the  balun  re- 
flection by  shortening  the  length  of  the  delay  cable  at  the  balun- 
antenna  connection.  In  this  case,  the  balun  reflection  would  occur 
in  the  time  region  where  the  target  signal  is  much  higher  in  amplitude. 

In  a later  section,  we  describe  a smaller  antenna  which  was  used  for  5 

improved  target  identification  performance.  This  small  antenna  was  | 

built  with  the  delay  cable  shortened  and  the  balun  structure  placed 

almost  at  the  antenna  feed  points.  The  balun  reflection  can  be  com-  I 

pletoly  eliminated  if  a balanced  pulser  is  made  available  in  the 

future,  l.l  iminatlon  of  the  balun  would  also  yield  a narrower  trans-  i 

in i t. ted  pulse  due  to  less  cable  loss  and  dispersion.  ; 

An  unprocessed  time-domain  waveform  is  present  in  the  receiver  ; 

for  target  characterization  and  identification.  The  structure  of  the  ; 

receiver  basically  consists  of  a sampling  oscilloscope  for  signal  j 

reception,  a signal -processing  unit  for  clutter  and  noise  reduction  ; 

and  .i  unit  for  target  characterization  and  identification.  In  this  j 

study,  because  of  the  flexibility  it  offered,  a general-purpose  digi-  j 

tal  computer  was  first  used  to  control  the  sampling  oscilloscope  and  “ 

to  implement  the  processing,  character i zati on  and  identification  units.  ] 

After  all  of  the  target  characterization  and  identification  procedures  J 

had  been  established,  much  of  the  flexibility  was  discarded  and  a re-  J 

latively  simple  system  designed  for  target  detection  and  identifica-  1 

tion.  Such  a system  was  implemented  with  a microcomputer  for  target  ] 

identification  in  real  time.  Discussion  of  the  microcomputer  system  I 

is  presented  in  Chapter  VIII.  i 
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C . The  Subsurface  Targets 

Five  targets  of  similar  size  were  buried  at  the  same  depth  of 
5 c.n  (2  inches,  measured  from  the  qround  surface  to  the  nearest  tar- 
get surface).  Figure  4 shows  the  geometry  of  the  targets. 

All  targets  were  buried  in  the  backyard  of  the  ElectroScience 
Laboratory  at  points  where  the  ground  is  known  to  be  relatively  un- 
disturbed (i.e.,  free  of  other  objects).  The  average  dc  conductivity 
within  30  cin  (12  inches)  of  the  ground  surface  measured  at  the  target 
sites  ranged  from  30  rnS/m  for  wet  ground  to  about  20  mS/m  for  dry 
ground  and  10  mS/m  for  icy  ground.  The  relative  dielectric  constant 
measured  at  approximately  100  MHz  ranged  from  25  for  wet  ground  to 
16  for  dry  ground  and  9 for  icy  ground. 

Since  a goal  of  this  study  was  to  achieve  separation  of  a mine- 
like target  from  other  (false)  targets  using  radar  data,  a major 
effort  was  directed  toward  the  study  of  the  mine-like  target.  The 
method  of  identification  developed  here  can  easily  be  adapted  to 
systems  in  which  other  targets  are  considered  as  desired  targets  or 
targets  to  be  separated. 

Backscattered  waveforms  were  obtained  using  the  subsurface 
pulse  radar  system.  Measurements  were  made  at  different  ground  loca- 
tions with  respect  to  the  various  targets.  Locations  of  the  antenna 
center  for  these  measurements  are  shown  as  dots  in  Figure  5.  At 
each  location  two  backscattered  waveforms  were  obtained  using  two 
different  antenna  orientations,  one  of  which  was  a 90°  rotation  with 
respect  to  the  other.  A standard  antenna  orientation  used  in  obtain- 
ing the  measurements  is  shown  in  Figure  5.  Measurements  were  obtained 
with  the  antenna  center  vertically  above  the  center  and  edges  of  the 
targets.  Beyond  the  target  edges,  measurements  were  made  at  the 
regular  interval  of  15  cm  (6  inches). 

Data  accumulation  was  started  in  early  June  19/7  and  continued 
through  early  April  1979.  During  this  period  the  ground  condition 
changed  from  wet  t.o  dry  and  to  icy.  Data  were  obtained  for  each 
ground  condition  to  gauge  the  effects  of  the  changing  ground  condition 
on  the  characterization  and  identification  of  the  subsurface  targets. 
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gurc  5.  Measurement  locations  and  antenna  orientation. 
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D.  Raw  Measured  Waveforms 


All  waveforms  collected  by  the  Terrascan-1 i ke  radar  used  in 
this  study  consist  of  256  samples  In  a time  window  of  50  ns.  The 
(hardware)  basic  sampling  period  Tg  is  0.2  ns,  giving  a sampling 
frequency  of  5.12  GHz*. 

There  are  three  possible  classes  of  signals  present  in  these 
raw  waveforms. 

1.  Noi se:  noise  refers  to  extraneous  signals  which  are 

not  in  any  way  related  to  the  radar  source  signal. 
Examples  are  thermal  noise,  interference,  etc. 

2.  Cluster:  clutter  refers  to  extraneous  signals  which 

are  related  to  the  radar  source.  Examples  are 
transmit-recei ve  coupling,  reflection  from  ground 
surface  irregularities,  and  echoes  from  objects 
other  than  the  desired  target. 

3.  Desi red  Signal : desired  signal  refers  to  echoes  of  the 

incident  source  energy  from  the  desired  target. 

for  shallow  targets  the  desired  signal  may  Include  direct  reflections 
from  the  target  and  multiple  reflections  between  the  target  and  the 
antenna.  Since  the  antenna  is  so  near  the  target,  the  antenna  radi- 
ation mechanisms  and  the  target  scattering  mechanisms  may  not  be  dis- 
tinct. 


Noise  and  clutter  are  the  extraneous  signals  that  the  signal- 
processing  unit  is  designed  to  suppress  under  certain  conditions. 

The  scattered  fields  from  the  targets,  both  desired  and  un- 
desired, plus  noise  and  clutter  produce  a signal  at  the  terminals 
of  the  receive  antenna.  In  most  conventional  receivers  noise  is 
reduced  by  the  introduction  of  filters.  For  broad-band  signals  this 
is  not  possible  but  noise  is  reduced  substantially  by  averaging  a 
number  of  received  waveforms.  In  general,  however,  for  this  antenna 
system  even  in  our  local  urban  environment  where  many  strong  inter- 
fering signals  exist,  the  voltage  pulse  caused  by  the  scatterer  is 
clearly  visible  on  the  oscilloscope.  In  the  subsurface  pulse  radar 
system,  due  to  the  low  sensitivity  of  the  antenna  system  to  above- 


*Waveforms  collected  by  the  microcomputer  system  consist  of  128 
samples  in  a time  window  of  25  ns.  The  microcomputer  system  is 
discussed  in  Chapter  VIII. 
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ground  disturbances,  the  noise  level  is  inherently  low.  Furthermore, 
it  was  found  that  a simple  arithmetic  averaging  process  is  an  effec- 
tive means  for  reducing  noise[22,24],  This  means  of  course  that  the 
noise  is  not  target-induced. 

Typical  average  raw  waveforms  are  shown  in  Figure  6.  The 
second  waveform  ( — ) in  each  figure  is  obtained  by  rotating  the 
antenna  by  90°.  If  the  observed  signal  is  caused  by  direct  coupling 
between  the  transmit  and  receive  antenna  (T-R  coupling),  the  signal 
would  not  be  changed  by  this  rotation.  On  the  other  hand  If  the 
signal  is  caused  by  any  external  scatterer,  a polarity  reversal 
is  observed.  These  raw  waveforms  illustrate  the  various  classes  of 
signals.  Figure  6-a  shows  a no-target  waveform  which  is  a received 
waveform  with  no  target  (desired  or  undesired)  present  within  the 
radar  range.  Such  a waveform,  after  averaging,  contains  clutter 
only.  Figure  6-b  shows  a waveform  from  the  desired  mine-like  target. 
In  this  waveform,  the  T-R  coupling  and  the  ground  surface  clutter 
occur  early  in  time  and  because  of  the  shal  ow  target  depth  a certain 
portion  of  the  desired  signal  is  overlapped  by  the  clutter.  It  Is 
these  extraneous  signals  that  various  signal-processing  techniques 
are  designed  to  suppress.  Figure  6-c  shows  a waveform  from  the  brass 
cylinder. 


E . Processed  Waveforms 

Refore  proceeding  to  the  target  identification  algorithms,  a 
certain  amount  of  preprocessing  of  the  data  is  desirable.  This  Is 
essential  here  since  the  goal  is  to  obtain  the  purest  scattering  data 
possible  for  target  characterization.  These  steps  are  made  possible 
by  the  availability  of  a computer  in  the  measuring  system.  The  data 
processing  will  not  all  be  essential  in  a field  system  to  detect  and 
identify  such  targets.  It  is  envisioned  that  any  such  preprocessing 
that  may  be  required  can  be  done  in  a microcomputer  which  will  be  a 
part  of  the  final  system.  The  preprocessing  here  included: 

1.  Arithmetic  averaging:  this  process  forms  the  arithmetic 

average  of  ten  raw  waveforms  at  the  same  antenna 
location  and  orientation. 

2.  Amplitude  shift:  any  dc  drift  in  the  waveforms  is 

corrected  by  adjusting  the  base  line. 

3.  Gating : the  time  regions  before  the  T-R  antenna  coupling 

and  after  the  first  balun  reflection  are  replaced  bv 
a straight  line  at  zero  level.  The  resulting 
“effective"  time  window  (same  for  all  waveforms)  Is 
36.5  ns  wide  ( 186  samples). 

4.  Time  shift:  the  effective  time  window  Is  shifted  to 

the  same  time  region  for  all  waveforms. 


1 
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5,  90°-rotation  difference  technique:  this  process  forms 

the  difference  between  two  average,  shifted  and 
gated  waveforms  from  the  same  antenna  location  but 
with  different  antenna  orientations  of  which  one 
Is  90°  rotation  with  respect  to  the  other. 

6.  Multistation  averaging:  this  process  forms  the 

arithmetic  average  of  all  the  difference  waveforms 
(5.  above)  from  the  antenna  locations  which  present 
identical  relative  geometry  between  antenna  and 
target. 

These  processed  waveforms  are  considered  to  be  relatively 
noise  and  clutter-free.  Furthermore,  they  contain  only  the  back- 
scattered  information  from  the  individual  targets  since  the  ground 
of  the  target  sites  was  undisturbed  (free  of  other  objects)  and 
relatively  smooth.  These  processed  waveforms  are  to  be  used  only  to 
establish  the  parameters  of  Identification  system.  They  are  not  used 
for  the  identification  of  an  unknown  target. 

A set  of  processed  waveforms  and  their  Fast  Fourier  Transforms 
(FFT) [64]  for  the  mine-like  target  at  various  antenna  locations  in  a 
wet  ground  is  shown  in  Figures  7 and  B.  Typical  waveforms  for  the 
other  targets  are  shown  in  Figures  9 and  10. 

The  following  important  generalizations  with  regard  to  these 
waveforms  are  made: 

1.  All  time-domain  waveforms  exhibit  transient  behavior 
in  the  late-time  region  where  only  the  natural  response 
of  the  target  exists.  This  transient  behavior  is  of 
prime  importance,  and,  as  will  be  shown  in  Chapter  III, 
dictates  the  characterization  and  identification  method 
for  these  taryets, 

2.  The  strong  peaks  in  the  FFT's  of  the  waveforms  Indicate 
the  possible  existence  of  resonance  behavior  In  the  back- 
scattered  waveforms.  These  peaks  may  be  a good  approxi- 
mate measure  of  the  imaginary  parts  of  the  complex  reson- 
ances of  the  targets  in  situ.  Note  that  while  the  time- 
domain  waveforms  from  different  antenna  locations  over 

the  mine-like  target  change  noticeably,  the  locations  of 
the  strong  peaks  in  their  FFT's  stay  relatively  unchanged 
(see  the  vertical  dotted  lines  in  Figure  8),  Indicating 
the  complex  natural  resonances  of  the  target  are  excita- 
tion invariant.  This  was  anticipated  and  Is  the  most 
attractive  feature  of  the  target  characterization  scheme. 


ANTENNA  LOCATION  « CENTER 
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Figure  7.  Processed  waveforms  from  the  mine-like  target  at 
different  antenna  locations  in  wet  ground. 


Figure  8.  FFT  of  fhe  processed  waveforms  from  the  mine-like  target. 


Processed  waveforms  from  the  other 
subsurface  targets. 


WOOD  BOARD 

Figure  10.  FFT  of  tne  processed  waveforms  from  the  other  subsurface  targets 


3. 


For  the  waveforms  given  in  Figures  7 and  9 the  signal 
level  of  the  cylinder  is  the  highest  while  the  signal 
level  of  the  sphere  is  the  lowest.  As  a reasonable  quanti- 
tative parameter  for  comparison  of  signal  and  clutter 
levels,  the  following  definition  of  signal -to-clutter  ratio 
( S/C ) was  used  in  this  study. 


s a £t__  „ em~ent 


where 

Ej  is  the  energy  of  the  target  signal, 

Em  is  the  energy  of  the  measured  waveform,  and 

Ent  is  the  statistical  mean  energy  of  the  ensemble  of 
clutter  or  no-target  waveforms  used  to  estimate 
Em j . Tn  this  study,  we  considered  only  single- 
target situations,  hence  a collection  of  51  no- 
target measurements  taken  at  various  locations  in 
the  vicinity  of  the  target  site  was  used  as  the 
ensemble  of  clutter  waveforms. 


The  Energy  of  a waveform  was  defined  and  estimated  as  follows. 
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whore  r(t)  is  the  waveform  under  consideration  and  ts,tp  are  the  start 
and  stop-time  of  the  interval  of  interest.  In  this  study  ts  was  taken 
to  be  the  time  at  which  the  absolute  maximum  of  the  waveform  occurred 
and  t,,  was  taken  to  be  the-  time  at  which  the  balun  reflection  occurred. 
The  reasons  for  these  choices  of  ts  and  te  are  given  in  later  chapters. 


The  mean  clutter  level  was  evaluated  as 
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Using  the  above  definitions,  the  signal -to-clutter  ratio  at 
various  antenna  locations  over  the  5cm  deep  targets  were  evaluated 
and  are  given  in  Figure  11.  It  was  found  that  the  signal -to-clutter 
ratio  for  the  mine-like  target  ranged  from  0.21  to  3.50  depending  on 
the  antenna  location  and  orientation*.  The  brass  cylinder  and  the 
wood  board  targets  had  the  highest  and  lowest  signal -to-clutter 
ratios,  respectively. 

The  three  items  mentioned  above:  the  transient  behavior,  the 
complex  natural  resonances  and  the  signal  level  of  the  backscattered 
waveforms  are  of  prime  importance  in  subsurface  target  identifica- 
tion. Their  dependence  on  changing  ground  conoitions  complicates  the 
identification  process.  Any  practical  subsurface  target  identifica-  '• 

tion  algorithm  must  be  able  to  adapt  to  this  changing  condition  and  j 

identify  the  desired  target  in  a wide  range  of,  if  not  all,  ground  j 

conditions.  A comparison  between  the  brass  cylinder  waveform  in  dry  • 

and  icy  ground  aiven  in  Figure  12  clearly  shows  the  effects  of  chang-  • 

ing  ground  condition.  Although  both  waveforms  exhibit  similar  trans- 
ient behavior,  the  time  intervals  between  the  zero  crossings  are 
different,  indicating  a shift  in  the  locations  of  the  target  reson- 
ances. The  amplitudes  of  the  two  waveforms  are  also  different. 

The  FFT's  of  the  backscattered  waveforms  indicate  that  the  ; 

strongest  frequency  concentration  is  at  about  70  MHz,  furthermore,  * 

almost  all  signal  energy  is  in  the  0-500  MHz  frequency  band.  This 
"system  bandwidth"  is  of  great  significance  in  subsurface  target  i 

identification  for  it  dictates  the  number  of  target  resonances  and  ' 

the  magnitude  of  the  correspond’ ng  residues  in  the  backscattered 
waveforms. 

In  the  next  chapter  we  attempt  to  characterize  the  various  rub-  ■ 

surface  targets  by  approximating  their  processed  backscattered  wave- 
forms with  d finite  complex  exponential  series  with  the  complex  expon- 
ents being  the  complex  natural  resonances  of  the  targets.  A method 
for  extracting  these  resonances  directly  from  the  time-domain  wave- 
forms will  be  presented.  Results  from  the  application  of  this  method 
are  given. 


*S/C  depends  also  on  the  ground  condition.  The  estimates  given  in 
Figure  11  was  based  on  i set  of  measurements  obtained  in  a relative- 
ly dry  ground  condition  over  a time  period  of  several  days.  Further- 
more, 5/C  at  symmetric  locations  are  assumed  equal.  S/C?0  when 
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i'j.jro  11.  Signal -to-clutter  ratio  estimates  of  the 
waveforms  from  the  subsurface  targets 
at  different  antenna  locations. 


CHAPTER  III 

CHARACTERIZATION  OF  SUBSURFACE  TARGETS  BY 
THEIR  COMPLEX  NATURAL  RESONANCES 


A.  Objectives 

This  chapter  summarizes  a study  of  the  characterization  of  sub- 
surface targets  by  their  complex  natural  resonances.  A method  for 
extracting  the  resonances  from  the  processed  backscattered  waveforms 
is  derived  for  completeness.  The  method  is  known  as  Prony's  method 
[12-15]  and,  when  applied  to  measured  data,  It  is  extremely  sensitive 
to  the  values  of  its  parameters.  An  approach  to  solve  certain  of 
these  problems  will  be  presented.  Prony's  method  is  applied  to  the 
processed  backscattered  waveforms  and  an  analysis  of  the  resulting 
resonances  is  focused  on  the  following: 

1.  The  excitation  invariance  of  the  complex  natural 
resonances  of  the  individual  targets  with  respect  to 
antenna  location. 

2.  The  degree  of  distinction  between  the  complex  natural 
resonances  of  the  different  targets. 

3.  The  effects  of  changing  ground  condition  on  the 
location  of  the  complex  natural  resonances. 


B . Comp 1 ex  Natural  Resonances 

The  concept  of  using  complex  natural  resonances  for  target 
characteri ration  is  developed  from  the  fact  that  all  finite-size 
objects  have  resonances  that  depend  on  their  physical  characteristics 
such  as  size,  shape  and  composition  as  well  as  the  medium  surrounding 
the  object.  These  resonances,  however,  are  independent  of  the  excita- 
tion[40].  As  a useful  but  inexact  analogy,  in  circuit  theory,  the 
form  of  the  transient  response  of  a lumped  linear  circuit  may  be 
determined  from  the  knowledge  of  the  resonances  and  the  corresponding 
residues  of  the  response  function  in  the  complex  frequency  plane. 

The  actual  transient  response  of  the  circuit  is  then  simply  a summa- 
tion of  all  the  residues  multiplied  by  the  inverse  transforms  of  the 
resonances.  In  1965,  Kcnnaugh  and  Moffat t[41 3 generalized  the  impulse 
response  concept  to  include  the  distributed  parameter  scattering  pro- 
blems and  suggested  that  a lumped  circuit  representation,  at  low 
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frequencies  or  long  time,  was  possible.  Later,  similar  and  more  for- 
mal representations  have  been  designated  as  the  Singularity  Expansion 
Method  (SEM)[7-8].  This  hypothesis  is  generally  supported  by  the  fact 
that  typical  transient  response  waveforms,  such  as  those  shown  in 
Chapter  II,  appear  to  be  dominated  by  a few  exponentially  damped  sinu- 
soids. Based  on  this  concept,  a subsurface  target  can  be  character- 
ized by  a set  of  complex  natural  resonances  which  is  Independent  of 
the  location  and  orientation  of  the  crossed-dipole  antenna.  These 
resonances.,  however,  are  dependent  on  ground  condition.  Such  a char- 
acterization is  attractive  for  it  catalogs  a target  by  a small  set 
of  complex  numbers. 

The  backscattered  waveforms  from  the  subsurface  targets  received 
by  the  pulse  radar  system  are  good  approx imat i ons  to  the  impulse  re- 
sponses of  the  targets.  Furthermore,  they  appear  to  be  dominated  by 
a few  exponentially  damped  sinusoids,  and  thus  can  be  represented  as 

N s-t  t 

r(t)  = l a^  e (4) 

n=l 

where  r(t)  is  the  received  transient  waveform,  s^'s  are  the  complex 
resonant  frequencies  or  pole  locations  in  the  complex  frequency  plane. 
These  have,  by  common  usage  in  this  representation,  become  designated 
as  complex  resonances  or  more  simply  as  resonances.  These  various 
terms  will  be  used  for  in  this  document,  ai's  are  the  correspond- 
ing residues  and  ii  is  the  number  of  complex  resonances  within  the  fre- 
quency band  of  the  radar  system.  The  corresponding  expression  in  the 
complex  frequency  domain  is 


.twin  = >;  fs,iT 


where  \.[  ] is  the?  Laplace  transform  operator[42]  and  s is  the 

complex  frequency.  Note  that  the  resonances  are  hot  dependent  on 
antenna  location  and  orientation  however,  the  residues  are. 

In  order  to  exploit  Equation  (4),  it  is  necessary  to  first 
determine  the  values  of  the  complex  natural  resonances  of  the  targets. 
The  method  used  here  extracts  the  resonances  of  a target  directly 
from  its  transient  response.  This  method  is  known  as  Prony's  method, 
which  was  first  derived  by  Prony  in  1795[12],  and  was  later  suggested 
by  Van  Blaricum,  et  al . , for  extracting  the  pole  singularities  of 
transient  waveforms  in  1975[14]. 
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Derivation  of  Prony 1 s Method 

In  discrete  form,  Equation  (4)  can  be  written  as 


N s-jKTg 

r(KTg)  = l fli  e K = 0,1,2- • • 

i = l 
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(6) 


where  K is  the  sampling  index  and  Tg  is  the  basic  hardware  sampling 
period  of  our  measurement  system  (see  Chapter  II).  For  an  exact 
solution  of  the  2N  unknowns  and  s-| , we  can  set  up  2N  (nonlinear) 
equations  by  using  2N  sample  values  of  r(KTg).  Prony's  method  uses 
2N  uniform  samples,  and 

N s.j  nT 

r(nT)  = V e ; n = 0, 1 ,2- • • ,M  » 2N-1  (7) 

i = l 


where  T,  the  Prony  interval,  is  the  Interval  between  the  samples  used 
along  the  waveform.  In  general,  N itself  also  represents  an  unknown 
which  is  usually  fixed  by  a trial  and  error  process.*  If  no  waveform 
interpolation  is  exercised,  T is  equal  to  integer  multiples  of  Tg. 
Writing  out  Equation  (7),  we  have  2N  equations 


ro  ' al  + a2  + aN 

n = *1*1  + a2z2  •••  + aN2N 

l ? 

r2  = a-j  i\  + a2z2+  “ ‘ + 

*m  , M , . M 

rM  = tUl  1 d2z2  + + l1NzN 


(8) 


where 
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*N  represents  the  number  of  target  resonances  which  are  excited  by 
interrogating  frequencies  within  the  "system  bandwidth". 
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equation  (H)  is  a set  of  nonlinear  equations  in  the  z^’s.  Let  z-j , z2 
be  the  roots  of  the  algebraic  equation 

n0  + + ugZ*"  + ~ 0 (9) 

so  that  the  left  hand  side  of  Equation  (9)  is  equal  to  the  product 


( Z-Z1 ) (Z-Z2 ) •••  (z-zN)  = 0 


(10) 


that  is, 

N N 

);  = n (z*Z-j)  = 0 • (11) 

m=0  1*1 

Titus,  if  we  can  evaluate  <tm,  then  z^  can  be  obtained  by  a simple 
factorization  of  an  Nth  degree  polynomial,  To  solve  for  am,  we 
obtain  from  Equations  (7)  and  (8) 


N 

n ''inrK+iii 
in-0 


m=0 


; K = 0,1 ,2* • *M-N 


Interchanging  the  order  of  the  summation  yields 


N 

,lllirK+m 

in=Q 


From  Equation  (11),  we  see  that  the  summation  inside  the  parenthesis 
of  the  above  equation  is  zero,  thus,  we  arrive  at  the  desired  linear 
homogeneous  difference  equation 

N 

nmrK+m  - ° 1 K * 0,1,2- -M-N  . (12) 

m=0 


Thus,  the  sample  values  of  r(t)  satisfy  an  Nth  order  linear  homogene- 
ous difference  equation.  This  difference  equation  is  commonly  refer- 
red to  as  the  Prony  difference  equation. 

The  Prony  difference  equation  is  linear  and  homogeneous,  and 
can  be  used  to  solve  for  the  N+l  coefficients,  l.e.,  um‘s.  In  the 
classical  Prony's  method,  these  coefficients  are  obtained  by  setting 
um=1  and  solving  the  resulting  matrix  equation  by  matrix  inversion 
tnat  is, 
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AB  - C 


where 


f’o  rl 


A - • 


[Vn  rM-N+1  rM-N+2  Vll 


'M 

B = . 


and  C = - 


Note  that,  for  M-2N-1,  A is  a square  symmetric  circulant  matrix  and 
is  readily  invertable.  Standard  computer  routines  such  as  GLLGL43J 
can  be  used  to  do  the  matrix  inversion.  Once  the  nm  s are  deter- 
mined, the  next  step  is  to  solve  for  the  N values  of  z\  . These  z^'s 
are  obtained  by  finding  the  roots  of  Equation  (11).  The  N roots  are 
complex  number's  and  because  r(t)  is  real,  these  complex  numbers  appear 
in  complex  conjugate  pairs.  The  polynomial  root  finding  process  can 
be  easily  performed  by  using  standard  routines  such  as  Mul ler[44 ,45J . 

It  is  now  trivial  to  obtain  the  poles  s-j . Since  the  toots  of 
Equation  ('ll)  were  defined  by  Equation  (8),  the  poles  are  simply 


f t'.n(Zj) 


The  final  step  in  Prony's  method  is  to  determine  the  value  of 
t|K.  residues  uj . To  do  this,  we  simply  solve  the  matrix  equation 
embodied  in  Equation  (8).  In  matrix  form  this  set  of  equations  is 
written  as 


where 
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where  now  the  only  unknowns  are  the  elements  of  the  residue  matlrx  E, 

The  above  derivation  of  Prony's  method  Is  valid  only  when  all 
natural  resonances  present  are  simple  poles.  For  multiple-order 
poles,  a slight  modification  Is  necessary  in  solving  for  the  residues. 
The  derivation  of  Prony's  method  for  multiple-order  poles  is  given  In 
Appendix  A for  completeness.  However,  In  this  study  we  have  not  found 
It  necessary  to  postulate  multiple-order  poles. 


f 


In  summary i Prony's  method  solves  for  the  complex  natural  reso-  j 

nances  (poles)  and  the  corresponding  residues  associated  with  the  back-  \ 

scattered  time-domain  waveforms  from  a system  of  nonlinear  equations 
(Equation  (8))  by  breaking  It  down  into  three  simple  steps: 

(1)  Solve  for  the  values  of  um's  of  the  linear  Equation  (13) 
by  matrix  inversion. 

(2)  Solve  for  the  poles  by  factoring  the  polynomial  of 
Equation  (11). 

; 

(3)  Solve  for  the  residues  from  the  linear  Equation  (15) 
by  matrix  inversion. 

The  derivation  of  Prony's  method  Is  simple  enough.  However,  Its 
application  to  the  measured  backscattered  waveforms  Is  a much  more  f 

complicated  process.  The  following  section  outlines  some  of  the  | 

difficulties.  f 


D. 


Clutter  and /or  Noise  in 

Prony's  Method 


Prony's  method  has  been  found  to  be  extremely  sensitive  to 
clutter  and/or  noise.  Its  ability  to  extract  the  complex  natural 
resonances  of  a waveform  accurately  Is  severely  inhibited  by  the  pre- 
sence of  clutter  and  no1se[46-48] . Since  Prony's  method  Is  an  Inter- 
polation process  (otherwise  referred  to  as  curve  fitting),  In  the  pre- 
sence of  clutter  and/or  noise,  It  will  give  a set  of  poles  which  fit 
the  noisy  transient  response  but  will  not  necessarily  represent  the 
complex  natural  resonances  of  the  target.  Various  signal-processing 
techniques  have  been  applied  to  reduce  the  effect  of  clutter  and/or 
noise  In  Prony's  inethod[46-48],  with  the  most  commonly  used  being  the 
least-square  error  technique.  With  It,  Equation  (13)  In  the  previous 
section  is  solved  in  the  least-square  sense.  In  this  case,  M samples 
are  used  in  lieu  of  2N  samples  where  M>2N.  Thus,  the  matrix  A becomes 
rectangular,  and  Equation  (13)  Is  solved  by  the  pseudo- Inverse  techni- 
que 


AtAB  ■ AtC 


(16) 


or 


4>B  = D 


(17) 


where 


and 


AtA 


0 ■ AtC 


Since  <!■  Is  the  signal  covariance  matrix,  It  Is  real,  symmetric 
and  positive  definite,  and  is  thus  readily  invertible  to  yield  the 
value  of  n|n'  s. 

A second  technique  applied  to  reduce  noise  In  Prony's  method 
was  brougnt  about  by  the  observation  that,  In  solving  for  the  N+l 
.im's  in  the  M homogeneous  Equations  (13),  we  can,  Instead  of  setting 
<»N=1,  require  that  the  Euclidean  norm  of  the  a vector  be  1 , i.e. , 


N 

l 

m=0 


(10) 


Such  an  approach  leads  to  the  eigenvalue  method[48 ,49,52] . 


33 


Instead  of  setting  the  loadiruj  coefficient  n jg  1 , we  can  of 

course  set  any  of  the  N+l  coefficients  to  1 in  solving  for  the  um's 
of  Equation  (10).  Such  a constraint  leads  to  the  interpolation 
version  of  Prony's  method[51 ,52]. 

The  classical  Prony's  method,  the  eigenvalue  method  and  the 
interpolation  version  of  Prony's  method  were  all  considered  in  this 
study.  For  completeness,  derivations  of  the  eigenvalue  method  and 
the  interpolation  version  of  Prony's  method  are  given  In  Appendix  B. 

Numerous  other  signal -processing  techniques  have  been  applied 
[46-47],  thus  far,  however,  no  completely  satisfactory  result  has 
been  reported  using  measured  data.  In  the  following  section  a syste- 
matic procedure  that  Is  giving  good  results  for  the  present  data  Is 
outlined.  This  procedure  was  used  in  extracting  the  complex  natural 
resonances  of  the  processed  backscattered  waveforms  In  Chapter  II  and 
yielded  our  best  results  to  date.  As  we  will  see,  the  procedure  does 
indeed  provide  satisfactory  target  separation. 

The  problem  Is  really  one  of  linear  prediction  and  Is  prominent 
in  many  diverse  d 1 sc i pi  1 nes[53] . No  general  solution  has  yet  emerged 
and  the  acceptability  of  a given  method  really  depends  upon  the  appli- 
cation. It  is  interesting  however  that  minimizing  the  total  squared 
error  (actual  vs  estimate)  in  some  sense  Is  a common  starting  point. 
Yet  this  in  no  way  optimizes  the  resonance  locations.  For  our  purpose 
because  the  resonance  locations  need  not  be  found  In  real  time,  pre- 
sent methods  are  adequate  If  not  completely  satisfying.  A real  need 
is  to  compile  and  translate  all  of  the  various  methods  already  being 
used.  Such  a tutorial  unified  review  would  be  Invaluable  but  is  be- 
yond the  scope  of  the  present  effort.  We  can  of  course  easily  Incor- 
porate any  new  techniques  Into  the  identification  scheme. 


E.  Applying  Prony's  Method  to  the  Processed 

Measured  B cc k scat ter eB  Waveforms 

In  applying  Prony's  method,  one  approach  Is  to  pre-determlne 
the  following  parameters: 

1.  N,  the  number  of  poles  to  be  extracted  from  the  waveform. 
Van  B1 ar1cum[14,15,47]  suggested  a method  which  relies  on 
the  fact  that  the  (N+l)th  eigenvalue  of  the  matrix  should 
equal  n2,  the  variance  of  the  additive  gaussian  stationary 
and  uncorrelated  noise.  Such  a method  does  not  seem 
practical  for  our  measured  data  In  which  the  clutter  seems 
to  be  nonstationary  (transient). 

2.  T,  the  Prony  interval.  Obviously,  undersampling  (T  too 
large)  will  almost  surely  bring  aliased  results.  It  was 
also  found  that  oversampling  produces  extraneous  high 
frequency  poles. 
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3.  ts,  te  the  start  and  stop  time  of  the  fitting  interval. 
ts  must  lie  in  the  time  region  where  the  forced  response 

portion  of  the  backscattered  waveform  has  ended  and  te 

must  lie  In  the  region  where  clutter/noise  effects  are  , 

not  dominant.  I 

4.  Mi  the  number  of  sample  points  used  in  the  fitting  process. 

M determines  the  amount  of  overspecification  on  the  system 

of  Equations  (13).  ■' 

In  the  presence  of  clutter  and/or  noise,  It  was  found  that  the 
accuracy  of  the  extracted  resonances  was  found  to  be  extremely  sensi-  i 

tlve  to  the  values  of  the  above  five  parameters[48].  For  the  method 
to  yield  an  accurate  solution,  we  have  to  find  the  "right"  set  of 
values  for  these  parameters.  The  approach  developed  In  this  study 
is  to  vary  these  parameters  (over  a reasonable  range)  and  assume  that 
the  "right"  values  of  the  parameters  corresponding  to  the  "desired" 
resonances  are  those  that  allow  the  closest  approximation  to  the  mea- 
sured waveform  in  the  time  domain.  That  Is,  a calculated  waveform  Is 
developed  from  the  resonances  and  residues  found,  and  this  waveform  is 
compared  to  the  original  waveform  point  by  point  over  the  fitting 
Interval  [tc,te]  and  the  iota!  squared  error  found.  The  solution  j 

which  affords  the  smallest  total  squared  error  is  considered  to  be  the 
"desired"  solution.  This  approach  is  used  to  find  the  coefficients  of 
the  difference  equation  but  once  done  for  a target  needs  not  be  re-  ! j 

peated  for  target-separation  purposes.  It  Is  conceivable  that  such  a 
searching  procedure  can  be  lengthy.  Furthermore,  ranges  of  the  para- 
meters are  dependent  on  the  waveform  being  processed.  However,  with  i 

a little  experience,  one  can  usually  minimize  them.  The  ranges  of 
these  parameters  in  this  sutdy  were  fixed  as  follows: 

1.  The  number  of  "significant"  peaks  In  the  Fast  Fourier 
Transform  of  the  waveforms  Is  usually  a good  measure  of 
N and  since  the  number  of  "significant"  peaks  is  between 
2 and  7 in  all  waveforms  considered,  the  range  of  N was 
chosen  to  be  from  4 to  14  (we  assumed  that  one  peak 
corresponded  to  ut  most  two  poles). 

2.  Shannon's  sampling  theorem  constrains  the  maximum  value 
of  T,  while  the  bandwidth  of  the  radar  system  (<500  MHz) 
constrains  the  minimum  value  of  T.*  The  values  of  T were 
chosen  to  be  3Ty,  5TB  •••  IOTm  corresponding  to  minimum 
and  maximum  Nyquist  ,requenc1es  of  256  MHz  and  768  MHz, 
respecti rely. 


♦Frequencies  beyond  the  system  bandwidth  contains  noise  only  unless 
we  attempt  some  spectral  estimation  techniques,  The  approach  may 
be  worthy  of  study  when  studying  "deeper"  targets. 
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3.  The  values  of  ts  were  chosen  to  be  tn,ax,(tmax+2Tg), 

(tmax+4Ts)  (tmax+8TB)  where  ^lax  is  the  time  at 
which  the  absolute  maximum  of  the  waveform  occurs.  It 
was  found  that  such  a choice  ensured  a decaying  nature 
in  almost  all  the  waveforms  considered,  t was  chosen  to 
to  be  the  time  at  which  the  balun  reflection  occurs 
(see  Chapter  II). 

4.  Since  a vastly  oversized  M would  result  In  an  unstable 
solution  from  Prony's  method[13,15,48],  the  values  of 
M were  chosen  to  be  2N,  3N,  •••  6N. 

Each  set  of  values  of  the  five  parameters  (N,T,ts,te,M)  will 
give  a set  of  complex  resonances  when  Prony's  method  is  applied  to  a 
waveform.  This  set  of  complex  resonances  maximizes  the  fit  between 
the  approximated  waveform  ra(t)  and  the  measured  waveform  r(t)  in  the 
Interval  [ts,ts+(M-l )T]  with  the  sampling  interval  of  T.  For  all 
complex  resonances  resulting  from  all  possible  sets  of  (N,T,ts,te,M) 

In  the  chosen  range,  the  "desired"  set  of  complex  resonances  is  e 
chosen  to  be  the  one  which  minimizes  the  total  normalized  point-by- 
point  squared  error  r over  the  error-calculating  Interval.  The  error 
r , which  will  be  henceforth  referred  to  as  mean-square  error,  is  defined 

C - (j.  (r(t)-rA(t))2jy^  (r2(t)+rj-(t)))  ; t-1TB  (19) 

In  this  study,  r(t)  was  taken  to  be  the  processed  measured  wave- 
form and  the  approximated  waveform  r.(t)  was  generated  via  the  method 
of  linear  predict1on[53] , where 

N -u 

rA(t+V)  H l ~vl~  rM(t+"lT)  . (20) 

m=0  in 
mtmQ 


In  Equation  (20),  r^(t)  and  r^t)  are  the  approximated  and  the 
processed  measured  waveform,  respectively,  the  um‘s  are  the  difference 
equation  coefficients  obtained  from  the  Prony's  method,  and  m0  is  an 
Index  chosen  for  suppression  of  clutter  and  noise  effects.  In  this 
study,  m0  was  chosen  to  be  the  coefficient  of  maximum  magnitude. 

With  Equation  (20),  the  mean-square  error  r.  can  be  expressed  as 
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(21) 


From  Equations  (12)  and  (21),  we  note  that  the  mean-square  error  ■. 
is  zero  when  the  measured  waveform  is  free  of  clutter/noise  and  is 
perfectly  characterized  b.y  Equation  (4).  The  error  should  be  small 
when  rjvjt t ) is  closely  approximated  by  Equation  (4). 

The  choice  of  the  above  error  criterion  is  related  to  the  form 
of  the  correlation  coefficient  osen  in  the  target  identification  al- 
gorithm of  Chapter  IV,  in  such  u way  that  minimization  of  c results 
in  the  maximization  of  the  correlation  coefficient. 

With  the  above  search  procedure,  Prony's  method  and  its  varia- 
tions were  applied  to  the  processed  waveforms  to  extract  their  complex 
natural  resonances.  Results  are  shown  in  the  next  section. 


F . The  Ex crae, ted  Resonances  of 
the  Subsurface  Targets 

The  locations  of  extracted  resonances  of  the  mine-like  target 
at  different  antenna  locations  in  icy  ground  are  plotted  in  Figure  13 
(here  only  poles  in  the  upper  left  half  s plane  are  shown,  details 
are  given  in  Appendix  C).  From  Figure  13,  we  make  the  following 
observations  : 


I.  The  extracted  resonances  Lend  to  form  "clusters".  Some 
possible  clusters  are  shown  in  Figure  13.  The  formation 
of  these  clusters  are  based  on  the  obviousness  of  cluster- 
ing of  the  resonances  and  the  known  fact  that  the  accuracy 
in  determining  the  real  part  of  the  extracted  resonance  is 
normally  poor.  A cluster  can  contain  at  most  one  pole 
extracted  from  a waveform.  Poles  with  residues  which  are 
three  orders  down  in  magnitude  compared  to  the  maximum 
residue  are  discarded.  Poles  which  are  remote  from  the 
clustered  groups  are  excluded.  Beyond  an  obvious  weighing 
dictated  by  the  actual  pole  locations  no  real  significance 
should  be  attached  to  the  shape  of  the  closed  contour 
surrounding  each  cluster. 

?,  Only  a small  number  of  clusters  or  resonances  are  present. 
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3.  The  variation  in  the  real  parts  of  the  resonances  within  a 
cluster  is  generally  greater  than  the  variation  in  their 
imaginary  parts.  There  is  at  least  one  more  major  factor, 
besides  the  ever-present  clutter/noise,  that  causes  such 
variations,  namely,  the  target-antenna  interactions.  At. 
the  shallow  depth  of  5 cm,  for  most  antenna  locations 
considered,  the  targets  are  in  the  near  field  of  the 
antennas  for  the  entire  bandwidth  (<  500  MHz)  of  our  radar 
sys  tem . 

4.  An  additional  factor  contributes  to  the  variations  in  the 
extracted  resonances  from  the  mine-like  target,  namely, 
its  complex  structure.  This  target  possesses  the  most 
complex  structure  of  all  targets  considered. 

5.  The  phenomenon  of  certain  resonance(s)  being  weakly  excited 
in  certain  radar  aspects  is  evident.  The  weakly  excited 
resonances  were  not  extracted. 

6.  As  expected,  the  residues  are  aspect  dependent  This 
becomes  evident  by  noting  the  variations  of  v magnitude 
of  the  residues  of  the  poles  in  the  clusters  {See 
Appendix  C). 

7.  The  mean-square  error  c is  small  (<_  .01,  see  Appendix  C) 
in  all  cases  considered,  meaning  that  the  finite  sum  of 
complex  exponentials  fits  the  measured  waveforms  well. 

This  is  a necessary  condition  for  our  identification 
algorithm  whose  correlation  coefficient  is  defined  to  be 
unity  minus  the  mean-square  error  e. 

In  tin's  study,  a subsurface  target  was  characterized  by  the  set  of 
average  extracted  resonances.  Averaging  was  performed  over  all  the 
extracted  resonances  in  each  cluster.  For  the  mine-like  target,  the 
average  extracted  resonances  are  shown  as  solid  dots  in  Figure  13. 
Parameters  such  as  the  variation  from  the  average  of  each  pole  with- 
in the  cluster  is  not  meaningful  because  of  its  causes  which  include, 
besides  the  effects  of  clutter  and  noise,  the  possible  variations  in 
the  pole  excitation  at  the  various  antenna  locations  and  orientations. 
Slight  pole  variations  due  to  the  variations  in  the  antenna  locations 
and  orientations  is  possible  for  the  finite  exponential  sum  represent- 
ation of  the  target's  transient  response  is  only  an  approximation  and 
that  the  targets  considered  are  located  in  the  closed  vicinity  of  the 
radar  system. 

The  extracted  resonances  shown  in  Figure  13  were  obtained  using 
classical  Prony's  method  (i.e.,  a^=l).  Classical  Prony's  method  was 
found  to  extract  poles  with  tighter  clusterings  among  the  results 
given  by  other  methods  under  the  constraint  am=l ,m=0 , 1 • • • ,N[52] . The 
eigenvalue  method  provided  results  similar  to  those  given  by  the 
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Classical  Prony’s  method.  Thus,  no  clear  cut  choice  of  method  was 
dlscernable.  Accordingly,  the  extracted  resonances  shown  In  this 
dissertation  were  the  results  of  either  of  these  two  methods. 

In  order  to  see  the  effects  of  the  changing  ground  condition  on 
the  location  of  the  extracted  resonances,  the  average  resonances  of 
the  mine-like  target  in  different  ground  conditions  are  tabulated  in 
Table  1 and  are  plotted  In  Figure  14.  From  Figure  14,  we  see  that 
there  were  five  (pairs)  extracted  resonances.  The  Imaginary  parts  of 
the  extracted  resonances  were  relatively  Insensitive  to  changes  In 
ground  condition.  This  seems  to  Imply  the  resonances  of  the  mine- 
like target  are  internal  resonances.  This  also  means  that  the  carget 
Identification  scheme  when  applied  to  this  target  will  be  relatively 
Insensitive  to  ground  conditions. 

The  Implications  of  the  fact  that  there  were  five  (pairs)  reson- 
ances extracted  from  the  mine-like  target  waveforms  Is  significant. 

It  means  that  this  target  can  now  be  characterized  by  a finite-order 
system. 

Not  all  the  extracted  resonances  are  related  to  the  scattering 
mechanisms  of  the  mine-like  target.  In  fact,  the  lowest  resonance 
was  found  to  be  the  antenna  resonance  of  the  system.  This  becomes 
particularly  clear  when  we  study  the  resonances  extracted  from  the 
mine-like  target  waveforms  collected  with  a 12m  long  antenna.  3y 
gating  out  the  latertime  portion  of  the  backscattered  waveform  from 
the  mine-like  target  received  by  the  12m  long  antenna,  we  effectively 
eliminated  the  resonance  of  the  antenna  created  by  the  finite  length 
of  the  antenna  arms.  Thus,  poles  extracted  from  these  waveforms  are 
all  target-related.  A typical  such  backscattered  waveform  In  the 
time  domain  Is  shown  In  Figure  15,  The  short  time  window  (compared 
to  the  0.6m  long  antenna  waveforms  of  Figures  7)  of  the  waveform 
indicates  the  absence  of  any  low  frequency  content.  The  average 
extracted  resonances  from  the  mine-like  target  waveforms  received  by 
the  12  in  long  antenna  at  various  locations  arc  shown  In  Figure  1C. 

A quick  comparison  with  the  extracted  resonances  from  the  0.6m  long 
antenna  waveforms  reveals  the  fact  that  the  resonance  with  imaginary 
part  of  approximately  60  MHz  is  the  antenna  resonance.  Thus,  for  the 
system  with  the  0.6m  long  antenna,  four  (pairs)  target  resonances 
were  present  in  the  received  waveforms.  The  antenna  resonance  was 
extracted  from  almost  every  waveform  of  all  targets  considered. 

In  contrast  to  the  case  of  the  plastic  mine-like  target,  the 
brass  cylinder  was  found  to  possess  external  resonances.  Table  2 
lists  the  average  extracted  resonances  of  the  brass  cylinder  in 
various  ground  conditions.  Locations  of  these  resonances  are  also 
plotted  in  Figure  17.  From  Figure  17  we  see  the  following  effects 
of  the  changing  ground  condition  on  the  extracted  resonances  of  the 
brass  cyMnder.  First,  the  antenna  resonance  is  Insensitive  to  changes 
in  the  ground  condition.  This  may  be  attributed  to  the  fact  that  the 
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arms  of  the  crossed-dipole  antenna  were  not  in  electrical  contact  with 
the  ground  surface.  Second,  the  imaginary  parts  of  ttie  three  higher 
order*  resonances  increased  significantly  when  the  ground  changes  from 
dry  to  icy.  This  is  to  be  expected,  because  the  resonances  of  the 
brass  cylinder  are  external  resonances.  The  Increase  in  the  Imaginary 
parts  of  these  external  resonances  indicated  a decrease  in  the  value 
of  the  dielectric  constant  of  the  ground.  Third,  the  real  part  of 
the  three  higher-order  resonances  generally  decreased  as  the  ground 
changed  from  dry  to  icy,  indicating  that  Icy  ground  in  this  case  was 
more  lossy.  The  increase  in  loss  seemed  to  be  the  reason  for  the 
absence  of  the  real  cylinder  pole  in  icy  ground. 

The  average  extracted  resonances  of  the  aluminum  sphere,  copper 
sheet  and  wood  board  are  tabulated  in  Table  3.  From  Table  3 we  see 
that  the  antenna  resonance  is  present  in  the  waveforms  of  all  targets. 
Note  that  the  extracted  resonances  of  the  five  targets  considered  lay 
in  the  same  general  region  of  the  complex  frequency  plane  and  are 
only  marginally  separated.  Such  is  expected  to  some  extent  because 
all  targets  considered  have  (again  marginally)  similar  sizes.  Such 
marginal  level  of  distinction  in  the  poles  is  expected  to  present 
difficult  tests  for  the  identification  algorithm.  Furthermore  the 
number  of  resonances  for  the  various  targets  are  also  close  (4  to  5 
pairs);  this  further  tests  our  identification  method. 

The  location  cf  the  target  resonances  are  related  to  the  scat- 
tering mechanisms  of  the  target.  For  subsurface  targets,  these  rela- 
tionships are  complicated  by  the  presence  of  the  air-ground  interface, 
the  ground  condition  and  the  characteristics  of  the  transient  antenna 
system.  For  shallow  targets  the  near-field  effects  and  the  target- 
antenna  interactions  further  complicate  the  picture.  In  this  disser- 
tation, we  do  not  Intend  to  explore  these  relationships.  Instead,  we 
proceed  to  use  the  extracted  resonances  for  identification  of  the 
various  subsurface  targets.  In  the  next  chapter,  a basic  Identifica- 
tion algorithm  will  be  given  and  identification  results  using  the 
extracted  resonances  as  the  discriminants  will  be  presented. 


*0rder  here  denotes  Increasing  imaginary  part. 
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Flquri1  14.  Locations  of  the  average  extracted  resonances 
of  the  mine-like  target  In  different  ground 
conditions. 
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‘ical  backscattered  waveform  received  by  the  12m 
antenna  from  the  mine-like  target. 


a l 10°  NEPERS /s  ) 

Figure  16.  Average  extracted  resonances  from  the 
mine-like  target  waveforms  received 
by  the  12m  long  antenna. 
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CHAPTER  IV 

THE  PREDICTOR-CORRELATOR  IDENTIFIER 


A . Objectives 

This  chapter  summarizes  the  target  identification  procedure 
based  on  the  predictor-correUtor  Identification  method[40,9*li], 
A detailed  analysis  on  the  predictor-correlator  is  presented. 
Identification  performances  based  on  real  radar  measurements  are 
y iven. 


Processing  for  target  identification  using  the  predi ctor-cor- 
relator  consists  of  comparing  a measured  waveform  from  an  unknown 
target  with  a calculated  waveform  produced  using  the  resonances  of 
a known  desired  target.  The  procedure  is  as  follows: 

1.  Preprocessing:  The  preprocessor  attempts  to  suppress 

clutter  and  noise.  For  "single-look"*  identification, 

the  following  preprocessing  steps  are  taken: 

a.  Arithmetic  Averaging  (see  Chapter  II), 

b.  Amplitude  Shift  (see  Chapter  II). 

c.  Time  Shift  (see  Chapter  II). 

d.  90'’-rotation  Difference  (see  Chapter  II). 

e.  Filtering:  It  was  found  that  with  the  system  under 
discussion,  almost  all  of  the  target  signal  energy 
resided  in  the  0-500  MHz  region  (see  Chapter  II), 
thus,  to  suppress  out-of-band  clutter  and  noise  a 
low-pass  trapezodal  filter  with  the  transfer  function 
shown  in  Figure  18  v;as  Inserted  Into  the  preprocess- 
ing unit.  Various  critical  frequencies  were  tested, 
the  ones  shown  in  Figure  18  yielded  the  best  identifi- 
cation performance. 


* i.e..  Identification  based  on  a single  radar  observation 
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Figure  18.  Transfer  function  of  the  pre-processor  filter 
used  in  the  0 . 6m  Innq  antenna  system  for 
target  Identification. 


For  design  flexibility,  the  filtering  processes  were 
performed  in  the  frequency  domain  with  the  Fast  Fourier 
Transform  (FFT)  package  available  In  the  digital  computer 
system  library.  The  only  goal  of  the  filtering  operation 
is  to  remove  the  out-of-band  frequency  contents  without 
adding  any  extraneous  frequencies  to  the  spectrum  of  the 
target  response.  The  trapezoidal  filter  structure  is  a 
digital  filter  and  was  used  here  because  of  Its  linear 
phase  characteristics  and  the  availability  of  the  FFT 
package  in  the  digital  computer.  The  FFT  package  allows 
tremendous  flexibilities  In  digital  filter  design.  It 
Is  important  to  note  that  the  use  of  the  trapezoidal 
filter  here  was  not  meant  to  be  "optimum".  More  extensive 
study  may  well  yield  a better  filter  structure. 

2.  Detection:  The  detector  performs  a screening  operation 
for  the  predictor-correlator  identifier  by  rejecting  as 
undesi red-target  waveforms  those  waveforms  whose  para- 
meters are  not  within  the  desired  ranges.  The  para- 
meters considered  in  this  study  were'- 


ajtui  


a.  Waveform  energy  (as  defined  In  Chapter  II). 

b.  Peak  timing  of  the  waveform. 

c.  Peak  amplitude  of  the  waveform 

The  desired  ranges  of  the  detection  parameters  were  determined 
by  studying  the  ranges  of  these  parameters  of  the  desired- target 
waveforms.  Examples  Illustrating  the  choosing  of  the  detection  para- 
meters are  given  In  Section  F of  this  chapter. 


3.  Prediction:  Predict  the  calculated  waveform  using  the 
Prony  difference, equation,. 

4.  Correlation:  Calculate  the  correlation  coefficient  for 
threshold  identification  by  comparing  the  processed  and 
the  calculated  waveforms.  The  correlation  coefficient  Is 
defined  to  be  unity  minus  the  normalized  total  squared 
error. 

This  approach  assumes  that  the  poles  and  the  difference  equation 
coefficients  for  the  desired  target  have  previously  been  obtained. 
Thus,  only  simple  algebraic  operations  are  Involved  in  calculating  the 
correlation  coefficient  for  an  Identification  decision. 


B.  The  Predictor 

Tho  predictor  generates  the  calculated  waveform.  With  the  com- 
plex natural  resonances  and  a chosen  value  of  T,  the  coefficients 
■«R1's  of  the  Prony  difference  equation  (see  Equation  (12))  can  be 
determined  via  Equation  (11).  Thus,  for  one  value  of  T,  we  can  gen- 
erate a calculated  waveform  by  one  of  the  following  methods; 

1 . One- step  prediction: 

N-l 

rc ( t+N r ,T)  = V -‘*n,rM(t+mT)  (22) 

m=0 


or 


2.  Interpolation: 

N 

rc<t+l"0T’T)  = >’  -■  r(  t+mT ) 

m=0  M 

mfni0 


(23) 
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i 


where 

t a ts+nTg  , n-0,1 ,2-  • • 
and 


(■ 


J 


* 


T«1Tb  , 

rc(t,T)  and  r^. ( t , T ) are  the  calculated  and  processed  measured  wave- 
form, respectively.  ts  is  the  start-time  (see  Chapter  II  and  III). 

The  parameter  m0  is  an  index  chosen  for  suppression  of  clutter  and 
noise  effects.  In  this  study  the  interpolation  method  was  used  and 
m0  was  chosen^ to  be  the  coefficient  of  maximum  magnitude.  By  com- 
paring Equations  (22),* (23)  and  (12)  we  note  that  when  the  measured 
waveform  is  from  the  desired  target,  rc(t)  and  ry(t)  are  equal  and 
hence  perfectly  correlated.  Note  that  one  calculated  waveform  is 
constructed  at  each  chosen  value  of  T and  the  same  measured  waveform 
is  used  to  construct  calculated  waveforms  for  all  chosen  values  of  T. 

The  one-step  prediction  method  was  first  used  by  Hill  who 
applied  it  to  the  detection  and  identification  of  targets  (above 
ground)  near  the  half  space[40].  This  method  was  then  modified  by  ] 

Moffatt,  et  al . [9]  to  become  the  interpolation  method  for  Improved 
performance  in  clutter/noise.  The  predictor  uses  past  values  of  the 
measured  waveforms  only,  while  the  Interpolator  uses  both  past  and 
future  values  and,  together  with  its  normalization  process,  was  found 
to  perform  better  in  the  presence  of  clutter/noise. 


C . The  Correlator 

The  correlation  coefficient,  formed  by  comparing  the  calculated 
and  the  processed  measured  waveform,  is  a function  of  T. 


(24) 


where  the  start-time  ts  is  taken  to  avoid  the  forced  response  of  the 
backscattered  waveform; while  the  stop-time  te  is  taken  to  avoid  the 
low-amplitude  signal  at  the  tail  end  of  the  waveform.  In  this  study 
ts  was  taken  to  be  the  time  the  absolute  maximum  of  the  waveform 
occurs,  while  te  was  taken  to  be  the  time  the  balun  reflection  occurs. 

Equation  (24)  can  be  rewritten  as 


.-(T)  ^ 1 


te-NT+m0T 

V 

te-NT+m0f 


[,c(t*T)-rM(t)]‘ 


)!  (^(t.T)^(t)) 

t-t  +m„T  c M 


s o 


te-NT+m0T 

l 2rc(t.T)rM(t) 

t=ts+,"0i 

te-NT+m0T 

I (r|(t.T)+rg(t)) 

t=ts+m0T 


We  note  that  the  numerator  of  the  quotient  in  Equation  (25)  is  the 
cross-correlation  between  the  calculated  and  processed  measured 
waveforms.  The  nominal  range  of  p(T)  is  -1  < p(T)  < 1,  and 


»•  ( T)  = *1  when  rc(t)  = * rM(t) 


The  quotient  term  in  Equation  (24)  is  identical  to  the  mean- 
square  error  >.  in  Equation  (21)  when  the  approximated  waveform  r.(t) 
in  the  Prony  process  is  equal  to  the  calculated  waveform  r (t).  rtThus, 

w 


p(T)  - 1 - 


The  importance  of  Equation  (27)  lies  in  the  fact  that  it  ties 
the  characterization  and  the  identification  processes  toqether  ‘ . 
such  a way  that  optimization  of  one  process  (e.g.,  minimum  e)  w i . I 
lead  to  the  optimization  of  the  other  (maximum  p(T)).  Therefore, 
the  evaluation  of  optimum  values  of  the  parameters  (N,T,ts,te,M)  in 
the  characterization  process  leads  also  to  optimum  values  in  the 
identification  process.  This  is  rather  significant  in  the  evaluation 
of  „(T). 

From  Equations  (24)  and  (25),  we  see  that  the  interval  in  the 
error  calculation  is  from  ts+m  T to  t^-NT+nuT,  givinq  an  error 
interval  size  of  te-ts-NT.  Thus,  the  size  of  the  error  interval 
decreases  as  T increases. 


D The  Predictor-Correlator  as  a Filter 

In  order  to  understand  the  operation  of  the  predi ctor-correl ator 
identifier,  it  is  helpful  to  consider  the  identifier  as  a linear 
t*me-1nvar1ant  filter. 

From  Equations  (21)  and  (24),  we  see  that,  an  instantaneous 
error  e can  be  defined  as 
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e(t+mQT)  = rM(t+m0T) 


rcU+m0T)  = L 


m=0  a' 


• rM(  t+mT) 


(28) 


Thus,  the  instantaneous  error  e can  be  interpreted  as  the  output  of 
a Finite  impulse  response  (FIR)  filter[ 54-56]  with  the  filter  coef- 
ficients a*  where 
m 


(29) 


hote  that  although  the  filtering  operation  Is  based  on  the  sampling 
interval  T,  the  Instantaneous  error  exists  at  the  sampling  period 
of  Tq,  where  T is  taken  to  be  multiples  of  TR.  In  the  Z domain  [54-56], 
the  instantaneous  error  is  given  by  u 


-N+fti  k-1 
E(Z)  - Z 0 l 
i-0 


1 

IT 


C(Z)Rm(Z) 


(30) 


where 


E(Z)  = Z transform  of  the  Instantaneous  error,  based  on  T 
N a. 


C(Z)  - l 


‘N-m 


m=0  am. 


Z , is  the  Z transform  of  the  filter 
coefficient  sequence. 


Rm(Z)  = Z transform  of  the  processed  measured  waveform 


and 


T - KTb,  1.2.-»  . 

The  product  C(Z)RmCZ)  in  Equation  (30)  is  the  filtering  opera- 
tion in  the  Z domain,  bince  the  Z transform  is  based  on  the  sampling 
interval  T,  thus  this  product  term  alone  results  in  samples  only  at 
the  sampling  intervals  <?f  T in  tiie  time  domain.  The  weighed  sum 
operation  (weighed  by  Z1/^)  fills  the  time  interval  between  these 
samples  with  samples  at  finer  sampling  interval  of  Tq,  and  thus 
assures  a better  error  waveform  resolution.  The  presence  of  the 
factor  Z'Nt|"o  is  to  make  sure  that  the  starting  point  of  the  instan- 
taneous error  be  consistent  with  that  stated  in  the  definition  of 
the  instantaneous  error  for  the  case  of  the  interpolation  method. 

Using  Equations  (11),  we  rewrite  Equation  (30)  as 


E(Z) 


Z 


-N+ri 

o 


(l-z-'zj) 


(31) 
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where  Zj  * e are  the  locations  of  the  desired- target  resonances 
in  the  l plane,  Sj  are  the  complex  resonances  In  the  s plane  and  N 
Is  the  number  of  complex  resonances  of  the  desired  target.  Note 
that  the  locations  of  the  complex  resonances  In  the  Z plane  are 
functions  of  the  sampling  Interval  T. 

From  Equations  (24-31),  we  make  the  following  Important  observa- 
tions: 

1.  The  mean-square  error  is  the  normalized  energy  of 
the  Instantaneous  error  e.  e Is  directly  related  to 
e as  follows. 


I e2(t) 

t 

>;7<7)  + 2 l rjjj(t)  - 2 £ rM( t)e(t) 

where  the  summations  In  Equation  (32)  are  taken  over  the  error 
Interval.  The  lower  and  upper  bounds  of  the  mean-square  error  >■ 
are  given  by  Equation  (33) 


1 

7T7~27 

* /x 


where 


: c 2 


/.'  e (t) 
A t 

x _ . .... 


).  rM<t) 

t 

and  the  summations  are  again  taken  over  the  error  interval. 


(33) 


The  error  bounds  are  plotted  as  a function  of  x In  Figure  19. 
From  Figure  19,  the  following  observations  are  made: 


a.  i-  is  bounded  as  follows  : 

2 -t  o - 

i = 2 when  e(t)  = 2r^(t)  and  (34) 

i : 0 when  e( t ) - 0 . 
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b.  The  lower  bound  of  e Increases  monotonlcally  as  the  energy 
of  the  Instantaneous  error  r'*  Increases,  c approaches 

1 as  e'  approaches  Infinity. 

c.  The  upper  bound  of  e increases  monotonlcally  from  0 to  2 
as  c*  Increases  from  0 to  the  value  of 


At  this  e'  value,  e attains  its  maximum  upper  bound  of  2. 
The  upper  bound  of  e decreases  from  2 to  1 as  e1  increases 
to  Infinity.  Although  the  upper  bound  of  c does  not  In- 
crease monotonlcally  as  Increases  from  0 through  its 
entire  range,  the  lower  bound  does.  Thus,  it  Is  reasonable 
to  say  that  large  c'  generally  means  large  e. 

2.  The  difference  between  the  one-step  prediction  and  the 
interpolation  method  lies  in  the  filter  gain  normaliza- 
tion factor  1/0^.  In  the  one-step  prediction  method, 

aN*l , thus,  the  possibility  of  larger-than-one  filter 
coefficients  am  exists.  Such  filter  coefficients  would 
of  course  amplify  clutter  and  noise  effects.  In  the 
Interpolation  method,  is  the  difference  equation 
coefficient  of  maximum  magnitude.  Thus,  no  filter  coef- 
ficient can  have  value  larger  than  one. 

3.  For  the  desired  target,  Ru(t)  is  given  by  the  l transform 
of  Equation  (4)  and  thus  Equation  (30)  can  be  written  as 


, -N-Hti  k-1 

e(z)  « 7y  l 0 l l k 

1=0 


1 


l»;i  oS\) 


1 /N+mo  M , IT  Q(Z) 

z z 


N 

11 

.1  = 1 
N 


1-0  "1  , j-1 

,w' V 


n-z-’zj) 


O-z-’Zj) 


(35) 


te-NT+m0T 


t=t,+m  T 
s o 


e2(t)/((te-ts)/TB) 


, A 

L ;; 
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where  N-.  Is  the  number  of  complex  resonances  present  In 
the  waveform  and,  N is  the  number  of  all  possible  complex 
resonances  of  the  desired  target  present  within  the  radar 
bandwidth  and  used  for  target  identification  purposes. 
Since  certain  target  resonances  may  not  be  excited  for 
some  radar  locations  and  orientations,  Ni  and  N are  not 
necessarily  equal.  Q(Z)  is  a polynomial  of  order  in 


From  Equation  (33),  we  see  that  the  predictor-correlator  Iden- 
tification process  Is  basically  a pole- removal  process.  Furthermore, 
the  process  of  pole  removal  is  performed  by  putting  zero's  at  the 
locations  of  the  poles. 

In  the  case  of  N-|=N, 

, -N+m.  k-1  i 

E(Z)  « l l Z K Q(Z)  . (36) 

V i-o 

o 


Thus 


e(nTB)  = ~ q(nTg)  - 0 for  nTg>jn  T (37) 

m. 


where 


q(nTB) 


k-1 

l 


1=0 


1 

z * Q(Z) 


where  r~  [ ] is  the  inverse  Z transform  operation.  Since  Q(Z)  Is  a 
polynomial  of  order  N In  Z”1 , q(nTg)  Is  zero  for  nTp  > m0T  [54-56]. 
Thus,  the  instantaneous  error  e(nT[j)  is  zero  for  nTp  > m0T  and  the 
mean-square  error  would  also  be  zero  In  this  time  region.  The 
above  derivation  is  based  on  the  assumption  that  the  resonance  time 
region  of  the  waveform  starts  at  t-0.  Since  the  resonance  region  of 
the  waveforms  considered  in  this  study  is  assumed  to  start  at  t-ts, 
the  mean-square  error  would  be  zero  for  t ^ t +m  T.  This  is  the 
lower  limit  of  the  error-calculation  intervaras0 defined  in  Equations 
(21)  and  (24). 


In  the  case  of  N^<N, 
1 


E(Z) 


-N+m.  k-1 

° l z 

1=0 


N-N,  , 

Q(Z)  u 1 (1-7.  'ZJ 

.1=1  J 


(38) 


57 


Thus 


e(nTg) 


q* (nTB) 


(39) 


where 


q'  (nTg)  = c 


-1 


■N-hn0  k-1 


7 N-Ni  l 

l Z k Q(Z)  n C\-rh,) 
1=0  j=l 


(40) 


Since  Q(Z)  Is  of  degree  Ni  In  Z , thus  q'(nTn)  Is  aqaln  zero  for 
nTBlmoT-  Thus  Pitting  the  number  of  all  possible  extracted  resonances 
of  the  desired  target  In  the  Identifier  assures  a zero  error  and  thus 
a unity  correlation  coefficient  for  the  desired  target  waveforms  from 
all  radar  locations  and  orientations.  If  not  all  possible  extracted 
resonances  are  used  for  target  Identification,  the  possibility  of 
N^N  exists  and  In  this  case, 


E(Z) 


1 


■N+m, 


k 

l z 

1-0 


Q(Z) 

FT  T“ 

n (1-Z  \) 


(41) 


1=1 


The  Instantaneous  error  has  poles  In  Its  7.  transform,  causing 
It  to  be  a waveform  of  Infinite  duration  [54-56].  Thus,  a large  mean- 
square  error  could  occur  for  a desired-target  waveform,  possibly  caus- 
ing the  correlation  coefficient  to  be  much  less  than  one. 

4.  For  an  undesired  target,  the  Instantaneous  error  depends 
on  the  processed  waveform  ry|(t)  as  well  as  the  locations 
of  the  zero's  in  the  identifier.  Small  error  could  result 
only  In  the  cases  where  the  undesi red-target  waveform  Is 
closely  characterized  by  a set  of  complex  resonances  which 
are  approximately  equal  to  the  desired-target  resonances 
(i.e.,  the  zeros  in  the  identifier).  In  the  one-pole  case 
where, 


Rm(Z)  = — 


1 


7T 


1-Z  '1 


C(Z)  « (1-Z. 


-1 


z2) 


(47) 


where  Ru(Z),  The  Z transform  of  the  measured  waveform  is  assumed  to 
have  only  one  pole  at  Z-j  with  residue  1.  C(Z),  the  Identification 

filter  Is  assumed  to  consist  of  a zero  at  Zj>.  The  error  can  be 
expressed  as 


„ ^ . ...In.  . . 
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(44) 


where  ||*||  denotes  complex  magnitude.  Hence  the  total  squared  error 
i ' is  directly  proportional  to  the  magnitude  square  of  the  difference 
between  the  locations  of  the  pole  and  the  zero. 

In  a two-pole  case  where 


Rm(Z) 


1-7.'  Z 


l-Z'^Z* 


c(z)  - (1-Z'1Z2)(1-Z'1Z*2) 


(45) 


where  [ 1*  denotes  complex  conjugation.  The  error  can  be  expressed  as 
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or 
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The  total  square  error  c1  Is  Independent  of  the  first  two  terms  on  the 
right  hand  side  of  Equation  (47).  Furthermore,  the  Inverse  Z trans- 
form of  the  last  two  terms  are  given  by  [54] 
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where  o-j  and  iu  are  the  arguments  of  the  complex  numbers  (Z]-Z2)  and 
Z2,  respectively. 

From  Equations  (47)  and  (48),  we  see  that 


e(t)  --  1—  j |Zt -Z2  t le 

"tn„ 


(t-T) 
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Thus,  the  instantaneous  error  depends  on  the  following  three  items: 

2 2a,  T 2a.,T 

a.  The  amplitude  factor  j|Z,-Z2|'  = e + e c - 

(*'l  +a2  )T 

-2e  cos  (a),  -atp  )T  i where  s?=«  ?+ju?  is  the 

location  of  the  zdro^of  the  filter  In  the  s plane, 

-ii,  T *(‘»n"oi,jT 

b.  The  factor  2e  ||Z2II  a 2e  ' d , 

c.  The  phase  term  w-|T+o2*(w.|+ii)2)T. 

All  of  the  above  three  items  are  dependent  on  the  sampling  Inter- 
val T.  Thus,  in  designing  the  identification  radar,  the  value(s)  of 
T should  be  chosen  to  maximize  t-.1  by  considering  Its  dependence  on  the 
above  three  Items.  In  maximizing  we  note  the  following  guidelines: 

a.  Small  values*  of  T result  in  ||Z,-ZJ|  approaching  zero. 

This  causes  > 1 to  approach  zero.  Thus,  small  T values 
offer  no  identification  capability. 

b.  Large  values**  of  T also  result  in  ||Z,-ZJI  approaching 
zero.  Thus,  large  T values  offer  no  identification 
capabil 1 ty . 

c.  Intel  mediate  value  of  T should  be  chosen  to  satisfy  the 
following  conditions  as  closely  as  possible: 

2u,T  uctpT  (<.i.+a,  )T 

i . e e and  e , approach  1 , 

1i.  (<u^ -uig )T  approaches  > n , 

ill.  (u^+i.i^)  approaches  2n  , 

iv.  To  prevent  severe  aliasing  effects,  we  require 
that  n ■ ui 1 1 and  t.^T  . 


* Small  T in  this  case 


means  that  the  quantity  (ui,-m„)T  approaches  0 
a,T  .ipT  (u,+Up)T 

, e and  e approach  1 . 


i “p  * "‘2  / 1 

and  the  quantities  e , e and  e approach  1. 

2u, T 2a?T 

**Large  T in  this  case  means  that  the  quantities  e , e 
(u, +n?)T 


and  e 


approach  ’0. 
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Conditions  1.  and  11.  maximize  the  amplitude  factor 
1 1 Z-t  -Zo  1 1 • Conditions  1.  and  111.  make  sure  that  the  second 
term  In  equation  (50)  Is  maximized  and  has  the  same  sign  as 
the  first  term.  Conditions  11.  and  111.  require  that 
2iuojjT*h.*  Condition  1.  requires  that  T be  as  small  as 
possible.  These,  toqether  with  condition  1v  Indicate  that  a 
reasonable  compromise  Is  to  choose  T such  that 

Although  the  above  analysis  is  based  on  a two-pole  system, 
for  system  of  more  than  two  polos,  the  above  results  are  expected 
to  be  reasonably  valid  when  u>|  Is  taken  to  be  the  Imaginary  part  of 
the  dominant  pole  of  the  measured  waveform  to  be  identified  and  m2 
Is  taken  to  be  the  imaginary  part  of  the  zero  of  the  Identification 
filter  which  is  closest,  to  the  dominant  pole.  In  practical  target 
Identifications  situations,  the  v.'aveforms  from  the  desired  as  well 
as  the  undesired  targets  almost  always  contain  more  than  two  poles 
(see  Chapter  III).  Furthermore,  the  dominance  of  the  resonances  of 
the  undesired  target  Is  generally  unknown,  Thus,  for  better  Identl 
ficatlon  performance,  It  Is  advisable  to  base  the  Identification 
decision-making  process  on  more  than  one  value  of  T.  Based  on  the 
analysis  In  this  section,  the  smallest  "effective"  value  of  T Is  In 
the  neighborhood  of  Tfj  where 


h__  „ J . 

“M 


(51) 


and  i.iM=2h f^  Is  the  maximum  Imaginary  part  of  the  desired-target 
resonances.  In  applying  the  criterion  given  by  Equation  (51)  to 
choose  the  values  of  T for  target  Identification,  to;  word  neighbor- 
hood  should  be  emphasized  and  taken  to  mean  that,  although  tFe  sTngT e 
value  of  T for  optimum  target-separation  performance  cannot  be  found 
without  the  complete  knowledge  of  the  undesired-target  waveform,  the 
value  T-Tfj  should  give  reasonable,  If  not  the  optimum  performance. 
Furthermore,  for  target  identification  based  on  more  than  one  value 
of  T,  the  region  of  T values  chosen  should  contain  T^. 

The  analysis  presented  in  this  section  arc  generally  supported 
by  the  Identification  result^  obtained  in  this  study.  Identifica- 
tion results  are  shown  In  Section  F of  this  chapter. 


* It  is  assumed  that  •<.*  j . Note  that  the  design  of  the  low-pass 
Miter  in  the  preprocessor  of  the  identification  radar  assures 
that  frequencies  higher  than  mg  will  be  highly  attenuated.  Thus, 
"'2  ''“I  d v,1^d  assumption, 
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The  Monti  f i cat  i on  Algorithm 


In  this  study,  a "single-look"  identification  algorithm  was  used. 
A n identification  decision  was  made  for  every  measured  backscattered 
waveform.  The  identification  algorithm  was  one  of  threshold  identi- 
fication based  on  the  values  of  p(T)  evaluated  at  different  values  of 
T,  and  the  average  value  of  p(T)  denoted  by  <p(Tq)>.  Thus,  the  decision 
algorithm  was 


"Desired  Target"  when  , 


n(Toi ) _ PjHi  for  i=1  and  2 

j n(Toi) 

Ufljl  ■ ^ ro)  • ‘ -~-r—  1 nTHA 


and  I 


"Undesired  Target"  when 


p(To1)  < nTHi  for  1=1  or  2 
’THA 


i ' ' oi ' 1 n- 

j ox  *;p(T0)>  < 


or  I 


(5?.) 


pfH-j  ’s  the  identification  threshold  for  p(T)  at  the  T value  of  T0^. 
T0j’s  are  the  values  of  T chosen  for  optimum  identification  perform- 
ance. p.„.  is  the  identification  threshold  for  the  average  of  p(T) 
taken  overall  the  values  of  Tpi.  Thus,  for  a single  measured  wave- 
form, we  evaluate  its  correlation  function  p(T)  at  all  chosen  values 
of  T0i  and  compare  the  values  p ( T0 -j ) and  <o(T0l>  to  pjm  and  pjha* 
respectively  for  an  identification  decision.  Figure  20  illustrates 
the  decision-making  process.  In  Figure  2C,  the  value  of  the  p(T) 
curve  represented  by  "dots"  at  Tg^ , as  well  as  <p(T0)>  are 

greater  than  their  corresponding  thresholds  hence,  the  targets  it  re- 
presents is  by  decision  a "desired"  target.  The  p(T)  curve  represented 
by  "boxes"  has  at  least  one  of  the  values  n(T0i)  or  •■p(Tq)  ■ below  its 
corresponding  threshold,  thus,  the  target  it  represents  is  by  decision 
an  "undesired"  target. 

T0-j  1 s arc  the  value  of  T chosen  for  optimum  identification  per- 
formance. The  values  of  T0i  are  closely  related  to  the  value  of  the 
resonances  of  the  desired  target  and  Equation  (51).  Details  concern- 
ing tne  values  of  T0j's  will  be  further  discussed  in  the  next  section. 
In  Equation  (52),  quotation  marks  are  used  to  signify  a decision,  i.e., 
"Desired  target"  is  by  decision  a desired  target  while  in  reality  it 
could  be  an  nudes  1 red  target  (in  case  of  false  alarm). 

The  identification  algorithm  given  in  Equation  (52)  reduces  to 
its  simplest  fowti  when  only  one  value  of  TQ.  is  considered.  In  this 
case,  tiie  identification  algorithm  is 


‘■Ifcku.hgSgifil 
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Desired  target."  when 

' 1 TH 

Undes i red  target"  when 

,.<v 

■ (,TH 

(53) 


In  the  fallowing  section,  the  predictor-correlator  identifi- 
cation method  is  applied  to  the  processed  Waveforms.  Identification 
performance  will  be  presented  and  discussed. 


F . Per f ortnance  j)_f  the  Predictcr- 

Correl ator  Identifier 

In  order  to  demonstrate  the  effectiveness  of  the  predictor- 
correlator,  a set  of  identification  performance  based  on  a single 
value  of  |.(T)  were  obtained  by  applying  the  predictor-correlator 
alone  (i.e.,  without  the  filter  and  the  detector)  to  the  differ- 
ence waveforms . Identification  results  based  on  multiple  values  of 
,,(T)  and  with  the  filter  and  detector  operations  insetted  will  be 
presented  late  in  this  section.  Identification  performance  was 
characterized  by  a ,.(T)  curve.  Typical  ,.(T)  curves  are  shown  in 
Figure  21.  In  this  case,  the  identifier  was  set  to  identify  the 
mine- like  target  in  a wet  ground  condition. 

From  Figure  21,  we  find  that  typical  r(T)  curves  have  the 
following  features: 

1.  The  value  of  ( (T)  is  approximately  1 for  small  values 
of  T.  Tiiir,  is  true  for  p(T)  of  both  the  desired  and 
undesired  targets.  This  region  of  T has  no  identifi- 
cation capability. 

2.  The  region  of  T.  where  "optimal"  identification  perform- 
ance occurs  is  in  the  immediate  neighborhood  of  Tjj 
where  tw  is  given  by  liquation  (51).  In  this  case  for 
identification  of  the  mine-like  target  in  a wet  ground 
condition.  The  resonances  of  the  mine-like  target  as 
given  in  Table  Ic  was  used.  Tiius,  f^  = 420  MHz,  corre- 
sponding to  the  rounded  T^  value  of  6 Tp  (or  1.2  ns). 

For  the  , • ( T ) curves  given  in  Figure  21,  the  region  of 
T0  is  from  6 Tp  to  11  Tp.  This  result  supported  the 
observations  made  in  Section  0 of  this  chapter.  The 
fact,  that  the  region  of  effective  values  of  T contained 
values  larger  than  T^  indicated  that  the  dominant  polc(s) 
in  the  copper  sheet  waveform  had  imaginary  |\.rts  lower 
than  420  MHz. 

3.  Tiie  value  at  e(T)  fluctuates  when  T is  much  larger  than 

Tm  (T  ■ 2Tm).  This  renion  of  T has  no  identification  capa- 
bility. lies  ides  the  reason  given  in  Section  D of  this 
chapter,  this  can  also  be  attributed  to  the  fallowing. 
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Mrst,  aliasing  of  the  high  frequency  zeros  in  the  iden- 
tifier makes  these  zeros  ineffective.  Second,  the  size 
of  the  error  interval  of  the  error  > decreases  as  T in- 
creases, for  the  number  of  instantaneous  error  samples 
used  in  the  calculation  of  i decreases.  Third,  because 
of  the  large  value  of  T,  only  the  tail  end  of  the  mea- 
sured waveform  is  used  in  the  calculation  of  t.  This 
region  of  the  measured  waveform  has  low  signal  level. 

The  above  features  of  the  e(T)  curves  were  generally  observed 
in  almost  all  cases  considered.  The  significance  of  these  features 
is  that  optimal  i denti f i cat i on  performance  occurs  in  a confined  range 
of  T values  in  the  immediate  neighborhood  of  T^.  Hence  one  knows  a 
priori  the  value  or  range  of  values  of  T on  which  to  base  the  imple- 
mentation of  the  radar  system.  These  features  are  again  demonstrated 
in  the  ,, ( T ) curves  given  in  Figures  22  and  23  for  the  identification 
of  the  mine-like  target  and  the  brass  cylinder  in  dry  ground,  respec- 
tively. 

In  this  study,  the  criterion  for  optimum  identification  perform- 
ance was  to  require  perfect  identification  of  desired-target  returns, 
and  thus  optimal  performance  was  the  one  with  the  lowest  false  alarm 
probability  obtained  by  varying  T.  Single-look  identification  an 
false  alarm  probabilities  are  estimated  in  a frequency  of  occurrence 

sense.  Thus, 


1 = NT 


FA  ' N' 


where 


Pj  Probability  of  identification 
Pc a = Probability  of  false  alarm 


- Total  number  of  desired-target  waveforms  within  the 
identification  range. 

N-,  • Numner  of  desi red-target  signals  which  are  correctly 

identified  within  the  identification  range. 


Total  number  of  other-target  waveforms. 
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Figure  23.  Typical  t.(T)  curve,  for  the  identification 
of  the  brass  cylinder  in  dry  ground. 
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N/j  ■'  Number  of  other-target  waveforms  which  are  mistaken 

to  be  desired-target  waveforms. 

Note  that  N3,  N/j  basically  determine  the  accuracy  or  "confidence  inter- 
val"^] of  the  above  estimates  of  the  identification  statistics. 

To  see  the  degree  of  separation  between  the  mine-like  target  and 
other  subsurface  targets,  consider  the  distribution  of  correlation  co- 
efficients shown  in  Figure  24.  In  this  case,  the  system  is  set  to 
identify  the  mine-like  target  in  wet  ground.  From  Figure  24,  we  see 
that,  for  Pj  = 100%,  and  T0  = 9Td,  the  minimum  separation  margin  be- 
tween the  mine- like  target  and  other  targets  Is  0.93  In  the  value  of 
i>(T0).  In  this  case,  for  an  identification  performance  estimate  of 
Pi  = 1 00% ,* Pfa”=  OX ‘one  can  set  the  identification  threshold  ptH  at  a 
value  ranging  from  -0.26  to  U.l(>7.  The  maximum  value  of  pj^  is  0.167. 
The  correlation  coefficients  for  the  mine-like  target  shown  in  Figure 
24  include  only  cases  measured  within  the  identification  range  RiD. 

In  this  study,  identification  range  is  the  radar  range  measured  from 
the  center  of  the  target.  Correlation  coefficients  for  the  other  tar- 
gets include  all  cases  measured. 

A summary  of  optimal  single-look  identification  performance  is 
given  in  Table  4 (for  details,  see  Appendix  E).  Estimates  of  Pj 
100%,  Ppfl  u 0%  identification  performance  were  obtained  for  the 
identification  range  of  30  cm  measured  from  the  center  of  the  mine- 
like target  and  the  brass  cylinder. 

This  set  of  identification  performance  basically  established 
the  fact  that  the  predictor-correlator  identification  method  worked 
in  its  simplest  form  and  with  a minimum  amount  of  preprocessing. 
However,  in  a practical  system,  one  would  prefer  to  have  the  filter 
and  detector  inserted  for  better  performance  in  more  severe  clutter/ 
noise  environment  provided  that  the  amount  of  time  taken  to  perform 
these  operations  are  not  too  largo  to  make  them  non-cost-effective. 
Furthermore,  tne  multiple-threshold  algorithm  would  also  provide 
better  performance. 

In  target  identification  with  detection  and  multiple-threshold 
operations,  the  ranges  of  the  detection  parameters,  i.e.,  waveform 
energy,  peak  timing  and  peak  magnitude  must  be  chosen.  Furthermore, 
the  values  of  T« -j  as  well  as  the  identification  thresholds  pjhI  anc* 
rjDA  must  be  determined.  In  thws  study,  «bhe>»Q« n—aroeters  were  .ei/fcari,, 
mentally  determined  via  the  following  procedure: 

1.  The  ranges  of  the  detection  parameters  were  chosen  to 
be  the  ranges  of  these  parameters  of  the  waveforms 
measured  from  the  desired  target  within  the  identifi- 
cation range. 

2.  Choose  the  range  of  T0j  values  based  on  the  conditions 
given  in  Section  D of  this  chapter. 


Figure  24.  ’‘Distribution  of  correlation  coefficients  for  the  identification  of  the 
; mine-like  target  in  wet  ground.  Rrr,=30  cm. 
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TABLE  4 

SUMMARY  OF  SINGLE-LOOK  IDENTIFICATION  PERFORMANCE 
OF  THE  SHORT-CABLE  SYSTEM  . 

P j = 1 00%  For  A1 1 Cases 


NUMBER  OF 
WAVEFORMS 

DESIRED  TARGET 

§ 

M 

a h- 

zd  a 
O 

Gd  O 
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To 

PTH 

RID 

(cm) 
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1/1  cc 
lli  «c 
Q t- 

a 

1 LO 

oe 

KG 

LU  O 
Q CC 
X < 
=3  h- 

MINE-LIKE 

TARGET 

WET 

7TB-9TB 

.183 

30 

0% 

12 

12 

MINE-LIKE 

TARGET 

DRY 

7Tb-9Tb 

.483 

15 

0 % 

5 

3 

MINE-LIKE 

TARGET 

ICY 

7TB-9Tb 

.423 

30 

0% 

9 

9 

BRASS  CYLINDER 

ICY 

5TB-7TB 

.675 

30 

0% 

9 

9 

♦PERFORMANCE  FOR  Rj d=30  cm  WAS  NOT  OBTAINED  IN  THIS  CASE. 
R|D  “ IDENTIFICATION  RANGE,  MEASURED  FROM  CENTER  OF  TARGET 


3.  Choose  the  identification  thresholds  to  be  the  minimum 
values  of  r(T0-j)  and  • t> (T0 ) ■ evaluated  from  the  desired- 
target  waveforms  measured  within  the  identification  range. 

An  example  is  given  below  to  clarify  the  above  procedure.  From 
the  waveforms  measured  from  a mine-like  target  in  a wet  ground  con- 
dition, we  evaluated  their  detection  parameters  (as  tabulated  in 
Table  5)  and  determined  the  minimum  and  maximum  of  these  parameters. 
The  ranges  of  these  parameters  were  used  as  the  desired  ranges  of  the 
detection  parameters  in  the  detector.  With  the  extracted  resonances 
from  the  measured  waveforms  (as  given  in  Table  6),  and  the  conditions 
given  in  Section  D of  this  chapter,  we  chose  the  range  of  T0^  values 
folio  from  (TN'-1p)  to*  2T^*.'*4  Tn  this’  cas£  ^=428  ^Hz,  thus 


L2x42Hx10<\ 


* [ J denotes  greatest  integer. 


TABLE  6 

AVERAGE  EXTRACTED  RESONANCES  OF  THE  MINE-LIKE 
TARGET  IN  A WET  GROUND  CONDITION 


73 


POLE 

(REAL)* 

POLE 

(IMAG) 

- 1 . 71 536397E8 

*6.1601933  E7 

-5.73125538E7 

• 1 . 32122486E8 

-2.87494683E8 

*2.270301 33E8 

-1. '82908509  E8 

* 3.0741 1043E8 

-9. 76454733E7 

*4. 28447900E8 

*Rea1  and  imaginary  parts  of  the 
extracted  poles  are  In  nepers/s 
and  Hz , respectively. 


With  the  values  of  Toi , we  evaluate  the  values  of  p(T(?i)  end  <pU'0)> 
for  the  all  measured  desi red- target  waveforms.  The  minimum  of  these 
values  are  chosen  to  be  the  identification  thresholds.  This  is  illu- 
strated in  Table  7. 

Table  8 summarizes  a set  of  single-look  identification  perform- 
ance based  on  the  additional  preprocessing  and  the  mul ti pie-threshold 
algorithm.  Note  that  the  size  of  the  ensembles  of  waveforms  consider- 
ed indicates  a better  estimate  of  the  identification  statistics.  In 
obtaining  this  set  of  identification  performance,  the  region  of  Tp 
values  was  chosen  to  be  from  (TN-T3)  to  2T^.  Details  concenring  the 
identification  results  are  given  in  Appendix  E. 

The  identification  statistics  discussed  so  far  were  obtained 
with  the  identifier  "tuned"  to  the  right  ground  condition,  e.g.,  the 
wet-ground  resonances  of  the  mine-like  target  were  used  to  Identify 
the  mine-like  target  in  that  wet  ground  condition.  If  the  Identifier 
is  not  tuned,  identification  performance  will  degrade,  e.g.,  for  the 
waveforms  ggn^ered  ip  .TaJjlp.  4v  w,lRn.  tt>«  d]\v-j0>;Qnnd.  resonances. .wpre  . . 
used  to  identify  the  mine-like  target  In  wet  ground,  performance  de- 
graded to  Pj=90*,  PpA“ 1 6 . 6% . Methods  for  on-location  calibration  of 
the  ground  condition  and  automatic  tuning  of  the  Identification  radar 
to  ground  condition  is  currently  being  pursued.  One  method  is  based  on 
the  backscattered  waveforms  from  thin  wires  measured  using  the  same 
pulse  radar  system.  On-location  calibration  of  the  ground  condition 
and  automatic  tuning  of  the  identifier  to  the  ground  condition  are 
discussed  in  Chapter  IX. 


Identification  is  degraded  when  the  radar  frequencies  do  not 
properly  span  the  target  resonances.  This  is  demonstrated  in  Chap 
ter  V. 


TABLE  B 

SUMMARY  OF  SINGLE-LOOK  IDENTIFICATION  PERFORMANCE  OF  THE 
SHORT-CABLE  SYSTEM  BASED  ON  ADDITIONAL  PREPROCESSING 
AND  MULTIPLE -THRESHOLD  ALGORITHM  . 

P]n100%  For  All  Cases 


DESIRED  TARGET 


NUMBER  OF 
WAVEFORMS 


a 

uy 

o CXL 
uj  H ►“ 

C£  UJ  CO  UJ 
m CD  UJ  CD 

£3  S3 

Qh  =>  h- 


MINE-LIKE  TARGET  WET.l 


4Tb  . 948 
5TB  .667 
6Tb  -.209 
7T8  .243 
8Tb  -.0615 
9Tb  -.118 
10Tr  .224 


.25937  30  0 % 


9 70 


.30050  30  5.75 
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CHAPTER  V 

EFFECTS  OF  RADAR  BANDWIDTH  ON  THE  CHARACTERIZATION 
AND  IDENTIFICATION  OF  SUBSURFACE  TARGETS 


A.  Introduction 

Th$  .target  resonances  present  in  the  backscattered  waveforms 
depend  on  the  bandwidth  of  the  radar  system.  Target  resonances 
residing  near  the  band  edge(s)  or  outside  the  bandwidth  of  the  radar 
system  are  either  weakly  excited  or  not  excited  at  all.  Thus,  when 
the  radat  bandwidth  is  reduced,  the  number  of  target  resonances  present 
in  the  backscattered  waveforms  will  decrease  accordingly.  In  this 
study,  a set  of  data  was  obtained  with  an  additional  ZOO  m of  con- 
necting cables(RG-8)  inserted  in  the  system  to  keep  the  equipment  in- 
doors during  inclement  weather.  These  cables  acted  as  a low-pass 
filter  with  attenuation  of  30  dB  at  400  MHz  and  thus  reduced  the  radar 
bandwidth.  Such  bandwidth  reduction  was  found  to  degrade  identifica- 
tion performance  for  certain  subsurface  targets  under  consideration. 


B .  Processed  W a veforms  Obtained  With 
the~~Long~-CaFTe~  System 

Typical  processed  waveforms  obtained  with  the  long-r.able  system 
are  shown  in  Figure  2b.  A quick  comparison  between  the  waveforms  In 
Figures  7,  9 and  25  indicates  the  severe  loss  of  high  frequency  in 
the  waveforms  of  the  mine-like  target  and  wood  board.  Furthermore, 
there  was  a significant  loss  of  signal  amplitude  for  all  waveforms. 


C.  Extracted  Resonances 

The  average  extracted  resonances  of  the  long-cable  waveforms 
from  the  subsurface  targets  are  given  in  Tables  9 and  10.  A compari- 
son between  the  extracted  resonances  in  Tables  3 and  9,  10  indicates 
a general  VoT.5^ff  high-frequency  resonances  iii"the  long-cable- wave-  • * « •, 

forms.  The  loss  of  high-frequency  content  causes  the  high  frequency 
resonances  to  be  either  weakly  excited  or  not  excited  at  all.  A 
comparison  between  the  resonances  extracted  from  the  two  sets  of 
mine-like  target  waveformr  (as  plotted  in  Figure  26)  shows  that,  in 
the  long-cable  system,  the  resonance  of  the  mine-like  target  at  300 
MHz  was  not  excited,  while  the  resonance  at  400  MHz  was  only  weakly 
oxc i ted . 
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Figure  26.  Average  extracted  resonanc 
short-cable  mine-like  targ 
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D.  Target  Identification  Performance 

Single-look  identification  performance  is  summarized  in  Table  11.* 
As  expected,  identification  performance  degraded.  The  drastic  degra- 
dation in  the  performance  for  the  identification  of  the  mine-like 
target  and  the  wood  board  clearly  demonstrated  the  importance  of  high 
frequencies  for  the  identification  of  these  two  targets. 


TABLE  11 

SUMMARY  OF  SINGLE-LOOK  IDENTIFICATION  PERFORMANCE  OF  THE 
LONG-CABLE  SYSTEM.  ADDITIONAL  CABLE  LENGTH*200  m . 
Pj  = 1 00%  For  All  Cases 


NUMBER  OF 
WAVEFORMS 

o 

O 

UJ 

DESIRED  TARGET 

GROUND 

CONDIT] 

To 

PTH 

rid 

(cm) 

PFA 

DESIRE! 

TARGET 

UNDESIF 

TARGET 

MINE-LIKE  TARGET 

DRY 

?tb 

.675 

30 

33.55% 

9 

152 

BRASS  CYLINDER 

DRY 

iotb 

.600 

30 

4.49% 

9 

156 

ALUMINUM  SPHERE 

DRY 

IOTq 

.867 

52 

9.29% 

12 

140 

COPPER  SHEET 

DRY 

9Td 

.964 

40 

4.35% 

11 

138 

WOOD  BOARD 

WET 

?tB 

.823 

45 

23.98% 

n 

138 

The  importance  of  the  high-frequency  resonances  of  the  mine-like 
taryet  stems  from  the  fact  that  its  high-frequency  resonances  are 
dominant  (see  residues  of  target  resonances  in  Appdenix  C).  An  attempt 
to  identify  this  target  without  strongly  exciting  its  dominant  reson- 
ances results  in  the  poorer  identification  performance.  To  improve 
target  identification  performance,  the  radar  must  be  able  to  transmit 
and  receive  more  high-frequency  energy  in  the  vicinity  of  these  high- 
order  resonances,  the  short-cable  system  discussed  previously  repre- 
sented an  improvement  over  the  long-cable  system.  However,  much  more 


*This  sot  of  identification  results  was  obtained  based  on  threshold 
identification  with  one  value  of  T0  and  without,  detection  and 
f i ltering. 


improvement  can  bo  nuulo  in  the  various  components  of  the  pulse  radar 
system  to  provide  a better  "match"  between  the  system  bandwidth  and 
the  bandwidth  in  which  the  dominant  target  resonances  reside.  One 
of  these  components  is  the  antenna  system. 

For  the  identification  of  the  plastic  mine-like  target,  the 
crossed-dipole  antenna  with  arm  length  of  0.6  m (2  feet)  did  not  pro- 
vide a "good"  match  between  the  radar  system  bandwidth  and  the  band- 
width of  the  target  resonances.  The  extracted  resonances  shown  in 
Figure  14  indicated  that  the  antenna  resonance  was  not  even  near'  Ute 
high-order  resonances  of  the  mine-like  target.  Such  mismatch  of 
bandwidths  results  in  the  low  level  of  correlation  coefficients  as 
tabulated  in  Table  8.  For  better  characterization  and  identification 
of  the  mine-like  target,  the  frequency  of  the  antenna  resonance  needs 
to  be  increased.  One  easy  way  to  accomplish  such  increase  is  to 
reduce  the  size  of  the  antenna.  In  this  study,  a smaller  crossed- 
dipole  antenna  with  arm  length  of  0.15  m (0.5  feet)  was  built  to  pro- 
vide a better  match  of  bandwidLhs  for  the  identification  of  the  mine- 
like  target.  Mine  identification  is  the  subject  of  Chapter  VII. 

In  the  next  chapter,  we  study  the  effects  of  target  depth  and 
size  on  the  characterization  and  identification  of  subsurface  targets. 
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CHAPTER  VI 

EFFECTS  OF  TARGET  DEPTH  AND  SIZE  ON  THE  CHARACTERIZATION 
AND  IDENTIFICATION  OF  SUBSURFACE  TARGETS 


A . Introduction 

To  study  the  effects  of  target  size  and  depth  on  the  character- 
ization and  identification  technique,  two  sets  of  targets  were 
buried  (see  Figures  27  and  28).  The  first  set  consisted  of  a series 
of  different-size  brass  cylinders  buried  at  different  depths.  The 
second  set  consisted  of  a series  of  different-length  0.3125  cm  (1/8- 
inch)  diameter  thin  brass  wires  buried  at  the  depth  of  5 cm  (2 
inches).  Backscattered  waveforms  were  obtained  using  the  subsur- 
face pulse  radar  at  various  antenna  locations.  Locations  of  the 
antenna  center  are  shown  as  dots  in  Figures  27  and  28. 
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13.  Processed  Waveforms 

Processed  waveforms  from  the  cylinder  and  the  wire  targets  are 
shown  in  Figures  29-31.  From  these  waveforms,  the  following  observa- 
tions are  made’. 

1.  All  these  waveforms  exhibited  some  transient  behavior, 
signifying  the  possible  existence  of  one  or  more  natural 
resonances. 

2.  A comparison  between  the  brass  cylinder  waveforms  shown 
in  Figures  9 and  29  indicates  that  the  signal  level  of 
the  waveform  from  the  brass  cylinder  at  5 cm  depth  was 
appraxinetely  12  dB  higher  than  that  of  the  waveform 
from  the  brass  cylinder  at  30  cm  depth.  This  signified 
a propagation  loss  of  over  10  dB  ner  30  cm  (1  foot). 

3.  A burst  of  signal  energy  appeared  in  the  early-time 
portion  of  the  waveforms  from  the  brass  cylinders  at 
90  and  150  cm  depth.  This  was  caused  by  a phenomenon 
known  as  the  "Trench  effcct"[22]. 

Rain,  snow  and  evaporation  changed  the  moisture  content 
of  the  ground  ar<d  the  distribution  of  moisture  was  per- 
turbed by  the  Trench  walls.  At  times  during  the  course 
of  a year  the  ground  in  the  trench  was  found  by  direct 
measurement  to  be  drier  and  at  other  times  wetter 


Figure  27.  Layout  of  i'..~  fi'-ied  cylinders 
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30  cm  LONG 
THIN  WIRE 


X ANTENNA 
ORIENTATION 


45  cm  LONG 
THIN  WIRE 


DEPTH  OF  ALL  THIN  WIRES  5 cm  ( 2 INCHES  ) 
DIAMETER  OF  ALL  THIN  WIP  .S  ■ 0.31  cm  ( 1/8  INCH  ) 


Figure  28.  Layout  of  the  buried  thin  wires. 
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than  the  ground  outside  the  trench.  Thus,  the  trench  be- 
came a scatterer.  The  trench  signal  obscured  the  cylinder 
signal  and  made  the  task  of  target  identification  difficult. 

4.  Careful  study  of  the  waveforms  of  Figure  31  from  the  thin 
wires  indicated  that  the  time  interval  between  zero  cross- 
ings increased  according  to  the  increase  in  the  wire  size. 
Such  phenomenon  was  not  observed  in  the  cylinder  waveforms. 


C.  The  Extracted  Resonances 

1 . Effects  of  target  depth 

To  see  the  effects  of  depth  on  the  complex  resonances,  the 
average  extracted  resonance  of  the  30cm  long  cylinders  at  different 
depths  are  shown  in  Figure  32,  From  these  resonances,  we  make  the 
following  observations: 

a.  The  high-frequency  content  in  the  backscattered  waveforms 
was  highly  attenuated  as  target  depth  increases.  The 
highest-order  resonance  which  was  present  in  the  5 and 
30cm  deep  cylinder  waveforms  was  absent  from  the  wave- 
forms of  the  deeper  cylinders. 

b.  All  resonances  with  imaginary  parts  smaller  than  400  MHz 
were  present  in  the  cylinder  waveforms  of  all  depths. 

c.  The  lowest-order  resonance  was  the  antenna  resonance 
(0.6m  long  antenna).  This  resonance  was  present  in  almost 
all  the  waveforms  of  all  targets  at  all  depths  collected 
using  this  antenna. 

d.  The  imaginary  part  of  the  target  resonances  of  the  brass 
cylinder  occurred  in  integer  multiples  of  the  imaginary 
part  of  the  lowest-order  target  resonance.  This  seemed 
to  suggest  that  the  extracted  resonance  of  the  brass 
cylinder  were  caused  by  the  multiple  scattering  mechanisms 
along  the  length  of  the  cylinder  only  (i.e.,  the  dipole 
mode) [21],  the  creeping-wave  modes  were  not  present  in  the 
waveforms.  The  absence  of  the  creeping  wave  type  reson- 
ances was  further  emphasized  when  we  compared  the  extracted 
resonances  of  the  same-size  thin  wire  at  the  depth  of  5 cm 
to  the  cylinder  resonances  (see  Figure  32).  The  same  lower- 
order  resonances  were  extracted.  This  proved  the  absence 

of  the  creeping-wave  modes,  for  the  creeping-wave  modes 
of  the  thin  wires  have  resonant  frequencies  which  are  too 
high  to  be  present  inside  the  bandwidth  of  our  radar  sys- 
tem. The  absence  of  the  creeping-wave  modes  is  attributed 
to  the  high  loss  suffered  by  the  creeping  waves  as  it 
traveled  around  the  circumference,  and  the  cross-polari- 
zation effects  of  the  crossed-di pole  antenna  system. 


91 


With  the  absence  of  the  creeping-wave  type  resonances, 
from  the  point  of  view  of  the  target  resonances,  the  thin 
wire  and  the  cylinder  are  indistinguishable.  Separation 
of  these  two  targets  is  possible  only  when  discriminants 
ether  than  the  lower  resonances  are  used.  One  such 
possible  discriminant  is  the  forced  response  of  the  back- 
scattered  waveforms  which  contains  information  about  the 
profile  area  of  the  target[48]. 

e.  The  fact  that  the  extracted  resonances  of  the  30cm  long 
cylinder  are  related  solely  to  the  length  of  the  cylinder 

suggests  the  use  of  these  resonances  to  estimate  the 
parameters  of  the  ground.  If  we  consider  these  resonances 
to  be  an  approximation  to  the  resonances  of  the  same-size 
cylinder  in  a homogeneous  medium  illuminated  by  a plane 
wave,  the  imaginary  parts  of  these  resonances  can  be 
approximated  by  Equati  on  (54)[21  ] 


fn  - n — c_  ; n = 1 ,2,3*  • • (54) 

2L/rr 

where  L is  the  physical  length  of  the  cylinder,  c is  the 
speed  of  light  in  free  space  and  rr  is  the  relative  die- 
lectric constant  of  the  medium.  Using  Equation  (54)  and 
the  extracted  resonances  of  the  30 cm  long  cylinder  as 
given  in  figure  32  the  relative  dielectric  constant  is 
estimated  to  be  approximately  20  at  the  resonant  frequency 
of  approximately  100  MHz.  This  estimate  agrees  closely 
with  the  estimate  from  other  methods  currently  being 
investigated  in  the  ElectroScience  Laboratory£58] . 

The  complex  resonances  of  the  30cm  long  cylirders  and  thin 
wire  relate  very  simply  to  their  physical  size.  This  simple  rela- 
tionship disappears  when  the  size  of  the  target  increases  beyond  a 
certain  threshold.  This  is  demonstrated  in  the  following  discussion. 


2.  Effects  of  target  size 

As  target  size  increases,  we  expect  the  target  resonances  to  be 
lower  in  frequency.  Such  expectation  turns  out  to  be  warranted  when 
we  consider  the  extracted  resonances  of  the  different-size  thin  wires 
at.  the  depth  of  5 cm  as  shown  in  Figure  33.  Again,  the  imaginary 
parts  of  target  resonances  occur  in  simple  multiples  of  the  imagi- 
nary parts  of  the  lowest -order  target  resonance,  and  as  tiie  length  of 
the  wire  increases  up  to  60  cm,  the  frequency  decrease  according  to 
its  size.  This  expectation  turns  out  to  be  unwarranted  when  the 
length  of  the  cylinder  increases  to  over  90  cm  at  the  depth  of  30  cm. 
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Extracted  resonances  of  the  different- size  cylinders  at  the  depth 
of  30  cm  are  shown  in  Figure  34. 

Although  resonances  of  the  large-size  targets  are  no  longer 
related  to  their  sizes  in  simple  form,  however,  there  are  still 
resonance:  pres°nt  in  the  backscattered  waveforms.  Hence  identi- 
fication of  these  targets  is  still  possible.  Some  identification 
results  are  presented  in  the  next  section. 


O 30cm  LONG 
□ 90  cm  LONG 
A 150cm  LONG 
V 300cm  LONG 


A y 


CYLINDER  DEPTH- 30cm 
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Figure  34.  Average  extracted  resonances  of  the  different 
size  cylinders  at  30  cm  depth. 
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D.  Target  Identification 

Some  single-look  identification  performance  results  are  given 
in  Table  12.  Based  on  these  results,  the  following  observations 
are  made: 

1.  All  false  alarm  probabilities  were  less  than  Q% . In  some 
cases,  estimates  of  Pi=100%,  Pp^=0%  were  obtained. 

2.  Identification  range  decreased  as  target  depth  increased. 

In  the  next  chapter  we  focus  our  attention  on  the  detection 
and  identification  of  mine-like  targets  in  more  extensive  and  practi- 
cal situations.  Improvements  over  the  Terrascan-1 i ke  pulse  radar  are 
made  for  the  implementation  of  a portable  on-location  real-time  sub- 
surface target  identification  radar. 


TABLE  12 

SUMMARY  OF  SINGLE-LOOK  IDENTIFICATION  PERFORMANCE 
FOR  THE  BRASS  CYLINDERS  OF  DIFFERENT  SIZES 
P i =1  OUVo  For  All  Cases 


BRASS  CYLINDER 


1 LNGTH 
(cm) 

DEPTH 

(cm) 

30 

JO 

90 

30 

150 

30 

300 

30 

300 

90 

300 

150 

rid 

(cm) 

South-  East-* 
North  West 

PFA 

NUMBER  OF 
WAVEFORMS 

DESIRED 

TARGET 

UNDESIRED 

TARGET 

30 

30 

3 , orvr- 

11 

65 

60 

30 

7.25m 

9 

69 

105 

30 

Om 

17 

cn 

150 

30 

0'.', 

9 

65  1 

150 

15 

7 . 9* 

7 

38 

150 

15 

OX 

9 

*See  Figure  27. 


CHAPTER  VII 

ELECTROMAGNETIC  MINE  DETECTION  AND  IDENTIFICATION 


A.  Objectives 

This  chapter  focuses  on  the  detection  and  Identification  of  mine- 
like targets.  First  a second  model  of  the  mine-like  target  was  buried 
at  the  same  depth  of  5 cm.  Measurements  were  obtained  over  the  two 
isolated  mine-like  targets  for  comparison  and  identification.  Second, 
to  simulate  a realistic  mine  identification  situation  a set  of  wave- 
forms was  obtained  over  a section  of  a rough  road  in  which  numerous 
false  targets  existed.  These  false  targets  were  caused  by  the  debris 
existing  In  the  old  country-style  road,  and  thus  represented  a set  of 
realistic  false  targets  in  a mine  field.  The  set  of  rough  road  measure 
men t was  used  for  evaluation  of  "typical"  false-target  discrimination. 
Third,  a small  antenna  of  arm  length  0,15  in  (0.5  foot)  was  built  and 
used  to  provide  better  characterization  and  identification  of  the  mine- 
like target.  Performance  of  the  smal 1 -antenna  system  is  presented  and 
discussed. 


B . A Second  Model  of  the 

Mine-like  Target 

Typical  processed  waveforms  from  the  two  mine-like  targets  at 
the  same  relative  antenna-target  geometry  are  shown  in  Figure  35, 

There  are  minor  differences  in  these  two  waveforms.  This  can  be 
attributed  to  the  minor  difference  in  the  structure  of  these  two 
targets,  clutter  and  the  difference*  in  the  ground  conditions  at 
the  two  target  locations.  Complex  natural  resonances  of  the  two 
models  of  the  mine-like  target  were  extracted  from  their  backscattered 
waveforms  using  Prony's  method.  The  average  extracted  resonances  of 
the  two  models  are  shown  In  Figure  3b.  Again  the  slight  difference 
in  the  locations  of  these  two  sets  of  resonance  is  attributed  to  the 
possible  slight  difference  in  the  structure  of  the  two  models,  clutter 
and  the  difference  in  the  ground  conditions  at  the  two  target  sites. 
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Figure  36.  Average  extracted  resonances  of  the  two 
mine-like  targets. 


f ( 10°  Hr  ) 


C . False  Targcl  Measurements 

i n the  Rough  Road 

A sot  of  53  (difference;  waveforms  was  collected  aver  53  loca- 
tions in  a 40m  section  of  a rough  road  adjacent  to  the  Electro- 
Science  Laboratory.  These  echoes  were  probably  a result  of  debris 
existing  in  this  old  country-style  road.  There  were  numerous  such 
false  targets  observed.  Pictures  of  the  rough  road  is  shown  in 
Figure  38- 

The  measurement  locations  were  chosen  to  ensure  significant 
false-target  signal  levels  for  a more  realistic  and  difficult  test 
of  the  identification  method.  Figure  38  shows  a typical  false-target 
waveform. 


D.  Mine  Identification 


The  predictor-correlator  identification  method  was  applied  to 
the  waveforms  from  the  two  mine-like  targets  and  the  rough  road-bed 
For  identification  of  the  two  mine-like  targets  and  false-target  dis- 
crimination. Single-look  identification  statistics  are  given  in 
Table  13. 

This  set  of  identification  statistics  basically  established  the 
fact  that  the  predictor-correlator  identification  method  did  Indeed 
successfully  Identify  different  models  of  the  mine-like  target  at 
different  locations.  Furthermore  the  method  successfully  separated 
the  mine-like  target  from  the  set  of  false  targets  which  were  typi- 
cal of  the  false  targets  found  in  a realistic  mine  field. 

The  pulse  radar  system  discussed  so  far  represented  a workable 
(and  successful)  mine  identification  system,  however  two  obvious 
improvements  could  be  made.  First,  there  was  the  bandwidth  mismatch 
problem  (as  discussed  in  Chapter  V)  that  resulted  in  low-level  values 
nf  the  correlation  coefficient.  The  bandwidth  mismatch  problem  could 
be  corrected  by  Introducing  an  antenna  system  with  a higher-frequency 
antenna  resonance.  Second,  the  size  of  the  crossed-dipole  antenna 
(1.2  m (4  feet)  maximum  dimension)  was  a bit  too  large  for  a practical 
mine  identification  system.  There  existed  a single  solution  to  the 
two  problems  mentioned  above.  This  solution  was  the  reduction  in  the 
antenna  size.  Since  the  0.6m  long  (arm  length)  antenna  had  a reson 
ance  at  65  MHz,  a 0.16m  (0.5  foot)  long  antenna  would  have  a resonance 
at  260  MHz  which  was  close  to  the  mid -frequency  of  the  resonance  band- 
width. A crossed-dipole  of  0.15  in  arm  length  was  therefore  built  for 
improved  mine  identification  performance[59] . Performance  of  this 
antenna  on  the  identification  of  the  mine-like  target  is  discussed  in 
the  next  section. 
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F.  A Small  An  ton  n<i  for  Improved  Performance 

1 1 Mm  I denti  fi cation 

A picture  of  the  small  crossed-dipole  antenna  is  shown  in  Figure 
39.  This  antenna  was  used  to  obtain  a set  of  measurements  over  the 
5cm  deep  targets.  Typical  processed  waveforms  are  shown  in  Figure  40. 
For  comparison  purposes,  Figure  41  shows  the  processed  waveforms  ob- 
tained using  the  0.6m  long  antenna  over  the  same  target  locations 
From  Figures  40  and  41,  we  note  the  following  improvements  obtained 
using  the  0.15m  long  antenna.  First,  the  0.15m  long  antenna  offers 
a gain  of  at  least  6 dB  in  target  signal  level.  Second,  the  level  of 
the  no-target  signal  obtained  using  the  0.15m  long  antenna  is  gener- 
ally lower.  Third,  the  time  interval  of  the  target  signals  obtained 
using  the  0.15m  long  antenna  is  only  half  as  long.  This  reduction 
in  interval  size  would  effectively  reduce  the  computation  time  for 
evaluating  the  correlation  coefficient  by  a factor  of  two. 

Due  to  the  minimization  of  the  length  of  the  inter-connecting 
cables  at  the  balun-antenna  connection  of  the  small  antenna,  the  balun 
reflection  exists  in  the  early-time  portion  of  the  waveform.  In  this 
region,  the  primary  target  signal  level  is  much  higher  than  the  level 
of  the  balun  reflection.  This  severely  suppresses  the  effects  of  the 
balun  reflection.  The  decrease  in  cable  length  is  also  a contributing 
factor  to  the  6 dB  gain  in  signal  level.  With  the  balun  reflection 
moved  in,  the  stop-time  te  of  the  error-calculating  interval  in  the 
Prony's  and  the  predictor-correlator  identification  processes  had  to 
be  changed.  In  this  study,  the  choice  of  to  for  the  small-antenna 
system  was  based  on  the  resonances  of  the  mine-like  target.  The  value 
of  te  was  chosen  to  yield  an  error  Interval  of  size  equal  to  five 
times  the  period  of  the  lowest  imaginary  part  of  the  target  resonances 
(125  MHz).  This  resulted  in  the  te  value  of  ts+100Tg. 

Figures  42-44  shows  more  waveforms  and  their  FFT's  obtained 
using  the  0.15m  long  antenna.  From  the  FFT's  of  the  waveforms,  it 
is  found  that  most  energy  resided  in  the  frequency  region  of  100  M 
to  500  MHz  with  the  most  dominant  frequency  in  the  region  of  about 
300  MHz  (as  compared  to  70  MHz  for  the  0.6m  long  antenna). 

The  set  uf  waveforms  obtained  with  the  small  antenna  over  the 
5cm  deep  targets  was  processed  for  pole  extraction  and  identifica- 
tion of  the  mine-like  target.  Average  complex  resonances  of  the  mine- 
like target  extracted  from  these  waveforms  are  shown  in  Figure  45. 

For  comparison  purposes,  the  extracted  resonances  from  the  mine-like 
target  waveforms  obtained  with  the  0,6m  long  antenna  in  a similar 
ground  condition  are  also  shown  in  Figure  45.  The  locations  of  the 
target  resonances  from  both  antennas  were  almost  identical.  However, 
the  imaginary  part  oF  the  antenna  resonance  increased  by  a factor  of  4 
when  the  length  of  the  antenna  arm  was  reduced  from  0.6  m to  0.15  ru. 


Figure  39.  The  small  crossed-dipo 
identification  perform 


Processed  waveforms  from  5cm  deep  targets  measured  using 
the  0.6m  long  antenna. 


Figure  42.  Addit 
5 cm  d 


antenna  waveforms  shown  in  Figure  40 
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Figure  45.  Average  extracted  resonances  from  the  mine-Hke 
target  waveforms  taken  using  the  0.15m  and  the 
0.6m  long  antennas  In  similar  ground  conditions. 
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Compl ox  resonances  of  the  various  5cm  deep  targets  are  shown 

in  Fiqure  46.  The  antenna  resonance  was  present  In  almost  all  target 
waveforms  measured.  The  complex  resonances  of  these  targets  were  in 
the  same  general  area  of  the  complex  frequency  plane  as  the  antenna 
resonance.  This  was  expected  to  present  a difficult  test  to  the 
identifier. 

Single-look  identification  statistics  obtained  using  the  small 
antenna  in  identification  of  the  mine-like  target  is  shown  In  Table 
14.  Typical  p(T)  curves  are  shown  in  Figure  47.  From  Table  14  and 
Figure  47,  the  following  observations  are  made: 

I.  Target  identification  performance  estimates  based  on  the 
extracted  resonances  of  the  mine-like  target  was  P]=100%, 
PFA~1 • 72%  within  a 30  cm  radius  of  target  center.  Within 
a 45  cm  radius,  identification  performance  was  Pj=l00%, 
Pfa^^O1^  This  set  of  statistics  was  obtained  with  an 
ensemble  of  13  mine-like  target  and  58  false-target  wave- 
forms. These  performances  represented  an  Improvement  over 
the  performance  obtainable  with  the  0.6m  long  antenna. 

The  0.6m  long  antenna  system  was  unable  to  identify  the 
mine-like  target  In  a distance  of  more  than  30  cm  from 
the  target  center. 

?.  In  addition  to  the  range  improvement,  a quick  comparison 
between  the  identification  statistics  given  in  Tables  8 
and  14  indicates  that  the  level  of  p(T)  values  for  the 
mine-like  target  were  generally  Improved.  It  is  expected 
that  this  will  produce  improved  Identification  in  the  pre- 
sence of  greater  amounts  of  clutter.  It  has  already  been 
observed  that  target  identification  in  the  presence  of 
clutter  has  been  improved. 

3.  The  region  uf  T0  values  for  optimum  identification  per- 
formance was  in  the  neighborhood  of  T^=6T g . However, 
this  region  of  T0  values  was  shifted  toward  small  T0 
values.  This  was  attributed  to  the  shift  to  dominance  of 
to  the  high-frequency  poles  due  to  the  shift  in  the 
antenna  resonance.  Note  that  the  region  of  To  Included 


Because  of  the  change  in  the  bandwidth  of  the  small -antenna 
system  a band-pass  filter  was  used  in  the  preprocessor  of 
the  identifier  to  suppress  out-of-band  clutter  and  noise. 
The  transfer  function  of  the  band-pass  filter  is  shown  in 
Figure  48. 

In  the  next  chapter,  we  discuss  the  implementation  of  the  small 
antenna  mine  identification  system  as  a microcomputer  system.  Target 
identification  performance  will  be  given. 


f ( 10®  Hz) 


TABLE  14 

INGLE-LOOK  IDENTIFICATION  PERFORMANCE  FOR  IDENTIFICATION 
OF  THE  MINE-LIKE  TARGET  WITH  THE  SMALL- ANTENNA  SYSTEM  . 
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Figure  47.  Typical  (.(T)  curves  for  the  Identification  of  the 
mine- like  target  using  the  smal 1 -antenna  system. 


CHAPTER  VIII 

A MICROCOMPUTER  SYSTEM  FOR  REAL-TIME  ON-LOCATION 
SUBSURFACE  TARGET  IDENTIFICATION 


A . Objectives 

This  chapter  discusses  the  Implementation  of  the  pulse  radar 
identification  system  as  a microcomputer  system  for  real-time  on- 
location  Identification  of  subsurface  targets.  The  identification 
radar  Is  Implemented  with  an  SDK-80  m1croconiputer[60]  In  conjunction 
with  the  Torrascan.  The  processing  algorithms  are  stored  In  the 
system  as  microprograms.  Arithmetic  operations  are  performed  by  a 
hardware  arithmetic  processing  unit  and  data  can  be  recorded  by  a 
cassette  recorder.  Thus  the  microcomputer  system  can  serve  either  as 
a real-time  Identification  system  or  as  a data-recordlng  system. 


B . Structure  of  the  Microcomputer 

targe t~'I don't i flcatlon  Sys tern 

The  basic  structure  of  the  microcomputer  target  identification 
system  Is  shown  in  Figure  49.  The  major  components  and  their  corre- 
sponding functions  are  the  following: 

1.  Terra  scan  System:  the  raw  backscattered  waveform  is 

“receTv’ed  ¥v"  the  crossed-dl pole  antenna  In  analog 
form.  This  analog  waveform  is  sampled  by  the 
Torrascan  system  which  essentially  serves  as  a sampl- 
ing oscilloscope  and  provides  various  flexibilities 
such  as  gain,  time  delay,  and  dc  offset  to  the  wave- 
form. The  sampled  raw  waveform  Is  displayed  on  the 
neon  display  and  then  fed  to  the  analog  board  as  a 
discrete-time  signal . 

2.  Analog  Board:  besides  converting  the  sampled  raw  wave- 

form"  digital  form  using  an  8-bit  A/D  converter, 
the  analog  board  provides  a keyboard  for  Inputting 
to  the  microcomputer  system.  Various  operating 
modes  of  the  microcomputer  system  are  Initiated 
through  the  command  codes  entered  from  the  keyboard. 
The  analog  board  also  provides  6 LED  display  lights 
f 0"*  displaying  echoes  of  the  input  command  codes, 
various  system  operating  information,  and  important 
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Figure  49.  U lock  diagram  of  the  microcomputer  system 
for  on-location  subsurface  target 
identification  in  real  time, 


information  concerning  the  received  waveform  (such 
as  peak  timing,  peak  va'ues,  etc.). 

DK-SO : this  is  the  microcomputer  unit.  All  operating 

system  programs,  signal-preprocessing,  detection  and 

identification  and  data-recordlng  operations  are 
stored  as  microprograms  in  the  SDK-80[60].  These 
microprograms  are  executed  when  the  proper  commands 
are  entered  from  the  keyboard.  The  SDK-80  is  an  8- 
bit  machine  running  at  a clock  rate  of  2 MHz. 

4.  Arithmetic  Processing  Unit  (APU):  the  hardware  arithmetic 
processimf  ulilt  is  interfaced  to  the  SDK-80  micro- 
computer for  high-speed  arithmetic  calculations.  The 
APU  provides  various  flexibilities  such  as  floating- 
point or  fixed-point  arithmetic,  single  or  double  pre- 
cision, function  (such  as  trigonometric  functions, 
etc.)  generations,  etc. 
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5.  Memory  Board:  8K  memory  bytes  (8  bits/byte)  for  data 

storage. 

6.  Cassette  Recorder:  used  for  data  recording 

Similar  systems  had  been  built  and  documented  under  different 
studies[67]  at  the  ElectroScience  Laboratory.  Except  for  the  addition 
of  the  hardware  arithmetic  processing  unit,  the  system  developed  and 
implemented  in  this  dissertation  has  basically  the  same  design  and 
hardware.  Thus,  only  the  details  concerning  the  APU  are  given  in 
Appendix  I.  Details  concerning  the  other  modules  of  the  microcomputer 
system  are  well  documented  in  various  reports[67]. 


C.  Implementation  of  the  Microcomputer  System 

for  Subsurface  Target  Identification 

Pictures  of  the  microcomputer  system  and  its  components  are 
shown  in  Figures  50  and  51,  Two  major  operations,  namely  target 
Iflonti Mention  m\  data  recording  wero  Implemented  in  the  microcomputer 
system,  The  target  identification  process  performed  by  the  microcom- 
puter is  outlined  in  Figure  52.  All  processes  except  the  automatic 
tuning  process  have  been  implemented.  The  tuning  process  is  discussed 
in  Chapter  IX.  In  short,  when  the  target  identification  mode  is  being 
executed,  the  Terrascan  samples  the  received  analog  waveform,  displays 
and  feeds  it  to  the  analog  board  for  A/D  conversion.  The  digital 
waveform  is  then  fed  to  the  memory  board  where  the  average  waveform 
is  formed.  The  average  waveform  goes  to  the  preprocessing  unit  for 
reduction  in  clutter/noise.  This  preprocessed  waveform  arrives  at  the 
detector  for  a screening  operation.  The  detector  screens  out  wave- 
forms whose  energy,  peak  timing  and  peak  amplitude  are  out  of  the  de- 
sired ranges.  The  "in-range"  waveform  would  then  be  fed  to  the  pre- 
dictor-correlator for  target  identification.  "Out-of-range"  waveforms 
are  considered  as  "undesired-target"  waveforms.  The  correlation  coef- 
ficients and  the  processed  waveforms  are  stored  in  the  8K  memory  which 
is  transferred  to  the  cassette  recorder  when  filled.  For  target 
identification  based  on  the  predictor-correlator  method,  only  the 

difference  equation  coefficients  corresponding  to  the  desired  target 
need  to  be  stored  in  the  microcomputer  for  the  evaluation  of  p(T),  the 
complex  resonances  are  not  stored.  For  multiple-threshold  identifica- 
tion, a set  of  coefficients  corresponding  to  each  chosen  value  of  T0,- 
needs  to  be  stored.  Note  that  these  coefficients  depend  on  ground 
condition,  thus,  for  on-location  subsurface  target  identification  in 
real  time,  an  automatic  tuning  process  is  necessary.  Automatic  tuning 
cf  the  identification  radar  to  the  ground  condition  is  discussed  in 
Chapter  IX. 

A list  of  executable  commands  and  their  detailed  descriptions 
are  given  in  Table  31  of  Appendix  I. 
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Figure  51.  The  microcomputer  and  Its  components. 
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The  following  notes  with  regard  to  the  Impl ementation  of  the 
microcomputer  identif Icatlon  system  are  important: 

1.  Due  to  hardware  constraints  in  the  Terrascan,  a waveform 
obtained  with  the  microcomputer  system  contained  only  128  data  points. 
The  waveforms  obtained  with  the  digital  computer  controlled  system 
contained  256  data  points.  Thus,  we  were  faced  with  possible  loss  of 
resolution  in  time  (i.e.,  a larger  sampling  period  must  be  used  in 
order  for  128  samp1 es  to  cover  a 50  ns  time  window)  or  frequency  (I.e., 
smaller  time  window  must  be  used).  Careful  study  of  the  srnal 1 -antenna 
waveforms  from  the  mine-like  target  revealed  the  fact  that  the  energy 
of  practically  all  waveforms  resided  in  a time  window  of  25  ns  (see 
Figure  40).  Hence,  by  choosing  a time  window  for  25  ns  and  128  samples 
for  the  mine  identifier,  both  time  and  frequency  resolution  were  main- 
ta  i ried . 


2.  Besides  frequency  and  time  resolution  of  the  waveforms,  the 
other  important  considerations  are  data  precision,  flexibility  and 
speed.  In  this  system,  double-precision  floating-point*  arithmetic 
is  used  to  perform  the  major  calculations  such  as  the  evaluation  of 
the  correlation  coefficient  p(T)  and  the  Fill  filtering  in  the  prepro- 
cessing unit.  The  floating-point  format  offers  large  dynamic  range, 
flexibility  and  ease  of  Implementation.  The  function-generation  capa- 
bility (such  as  trigonometric  functions,  transcendental  functions,  etc., 
arc  available  ftr  floating-point  numbers  only)  of  the  API!  provides  tre- 
mendous flexibilities  In  the  Implementation  of  digital  signal  processing 
techniques.  These,  together  with  the  fact  that  scaling  operation  is  not 
needed  for  floating-point  operations,  makes  t.he  floating-point  format 
an  ideal  choice  for  a "first-generation"  computer  system.  The  tremend- 
ous amount  of  scaling  operations  in  the  fixed-point  calculations  makes 
the  calculations  difficult  to  track  and  furthermore,  the  great  amount 
of  scaling  often  causes  loss  in  data  precision.  The  trade-off  in  us- 
ing the  floating-point  format  is  the  computing  speed.  Floating-point 
operations  usually  take  more  computing  time  than  the  corresponding 
fixed-point  operations.  Furthermore,  the  25-bit  mantissa  offers  less 
"nominal"  data  precision  than  the  32-bit  fixed-point  formated  number. 
With  trie  choice  of  data  format,  we  now  estimate  the  computing  time 
needed  for  an  identification  decision.  The  major  areas  requiring 
significant  amount  of  computations  are  the  following: 


*FToa ting-point  number  in  this  system  is  represented  by  a 7-bit 
exponent  and  a 25-bit  mantissa. 
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a .  Arithmetic  averaging 

Titis  process  forms  the  arithmetic  average  of  ,’Ji  raw  waveforms  from 
the  same  antenna  location  and  orientation  for  reduction  of  noise.  Here, 
the  major  part  of  the  computation  time  goes  to  the  sampling  of  the  N] 
waveforms.  The  sampling  time  At  is  inversely  proportional  to  the  repeti- 
tion frequency  of  the  radar  source,  i.c.,  the  Impulse  generator,  and 


where 

Ni  = number  of  waveforms  required  t.o  form  the  average, 

Nj  » number  of  samples  in  a waveform,  and 

Rp  ° repetition  frequency  of  the  impulse  generator. 


(55) 


In  the  present  system,  Nj-10,  N2»128  and  Rpa256.  Thus  At»5  seconds. 
If  the  repetition  frequency  can  be  increased  to  in  KHz*,  At  will  be 
decreased  to  0.125  second. 


b .  9Qu-difference 


This  process  forms  the  difference  between  two  waveforms  from  two 
antenna  orientations  (one  of  which  is  a 90°  rotation  from  the  other) 
at  the  same  antenna  location.  The  differencing  operation  basically 
increases  the  sampling  time  by  a factor  of  2.  In  the  final  system 
applications,  this  process  Is  probably  not  needed.  At  any  rate,  the 
antenna  rotation  process  will  represent  a major  time  factor. 

c.  Calculation  of  n(T) 

This  process  evaluates  the  correlation  coefficient  p(T)  at  all 
chosen  values  of  Toi's  for  multiple-threshold  identification,  Note 
that  with  the  choice  of  floating-point  format,  the  expression  for  p(T) 
given  in  Equation  (24)  can  be  manipulated  for  minimum  computation  time. 
In  floating-point  calculations,  due  to  the  shifting  of  exponents  re- 
quired for  the  add/subtract  operations,  the  amount  of  time  required  to 
perform  an  add/subtract  operation  can  be  larger  than  that  required  for 
a multiply/divide  operation.  Table  15  lists  the  computation  time  re- 
quired by  the  APU  to  perform  each  of  the  four  basic  floating-point 
arithmetic  operations.  The  maximum  difference  In  The  average  computa- 
tion time  for  the  four  operations  are  less  chan  23%  of  the  minimum 
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average  computation  time.  Thus  it  seems  reasonable  to  manipulate  the 
expression  for  P(T)  to  minimize  the  total  number  of  operations.  Based 
on  the  above  criterion,  the  following  expression  for  p(T)  is  implement 
ed  in  the  predictor-correlator  of  the  microcomputer  system  for  target 
identification. 


i>(T) 


where 

e(t) 


The  number  of  arithmetic  operations  required  to  evaluate  p(T)  at  a T 
value  of  T . is  given  in  Table  16. 


VNT4V  , 

V e2 

t=t  +m  T 
s o 


e^(t) 


tcFf+m0T''  ' tc-NT+m0T 

)’  rfl(t)(rM(t)-e(t))  + )[  e (t) 

t +m„T  t=tc+mnT 


“ rM(t)-rc(t 
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TABLE  16 

NUMBER  OF  ARITHMETIC  OPERATIONS  REQUIRED 
TO  EVALUATE  p(T)  OF  EOUATION  (56)  AT 
THE  VALUE  OF  T=Toi 


OPERATION 

NUMBER  OF  OPERATIONS 

ADD 

(N+2)(M-NIo1)-l 

SUBTRACT 

M-NIQl 

MULTIPLY 

(N+l)(M-NIoi)+l 

DIVIDE 

1 

In  Table  16, 

N = number  of  resonances  of  the  desired  target  used 
for  Identification  purposes, 

M = (te-ts+Tg)/Tg,  and 
'oi  s To1 ’ 

Note  that  M-NI0*  Is  the  number  of  Instantaneous  error  samples  In 
the  error  i nterval . 

With  the  small  antenna  system  and  for  mine  Identification,  N=10, 
M=101,  loi=4,5,6,7  and  8 (see  Chapter  VII).  The  total  computation  time 
ranges  from  0.26  seconds  to  0.64  seconds  with  an  average  of  0.45  seconds. 


d.  The  filtering  operation 
1 n The  preprocessor 

Because  of  the  unavailability  of  the  FFT  package  in  the  micro- 
computer system,  the  digital  filtering  operation  In  the  preprocessor 
was  performed  in  the  time  domain  via  the  following  relationship^, 56] 

N' 

rM(nTB)  » E h(1TB)ri(nTB-1TB)  (57) 


where 
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r^(nTg)  = the  Processed  filtered  waveform, 

r 1 ( nig ) = the  Processed  unfiltered  waveform, 
h(1Tg)  = impulse  response  of  the  digital  filter,  and 
N‘  = order  of  the  digital  filter. 

The  impulse  response  of  the  digital  bandpass  filter  given  In 
Figure  47  was  obtained  by  taking  the  inverse  FFT  of  the  digital  filter 
transfer  function.  This  impulse  response  and  transfer  function  are 
shown  in  Figures  53a  and  b,  respectively.  Using  this  time-domain  fil- 
ter, a set  of  identification  results  was  obtained  using  the  In-house 
computer  for  the  Identification  of  the  mine-like  target  based  on  the 
small-antenna  waveforms  previously  obtained  (see  Chapter  VII).  Identi- 
cal estimates  of  identification  performance  were  obtained  (l.e.,  P,a 
100%,  Pp^=l . 7 2%  for  R j p=30  cm,  and  P j*l 00% , Pp^=6.90i  for  RjD”45  cm). 

The  time-domain  filter  shown  In  Figure  53b  Is  a FIR  filter  with 
a large  order  of  N‘*128.  Thus,  tremendous  computation  time  will  be 
consumed  in  the  filtering  operation  if  it  is  to  be  impelemented  in  the 
microcomputer  system.  To  save  computation  time,  the  order  of  this  FIR 
filter  was  lowered  to  37  by  using  a "windowing"  technique  with  a Kaiser 
window[55,56,61].  The  Kaiser  window,  the  "windowed"  impulse  response 
of  the  filter  and  the  FFT  of  the  windowed  impulse  response  are  shown 
in  Figure  53c-e.  Identification  performance  estimates  obtained  using 
the  In-house  computer  with  this  windowed  filter  was  Pt=10(M,  Pp/\=8 . 62% 
for  Riq-45  cm.  Thus,  performance  only  slightly  degraded.  This  filter 
was  Incorporated  into  the  microcomputer  system  for  real-time  identifi- 
cation of  the  mine-like  target.  The  windowed  time-domain  filter 
has  37  coefficients  and  is  symmetric  about  its  center.  Based  on 
Equation  (57)  and  Table  15,  the  average  time  required  for  the  filtering 
operation  on  a smal 1 -antenna  waveform  is  0.5  seconds  (vs  2 seconds  for 
the  unwindowed  filter).  The  digital  filter  can  be  replaced  by  an  ana- 
log filter  placed  at  the  front  end  of  the  receiver  in  the  final  identi- 
fication system. 

e.  Detection 

Operations  performed  in  the  detection  include  peak  detection, 
energy  calculations,  etc. 

f.  Execution  of  computer  Instructions 

The  amount  of  time  taken  for  execution  of  the  computer  instruc- 
tions depends  on  the  efficiency  of  the  microprogram  written. 

In  summarizing,  the  amount  of  time  required  for  an  Identification 
decision  in  the  Identification  of  the  mine-like  target  with  the  exist- 
ing smal 1 -antenna  microcomputer  system  is  estimated  to  be  approximately 
11  seconds  with  0.4E)  seconds  of  average  predictor-correlating  time. 


BANDPASS  FILTER 
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The  decision  time  goes  down  drastically  to  about  0,7  seconds  if  the 
repetition  frequency  of  the  pul  ter  can  be  increased"  to  1'6  <Hz,  and  the 
filtering  operation  In  the  preprocessor  can  be  performed  with  an  ana- 
log filter.  Using  state-of-the-art  technology,  both  these  conditions 
can  be  easily  met. 

Based  on  the  simulated  identification  results  obtained  with  the 
microcomputer  conditions  (l.e.,  128  samples  per  waveform,  FIR  filter- 
ing), the  microcomputer  identification  system  is  expected  to  identify 
the  mine-like  target  with  the  same  level  of  performance  as  the  results 
previously  presented  in  this  dissertation. 


D . CaJJ bratlon  of  the  M i c rocom outer 

fa  rg'e  t'~I(J~en"t  IT I catforTSystem 

Calibration  of  che  microcomputer  system  Includes  the  following: 

1 . Time  scale 

The  time  scale  was  calibrated  by  comparing  two  waveforms  from  the 
same  target.  One  of  the  waveforms  was  obtained  with  a fixed  length  of 
cable  added  to  the  system.  The  time  window  was  calibrated  at  25  ns. 
However,  the  impedance  mismatch  at  the  trigger  pick-off  circultory  of 
the  Terras~an  limited  the  reflectionless  time  window  to  20  ns  (see 
figure  04).  The  secondary  reflections  caused  by  this  Impedance  mis- 
match occurred  at  a time  which  was  20  ns  from  the  transmlt-recelve 
coupling.  This  reduction  In  the  size  of  the  reflectionless  time  window 
was  later  found  to  limit  the  Identification  range  of  the  radar  system. 

2 • Amp! i tu.de  scale 

Operations  performed  in  the  detector  required  accurate  calibration 
of  the  peak  amplitude,  peak  timing  as  well  as  the  energy  of  the  mine- 
liko  target  waveforms.  These  calibrations  were  performed  by  collecting 
a set  of  waveforms  aver  the  mine-like  target  and  tabulating  the  various 
detection  thresholds  based  on  the  set  of  waveforms 

A . On- location  ground  condition 

On-location  ground-condition  calibration  is  discussed  In  Chapter 
IX.  At  the  present  time,  the  in-situ  complex  resonances  of  the  mine- 
like  target  waveforms  were  obtained  by  analyzing  the  waveforms  collected 
in  ?.,  using  Prony's  method.  The  difference  equation  coefficients  asso- 
ciated with  these  resonances  were  then  used  for  real-time  target  iden- 
tification. Real-time  identification  performance  of  the  microcomputer 
system  is  discussed  In  the  next  section. 
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E . Real -T Ime  I dent  if lea t i on  Performance 

o f theJlTcroc oiifpuTe r~T arge  t 

Tcfen t i f fca tTb'n  System 

1 • Identification  performance 

Two  sets  of  identification  data  were  obtained  for  the  Identifi- 
cation of  the  mine-like  target  in  two  different  ground  conditions. 

The  first  set  was  obtained  using  the  poles  and  thresholds  as  given  in 
Tables  27  and  30  of  Appendix  G for  real-time  identification  of  the 
mine-like  target.  Real-time  identification  performance  of  Pi“l 00%  and 
PpA»0fl  was  obtained  for  Rig*30  cm  (based  on  9 mine-like  target  wave- 
forms and  21  other-target  waveforms).  It  was  found  that,  for  Ride45  cm 
identification  performance  degraded  to  Ptj12'/,,  Pp/\=0'/;,  This  reduction 
In  identification  range  was  attributed  ti  the  presence  of  the  unwanted 
reflections  from  the  trigger  pick-off  in  the  latter  portion  of  the  time 
window. 

I 

The  second  set  of  Identification  data  was  obtained  under  a dlf-  a 

ferent  and  rapidly  changing  ground  condition.  A set  of  measurements 
over  the  mine-like  target  was  collected  over  a lightly  wet  ground  for 
calibrating  the  ground  condition,  target  resonances,  and  the  various 
detection  and  Identification  thresholds.  However,  there  were  rapid 
weather  changes  between  the  time  of  data  analysis  (for  calibrations) 
and  real-time  identification.  There  was  heavy  rainfall  during  the 
interim  period,  and  In  one  day,  even  a moderate  snow  storm.  These 
drastic  weather  changes  altered  the  ground  condition,  and  with  the 
system  tuned  to  the  previously  calibrated  conditions,  It  was  found  that 
the  level  of  the  correlation  coefficients  for  the  mine-like  target  wave 
forms  was  generally  lowered.  Hence  to  ensure  a 1 00%  identification 
probability,  the  Identification  thresholds  had  to  be  lowered.  Thus, 
with  an  untuned  system  and  using  lower  threshold  values  (see  Appendix 
1 for  details),  a set  of  real-time  identification  data  was  obtained. 

Identification  data  based  on  17  mine-like  target  and  13  other-target 
waveforms  were  PplGOT  and  Pp/\-07.  for  Rjd-30  cm.  Thus,  lowering  the 
threshold  values  seems  to  be  an  effective  means  to  counter  the  problem 
of  uncertain  ground  conditions  for  this  penetrable  "desired  target".  , 

This  aspect  of  the  system  warrants  a detailed  future  investigation.  j 

tj 

2 . Speed  | 

The  two  sets  of  identification  d*ta  discussed  1 r.  the  above  sub- 
section were  obtained  with  N-|  (the  number  of  raw  waveforms  required  to  j 

form  an  average  waveform)  set  equal  to  8.  Thus,  for  a pulse  repetition  j 

rate  of  286  Hz,  the  amount  of  time  required  for  the  microcomputer  sys- 
tem to  form  an  average  difference  waveform  was  8 seconds.  The  amount  ) 

of  time  taken  for  the  correct  identification  of  tho  mine-like  target  ! 

was  clocked  at  approximately  2.5  seconds,  Thus,  the  total  real-time  ! 

for  a correct  identification  of  the  mine-like  target  was  approximate-  \ 

ly  10.5  seconds.  The  amount  of  time  required  for  fal  se-ta.  get. 
or  no-target  discrimination  was  less  than  10.5  seconds.  In  the 
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implementation  of  later  generations  of  this  identification  system,  it 
is  expected  that,  the  identification  speed  will  drastically  increase 
with  the  increase  in  pulse  repetition  fate,  decrease  In  N] , elimina- 
tion of  the  90°  difference  operation,  and  more  efficient  microprograms. 

3 . Precis! on 

To  find  the  precision  of  the  correlation  coefficients  calculated 
in  the  microcomputer  system,  correlation  coefficients  were  calculated 
for  a waveform  using  the  microcomputer  system.  These  correlation 
coefficients  were  then  compared  to  those  calculated  by  the  in-house 
computer  system  for  the  same  waveform.  It  was  found  that,  the  two  sets 
of  results  differed  In  the  10th  bit  of  the  mantissa  of  the  floating- 
point formatted  numbers.  No  error  was  observed  in  the  exponents.  Thus, 
using  the  25-bit  floating-point  format  of  the  API!  resulted  in  a G.l« 
truncated  error  In  the  correlation  coefficients.  This  is  expected  to 
have  little  effect  on  the  identification  performance. 


F.  Possible  Future  Improvements  on 

the  Mi  croc onlpu  ter  Sy stem 

In  addition  to  the  impedance  mismatch  at  the  trigger  pick-off 
circultory  of  the  Terrascan,  a few  other  aspects  of  the  microcomputer 
system  can  be  modified  to  improve  identification  performance. 

1 ■ The  sampling  system  in  the  Terrascan 

The  Terrascan  system  was  originally  designed  for  low-frequency 
operation  (i.e.,  frequencies  less  than  200  MHz),  as  such  its  sampling 
system  is  therefore  slow.  The  slow  sampling  system  causes  loss  in  high 
frequency  content  of  the  received  waveforms,  Figure  55  shows  two  wave- 
forms, one  received  by  the  Terrascan  and  the  other  by  a faster  sampling 
system,  From  the  same  target  at  the  same  antenna  location  and  orienta- 
tion. The  loss  of  high-frequency  content  in  the  waveform  received  by 
the  Terrascan  system  Is  quite  apparent. 

With  the  Identification  performance  obtained  thus  far  using  the 
microcomputer  system,  the  high-frequency  loss  In  the  Terrascan  sampling 
system  does  not  seem  to  affect  identifi cation  performance  for  the 
identification  of  the  mlnc-like  target  considered  in  this  study.  How- 
ever, such  hiyh-fr* quency  loss  will  almost  surely  degrade  performance 
when  the  system  is  used  to  Identify  smaller  mines. 

2 • The  dynam i c range  in  signal  level 

The  microcomputer  system  currently  uses  an  8-bit  A/D  converter 
to  convert  the  received  analog  waveforms  to  digital  form.  This  results 
in  a 48  dB  dynamic  range  in  signal  level.  For  the  identification  of 
low-level  signals  the  dynamic  range  needs  to  be  increased. 


MICROCOMPUTER  SCALE 


Fiyure  Db,  Processed  mine-like  target  waveforms  received 
by  the  microcomputer  and  the  fast  sampling 
system. 


FAST  SAMPLER  SCALE  ( mv  ) 
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3 . Gain  control  of  the  Terrascan 

The  gain  control  of  the  Terrascan  system  offers  a limited 
range  of  gains  in  large  discrete  steps.  Such  large  discrete  steps 
may  result  in  large  truncation  error. 

In  summarizing,  a microcomputer  system  has  been  implemented  for 
real-time  subsurface  target  identification.  Real-time  single-look 
Identification  performance  of  Pj=>100%  and  PpAs0 % for  RjD“30  cm  was 
obtained  with  a total  of  26  mine-like  target  waveforms  and  34  other- 
target  waveforms  for  the  identification  of  the  mine-like  target  in 
different  ground  conditions.  The  amount  of  time  required  for  a correct 
identification  of  the  mine-like  target  was  approximately  10.5  seconds. 
The  amount  of  time  required  for  a correct  other-target  or  no-target 
discrimination  was  less  than  10.5  seconds. 

The  microcomputer  system  built  in  -his  study  was  intended  to  be 
a "first-generation"  system.  Many  possible  modifications  of  the  system 
were  recommended  in  this  chapter  for  better  identification  speed,  data 
precision,  and  identification  performance.  The  most  important  modifi- 
cation Is  the  process  of  automatic  tuning  of  the  radar  system  to  the 
on-location  ground  condition.  This  is  the  topic  of  Chapter  IX. 
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CHAPTER  IX 

A METHOD  FOR  REAL-TIME  ON -LOCATION  TUNING  OF  1 HE 
IDENTIFICATION  RADAR  TO  THE  GROUND  CONDITION 


r. 
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A.  Objectives 

The  importance  of  tuning  the  radar  system  to  the  ground  condition 
was  illustrated  in  Chapter  IV.  An  un-tuned  subsurface  identification 
raaar  system  yielded  degraded  performance.  This  chapter  discusses 
a possible  method  for  automatic  tuning  of  the  radar  system  to  the 
ground  condition.  Emphasis  is  on  the  simplicity  of  tin-  method  and  the 
possibility  of  incorporatin'!  such  method  into  the  microcomputer  system 
for  on-location  target  identification  in  real  time. 


B.  The  Use  of  Backscattered  Waveforms 

FronTTh'in  "Wires  to  Estimate  Tne 

Ground  Parameters 

The  backscattered  waveforms  from  the  buried  thin  wires  as  shown 
in  Chapter  VII  possess  the  following  properties  which  are  useful  in 
the  calibration  of  ground  conditions: 

1.  The  waveforms  contain  complex  natural  resonances  which 
are  related  to  the  ground  parameters  in  a simple  fashion. 

?.  1 he  wire  waveforms  arc  of  high  amplitude,  thus  the 

estimate  of  the  resonances  will  be  accurate. 

3.  The  lowest-order  target  resonance  of  the  thin  wires  are 
the  most  dominant.  (See  Appeidix  J).  This  implies  fli.it 
t.he  thin  wire  waveforms  can  be  closely  characterized  o> 
a single  resonance.  Thus, 

“1 C t-t‘ ) 

rM(t)  = 2a|  e cosi,^  ( t - 1 ' ) (58) 

where 

r^(t)  - processed  waveform  , 
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S1  = L'li+J|a'i  " complex  resonance, 

=>  residue  associated  witn  c.,  and 

t1  - start  time  of  the  transient  signal  from  the  wire. 

For  a thin  wire,  the  first  resonance  can  be  approxlniuted  by 
Equation  (58) 


3 - y-  - 0.0828  -~£-. 
1 K 


"1 

where 


0.9251 


u = ground  conductivity  (homogeneous  medium  assumed), 
» ground  permittivity  , 
c ■ speed  of  light  in  free  space  , 
r = relative  dielectric  constant  of  the  ground,,  and 
L = length  of  wire  . 


(59) 


The  complex  natural  resonance  given  by  Equation  (59)  is  a good 
approximation  to  the  first  natural  mode  of  a thin  wire  with  a length/ 
diameter  ratio  of  100  buried  in  a homogeneous  medium  of  good  dielec- 
tric (i.e.,  medium  with  a small  loss  tangent  of  o/m <<] ) [52 ,69] . 

Here,  it  is  used  to  estimate  the  ground  parameters  based  on  the  wave- 
forms from  the  buried  thin  wires  and  is  found  to  yield  reasonable 
results. 


Since  the  waveform  is  dominated  by  a single  resonance,  the 
values  of  .u  and  i,q  can  be  easily  estimated  from  the  decaying  envelope 
and  the  peak  timing  at  two  samples  of  the  waveform,  thus,  the  ground 
parameters  can  be  estimated  via  the  following  equations: 


0.9251 ( t1 -tq)C 
2L 


i ( t1  — 10)_1  tn 


- 0.0828 


and 


2 
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where  r(t„)  and  r(t,)  are  the  sample  value,  of  the  measured  waveform 
at  time  tg  and  t^ , respectively.*  All  other  parameters  are  as  pre- 
viously defined. 

Using  Equation  (60),  the  ground  parameters  represented  by  the 
waveform  shown  in  Figure  56  are  estimated  to  be  r.  =14  and  o=24  mS/m.** 
These  estimates  are  found  to  be  very  close  to  therestimates  obtained 
suing  different  techniques  currently  being  investigated  in  another 
study  at  the  El ectroSc ierce  Laboratory!! 58  J.  The  method  outlined  In 
Equation  (60)  is  attractive,  for  it  uniquely  characterizes  the  ground 
condition  by  two  sample  values  of  a thin-wire  waveform. 

C.  Automatic  Tuning  of  the  Identification 

HadaVTo,  the  Ground  Condi  ti  or7 

The  knowledge  of  the  ground  parameters  basically  solves  the 
automatic  ground-condition  tuning  problem  for  the  real-time  identifi- 
cation of  simple  targets  such  as  the  brass  cylinder  and  the  aluminum 
sphere  whose  resonances  are  related  to  the  ground  parameters  In  a 
known  analytical  fashion[21 ,63].  However,  it  does  not  solve  the  tun- 
ing problem  for  the  real-time  identification  of  plastic  mines,  for 
the  analytic  relations  between  the  complex  natural  resonances  of  the 
plastic  mine-like  target  and  the  ground  parameters  are  not  known. 

To  date,  the  characterization  of  t.rgets,  should  it  be  free-space  or 
subsurface  targets.  Is  a relatively  new  research  problem.  Further- 
more, research  efforts  have  mostly  been  concentrated  in  the  character- 
ization of  perfectly  conducting  targets  only.  Characterization  of 
dielectric  targets  such  as  the  plastic  mine-like  target  with  complex 
natural  resonances  «ia  analytical  method  has  not  been  previously 
treated.  This  could  be  achieved  but  it  represents  a substantial 
research  effort  beyond  the  scope  of  the  present  study.  In  this  study 
we  suggest,  solving  the  automatic  tuning  problem  for  identification 
of  mines  experimentally  with  the  procedure  outlined  as  follows: 

Step  1:  Estimate  the  ground  parameter  by  using  the  method 
outlined  by  Equation  (60). 


★ 


Equation  (GO) 
zero  and  that: 


is  val id  only  when  r 
( 1 1 - tQ ) = 2r»n  where 


are  both  non- 


**  Note  that  these  estimates  are,  of  course,  accurate  only  at  t.i, . 
Furthermore,  the  loss  tangent,  based  on  the  estimated  values  of 

i and  o,  is  found  to  be  0.252. 
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Step  2:  Determine  the  loci  of  the  complex  resonances  of  the 
mine-like  target  experimentally  as  the  ground  para- 
meters vary.  Ground  parameters  can  be  controlled 
by  using  different  kinds  of  soils,  or  by  Increasing 
the  moisture  and  salt  content  of  the  ground.  Thus, 
the  loci  of  the  natural  resonances  (or  the  coeffici- 
ents Ufp 1 s ) can  be  determined  as  a function  of  t^-tg, 
tmUo),  and  r^ti)  and  stored  in  ROM's  In  the  micro- 
computer identification  system  for  on-location  target 
identification  in  real  time. 

To  make  the  above  method  even  more  attractive  and  practical , the 
thin  wire  does  not  have  to  be  buried  for  ground-condition  characteri- 
zation. Figure  57  shows  waveforms  obtained  ,from  a thin  wire  laying 
at  different  locations  on  the  surface  of  the  ground.*  Again,  the  wave- 
form Is  dominated  by  one  (pair)  resonance  and  thus,  the  ground  condition 
is  uniquely  characterized  by  the  values  and  timing  of  two  sample  points 
of  the  waveform.  In  the  case  of  the  on-surface  wire  however,  because 
of  the  fact  that  part  of  the  wire  Is  in  the  air,  the  ground  parameters 
are  no  longer  estimated  accurately  by  Equation  (60).  Nevertheless,  the 
unique  Information  about  the  ground  condition  Is  assumed  to  be  contain- 
ed in  the  two  samples  of  the  waveform.  Thus  far,  this  assumption  is 
found  to  be  valid.  The  on-surface  wire  waveforms  given  In  Figure  55 
Indicate  different  decaying  envelopes  and  zero-crossing  Intervals  for 
different  ground  conditions. 

The  method  of  ground-condition  tuning  outlined  above  amounts  to 
an  extensive  experimental  effort  and  is  planned  for  the  future. 


* The  wire  has  to  be  In  good  contact  with  the  ground. 


CHAPTER  X 

SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  Summary 

At  the  time  research  In  this  report  was  Initiated,  the 
claim  that  electromagnetic  signals  were  impractical  for  subsurface 
exploration  because  of  attenuation  had  already  been  refuted.  Refer- 
ences In  this  report  contain  additional  evidence  refuting  this  claim. 

We  demonstrate  that  the  received  signals  from  objects  beneath  the 
surface  of  the  earth  can  also  be  interpreted  Intelligently.  That  Is, 
that  much  more  than  simple  ray  tracing  with  a postulated  reflection 
coefficient,  at.  some  depth  is  possible.  We  show  conclusively,  using 
real  subsurface  probing  data,  that  the  target  can  be  identified. 

It  is  emphasized  that  the  Identification  is  not  simply  a seismic-type 
map  subject  to  Interpretation  and  quail ficatTorf  by  the  observer  but 
an  actual  processing  scheme  which  inquires  as  to  the  presence  of  a 
particular  target.  We  also  stress  that  the  methods  developed  In  this 
dissertation  are  field-oriented  and  operate  In  real  time.  The  para- 
graphs below  Itemize  the  specific  progress  which  has  been  made. 

The  predictor-correlator  method  for  characterizing  and  identi- 
fying subsurface  targets  was  studied  and  extensively  tested  using  real 
radar  measurements  obtained  with  a Terrascan-type  subsurface  pulse 
radar.  Measurements  were  obtained  over  three  sets  of  targets.  The 
first  set  consisted  of  five  similar-size  targets  Including  a plastic 
mine-like  target,  a brass  cylinder,  an  aluminum  sphere,  a copper 
sheet  and  a wood  board.  These  targets  were  burled  at  a depth  of  5 
cm  (2  Inches).  The  second  set  consisted  of  a series  of  different-size 
(maximum  dimension  varies  from  30  cm  to  300  cm)  brass  cylinders  buried 
at  different  depths  (depth  varies  from  30  cm  to  150  cm).  The  third 
set  consisted  of  a series  of  thin  wires  burled  at  the  depth  of  5 cm 
(2  Inches).  The  method  of  characterization  and  identification  uses  the 
complex  natural  resonances  as  the  discriminants  for  the  identifier  and 
the  method  of  linear  prediction  for  evaluation  of  the  correlation  coef- 
ficients for  threshold  Identification.  The  complex  natural  resonances 
of  the  desired  target  are  determined  a priori  and  the  difference  equa- 
tion coefficients  related  to  these  resonances  are  stored  in  the  Identi- 
fication system  for  real-tine  on-location  target  identification.  The 
Identification  process  is  simple  and  Involves  only  simple  algebraic 
operations . 
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The  characterization  and  Identification  methods  were  found  to  be 
successful  and  a "first-generation"  microcomputer  system  was  built  for 
on-location  Identification  of  mines  In  real  time.  A simple  method  for 
automatic  tuning  of  the  Identification  radar  to  the  ground  condition 
was  suggested.  This  method  Is  simple  to  use  and  can  easily  be  Incorpor- 
ated Into  the  microcomputer  system  for  real-time  target  Identification 
purposes.  The  significant  findings  are: 

1.  The  identification  method  was  extensively  tested,  with 
measurements  over  the  subsurface  targets,  and  found  to 
be  extremely  successful.  Single-look  identification 

s£a.t.1sjLi<;s  fo.r,  the, Identification  of  the  mine-like  target 

In  different  ground  conditions  were  est'l mated  t*o  be'  * ' 
Pt-100%,  PcA-2.25iu  for  R];qb30  cm  over  an  ensemble  of  55 
mine-like  target  waveforms  and  222  other-target  waveforms.* 
The  S/C  of  the  ensemble  of  mine-like  target  waveforms 
ranged  from  0.21  to  3.5. 

2.  A microcomputer  system  was  Implemented  for  real-time  sub- 
surface target  identification  using  the  techniques  develop- 
ed In  this  study.  Single-look  Identification  statistics 
for  the  identification  of  the  mine-like  target  were  esti- 
mated to  be  Pt“100%,  Pp/\*0%  for  Rid“30  cm  over  an  emsenible 
of  30  mine-like  target  waveforms  and  30  other- target  wave- 
forms. The  amount  of  time  required  for  a correct  identifi- 
cation of  the  mine- l Ike  target  was  10.5  seconds,  a correct 
discrimination  of  an  other-target  or  no-target  required 
less  time. 

3.  Identification  performance  degraded  when  the  radar  system 
was  not  tuned  to  the  right  ground  condition 

4.  Identification  performance  degraded  when  the  radar  fre- 
quencies did  not  properly  span  the  target  resonances. 

5.  Single-look  Identification  performance  was  characterized 
by  a correlation  coefficient  vs.  sampling  interval  ( p (T ) ) 
curve.  Optimal  Identification  performance  occurred  when 
the  sampling  Interval  T was  In  the  Immediate  neighborhood 
of  Tn  (see  Equation  (51)).  Thus,  one  knows  a priori  the 
value  or  the  region  of  values  of  the  sampling  interval 
needed  for  design  and  implementation  of  the  Identification 
radar. 


* Identification  statistics  given  Is  estimated  from  the  ensemble 
of  all  mine-like  target  waveforms  collected  using  the  Terrascan- 
llke  systems. 
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6.  All  subsurface  targets  considered  were  characterized  by 

a small  and  finite  number  (5  pairs  or  less)  of  resonances. 
These  resonances  were  Invariant  with  respect  to  radar 
location. 

7.  The  resonances  of  the  plastic  mine-like  target  were  found 
to  be  Internal  resonances,  with  their  Imaginary  parts 
Independent  of  ground  condition.  The  resonances  of  the 
brass  cylinder  were  found  to  be  external  resonances; 
both  their  Imaginary  and  real  parts  were  dependent  on 
ground  conditions. 

8.  The  extracted  resonances  of  the  brass  cylinders  were 
found  t.o  be  the  dipoles  modes  along  the  length  of  the 
cylinder.  Furthermore,  they  closely  approximated  the 
resonances  of  the  same  target  burled  In  a homogeneous 
medium  with  plane-wave  Illumination.  This  discovery 
suggested  a good  way  to  estimate  the  ground  parameters 
(see  Section  VI1-B), 

9.  The  hlcjh-frequc-icy  content  of  the  backscattered  waveforms 
was  highly  attenuated  as  target  depth  Increased,  Identi- 
fication range  decreased  as  target  depth  Increased. 

10.  Up  to  a certain  threshold,  the  frequency  of  the  target 
resonances  decreased  according  t,o  the  Increase  In  target 
size,  and  target  resonances  depended  solely  on  the  scatter- 
ing mechanisms  of  the  targets.  Beyond  the  threshold,  In- 
crease In  target  size  did  not  warrant  the  decrease  In  the 
frequency  of  the  resonances.  In  this  case,  target  reson- 
ances would  depend  on  the  scattering  mechanisms  as  well  as 
other  quantities  such  as  the  antenna  pattern,  etc. 

The  significant,  contributions  of  this  report,  are: 

1,  The  problem  of  applying  Prony's  Method  to  real  radar  measure- 
ments for  extraction  of  the  complex  resonances  was  consider- 
ed as  one  of  parameter  optimization.  This  approach  yielded 
legitimate  target  resonances  from  waveforms  with  S/C  as  low 
as  0.P1  (see  Chapter  III). 

d.  The  predictor-correlator  method  was  extensively  tested  with 
real  radar  measurements,  and  found  tc  yield  practical 
single-look  Identification  performance.  Furthermore,  It 
was  found  that  there  existed  a limited  range  of  T values  in 
which  optimum  identification  performance  occurred.  (See 
Chapter  IV). 


■ :%i  •:  ,N*tj : lil 


142 


3.  A microcomputer  system  was  implemented  for  real-time  on- 
location  subsurface  target  identification  using  the  tech- 
niques developed  in  this  dissertation.  The  system  was 
found  to  yield  practical  single-look  identification  per- 
formance with  a 10.5  second  identification  time  (see 
Chapter  IX). 


B.  Conclusions 

We  have  now  provided  a successful  method  for  the  identification 
of  subsurface  targets  at  shallow  depth.  The  method  is  based  on  single 
radar  returns  and  is  simple  to  implement^  Based  on  this  method  a micro- 
computer system  has  been  implemented  for  on-location  identification  of 
subsurface  targets  in  real  time.  Some  modifications  have  to  be  made 
for' the  microcomputer  identification  radar  system  to  be  practical,  the 
most  important  being  a method  to  calibrate  the  ground  parameters  (e,u, 
o)  in  real  time  and  to  adjust  the  resonances  of  the  desired  target 
accordingly  for  real-time  subsurface  target  identification  in  different 
ground  conditions.  This  is  now  being  implemented  and  is  expected  to  be 
incorporated  in  the  final  system. 


C . Recommendations  for  Future  Work 

The  subject  of  subsurface  target  characterization  and  identifica- 
tion is  in  its  infancy.  The  results  obtained  in  this  study  represents 
a significant  step  forward,  but  much  remains  to  be  done.  The  following 
items  are  highly  recommended  for  future  research: 

1.  The  taks  of  automatic  on-location  tuning  of  the  subsurface 
radar  to  the  ground  condition  is  crucial  to  the  problem  of 
subsurface  target  identification.  The  consideration  of  the 
waveform  values  of  the  on-surface  wire  waveforms  suggests 

a feasible  way  to  solve  this  problem. 

2.  The  relationships  between  the  extracted  resonances  of  the 
targets  and  the  various  system  parameters  such  as  target 
size,  depth,  antenna  pattern,  air-ground  interface,  etc., 
need  to  be  exploited  for  a better  understanding  and  inter- 
pretation of  the  extracted  resonances. 

3.  Application  of  the  characterization  and  identification 
technique  to  deeper  targets  needs  to  be  expanded.  In  this 
regard  very  little  deep  electromagnetic  probing  has  been 
done. 

4.  Identification  of  subsurface  targets  in  the  presence  of 
other  objects. 
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APPENDIX  A 

DERIVATION  OF  PRONY'S  METHOD  FOR  TRANSIENT  WAVEFORMS 
WITH  MULTIPLE-ORDER  POLES 


This  appendix  derives  Prony's  method  for  the  extraction  of 
resonances  from  waveforms  with  multiple-order  poles. 

With  multiple-order  poles  the  transient  waveform  can  be 
expressed  as 

N f '^1  i ,\  Sjt 

Kt)  -l  1 + l b iPtJ*'  a.e  1 (( 

i-1  V j=2  / 1 

where 

p - 0,  if  Mi  * 2 
p * 1,  if  M1  > 2 

and  where  Mi  is  the  multiplicity  of  the  ith  pole. 

The  corresponding  expression  of  Equation  (61)  in  the  complex 
frequency  domain  is 
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In  discrete  form,  Equation  (61)  can  be  written  as 
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With  the  following  2-transform  pairs [54 3 
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where  ?,[]  Is  the  z transform  operation.  We  can  transform  Equation 
(63)  into  the  z domain,  viz., 
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Each  term  in  the  right  hand  side  of  Equation  (65)  is  a rational 
function  of  z"’  with  its  denominator  one  degree  higher  in  z'1  than  the 
numerator.  Thus,  when  the  series  is  summed,  it  becomes  a rational 
function  of  z_1  in  the  form  of 
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where  the  polynomial  of  z"  in  the  denominator  is  one  degree  higher 
than  that  in  the  numerator.  The  denominator  polynomial  is  commonly 
referred  to  as  the  characteristic  polynomial  and  its  related  to  the 
pole  locations  s^  by 
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Multiply  both  sides  of  Equation  (67)  by  z*-,  we  obtain 
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where 
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Equation  (68)  is  analogous  to  Equation  (11). 

We  now  proceed  to  write  Equation  (66)  as 

dQ+d1z"1  + dLz"Lk[r(nT)]  = Cg+^z"1  + •••  CL_-j  z“L+1  .(70) 


With  the  relationship[54] 

■•■'1[2'J(r(nT))]  = r(nT-IT)  (71) 

where  ] is  the  inverse  i transform  operation.  We  obtain  the 

inverse  z transform  of  Equation  (70) 
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Since  the  right-hand  side  of  Equation  (72)  is  zero  for  n^L.  hence 
L 

l dj r(nT-IT)  = 0 ; n > L . (73) 

1=0 

Now  with 

n - L + K;  K - 0,1,2  — 
m = L-i  and 
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we  obtain  the  desired  difference  equation, 
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where  r^..--r(KT+niT).  This  difference  equation  is  identical  to  the 
Prony  difference  equation  in  Equation  (12)  except  for  the  change  of 
order.  Thus,  the  sample  values  of  a transient  waveform  with  the  pre- 
sence of  multiple-order  poles  satisfies  a Prony  difference  equation 
of  order  L,  where 

N 

L •-  V M1  . (75) 

1*1 


With  Equation  (74)  wo  can  solve  for  the  coefficients  u^'s,  and  sub- 
sequently solve  for  the  s^  by  Equation  (67),  We  can  also  solve  for 
the  poles  5i,  as  was  done  in  the  simple-pole  case,  by  solving  Equations 
(12)  and  (11)  with  N changed  to  L.  Thus,  the  same  procedure  taken  to 
solve  for  the  polo  locations  in  the  simple-pole  case  can  also  be  used 
in  the  multiplo-order-pole  case.  However  the  procedure  for  the  calcu- 
lation of  the  residues  requires  a slight  modification. 

The  calculation  of  the  residues  is  done  by  solving  Equation  (63), 
which  differs  from  Equation  (H)  because  of  the  presence  of  the  terms 
Involving  (nT),  and  the  fact,  that  there  are  L unknowns  rather  than  N. 
With  tiie  assumption  that  the  ith  pole  is  of  order  , Equation  (63) 
can  be  written  as 
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l-rom  liquation  ( 75} , we  see  that,  In  solving  for  the  residues  In 
the  presence  of  multiple-order  poles,  the  matrix  Equation  (15  must 
be  modified  as  follows 
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APPENDIX  R 

DERIVATION  OF  PRONV'S  METHOD  AND  ITS  VARIATIONS 


This  appendix  derives  the  variations  of  Prony's  method  aue  to 
the  constraints  u-|=l  (Interpolation  method)  and 


(Eigenvalue  method).  Computer  codes  (in  Fortran  Language)  for  these 
methods  are  also  included. 


A • The.  Interpolation  Method 

In  the  Interpolation  method,  the  Prony  difference  equation  can 
be  written  as 


“or < t ) +.«i r ( t+T ) + - ••a1_1r(t+(1-l)T)  + "-u1  + 1r(t+(1+1)T)  + "> 

+ ,xNr(t+NT)  » - air(t+il)  . (76) 


Thus  the  matrices  A,  B and  C shown  in  Equation  (13)  need  to  be 
modified  as 
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liquation  (77)  can  bn  used  to  solve  for  the  coefficients  u ,(!,,••• 
a-j  _ | .1  j + 1 , • * • in  the  interpolation  method.  Once  the  coefficients 

.ire  ftnown,  we  can  solve  for  the  poles  and  residues  by  following  the 
Identical  procedure  outlined  in  the  Classical  Method  (Chapter  III), 


13.  The  eigenvalue  Method 
Under  the  constraint 


the  instantaneous  error  in  the  error-evaluation  process  of  Prony  is 
modified  follows 

N N , 

e(thNT)  - V '»  r (t+mT)  ; V •/  = ! . (78) 
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In  matrix  form,  the  instantaneous  error  can  be  written  as 
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The  total  square  error  over  the  fitting  Interval  (ts,ts+MT-T)  is 


L1  = clTrl 
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where  l i is  the  total  squared  error  over  the  fitting  Interval.  [ ] 

denotes  the  transpose  operation.  T = A 1 A » is  the  data  covariance 
mati x [53] . 

Solving  for  Matrix  B by  minimizing  under  the  constraint 

N y 

\ «£•' 

ir.-O 

Is  a standard  eigenvalue  problem.  Via  the  eigen-analysis,  the 
solution  matrix  B is  the  eigen-vector  corresponding  to  the  minimum 
eigenvalue  of  the  covariance  Matirx  . 

Thus,  Equation  (80)  can  be  used  to  solve  for  the  coefficients 
, , Once  these  coefficients  arc  known  the  identical  proce- 

dure outlined  in  the  Classical  Prony's  method  can  be  used  to  solve 
for  the  poles  and  residues. 

Computer  codes  for  the  Prony's  method  and  its  variations  are 
given  below.  The  program  names  "SEMI"  implements  the  Prony's  method 
under  the  constraints  of  u,|(=i,  m~0, 1 ,2  • • • N (See  Chapter  III).  The 
program  " SemZ"  implements  the  Prony's  method  under  the  constraint  of 
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Both  programs  are  well  commented  and  user  oriented. 
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APPENDIX  C 


This  appendix  tabulates  the  extracted  resonances  of  the  mine- 
like target  at  different  antenna  locations  in  icy  ground.  The  poles, 
residues  as  well  as  the  minimum-error-case  parameters  associated 
with  the  Prony  Process  are  given.  These  parameters  are  N,  IBS,  III, 
and  M and  are  defined  as  follows. 

N - Number  of  poles 

IBS  = Interval  between  samples,  i.e.,  T=IBSxTg 

III  = The  start  time  of  the  fitting  interval  in  the 
Prony  process  is  ts+( 1 1 1-1 )x2xTg 

M = Number  of  samples  used  in  the  fitting  interval  is 
MxN . 


TABLE  17 

EXTRACTED  RESONANCES  OF  THE  MINE-LIKE  TARGET 
AT  VARIOUS  ANTENNA  LOCATIONS  IN  ICY  GROUND 
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ANTENNA  LOCATION  = CENTER 

N-l 1 , I8S=5 

, I I 1=5,  M=3 

— 

E = 0. 

632E-2 

POLE 

POLE 

RESIDUE 

RESIDUE 

(REAL  PART) 

(I MAG  PART) 

(REAL  PART) 

(IMAG  PART) 

-.6165323E8 

±.  5980379E8 

. 1627237  EO 

i.22C3193E0 

-.9816088E8 

+.1153422E9 

-. 1200819E0 

•i.  4662508E-1 

- . 2972145E9 

i.  31 99646E9 

-.6727021 EO 

±.2528001E-1 

- . 2574395E9 

■E.4525337E9 

-. 2245525E0 

i.9734140E-l 

15  cm  EAST 

OF  CENTER 

N=ll,  IBS=7 

i;=0. 

, I I 1=1,  Ma2 
604E-3 

POLE 

POLE 

RESIDUE 

RESIDUE 

(REAL  PART) 

(IMAG  PART) 

(REAL  PART) 

(IMAG  PART) 

- . 7955556E9 

+.7207162E8 

-.3441590E0 

?5773295E0 

- . 9747022E8 

+.1 32751 7E9 

- . 29 1 4887  E - 1 

i. 4185661 EO 

-.2185998E9 

1. 1964866E9 

. 17861 04 El 

6787741 EO 

- . 2235657E9 

+. 2570622E9 

.1210215E1 

±..  1496749E1 

- , 2108089E9 

+ . 3265433E9 

- . 291 4887E-1 

T,  41 85661 EO 

t i e..i  , it  ti  14  1 

t ll,  r ,la  4 i u.-,  r .*  n 

t-  (t  -Milt  V<i  .*  :i  •:  T r.tt  *• 

N=10,  IBS=4 

, 1 1 1 =2 , M=2 

15 .cm  SOUTH  OF  CENTER 

r=0. 

375E-2- 

POLE 

POLE 

RESIDUE 

RESIDUE 

(REAL  PART) 

(IMAG  PART) 

(REAL  PART) 

(IMAG  PART) 

-.4419169E8 

+.5302244E8 

-. 1283107E0 

±.  266R165E0 

1687482E9 

+.1862940E9 

. 1 938945E0 

.t.  4774533E0 

- . 1 866233E9 

+-.2797166E9 

. 2234024E-1 

±.  551 2929E0 

-.3106326E8 

-t.  3993629E9 

- . 5570095E- 1 

i..  481 6288E0 

-.4141 82 1 E9 

+.55 7 1886 EE 

. 7225376E0 

T-,  9451689E0 



i a.  ■ x.-i  . i 

* Real  and  Imaginary  parts  of  the  extracted  resonances  arc-  In 
Nepers/s  and  Hz,  respectively. 
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TABLE  17  (Cont.) 


15  cm  WEST 

OF  CENTER 

N=7,  1BS*9, 

111=3,  M=3  ■ | 

e*0.448E-2  i 

POLE 

(REAL  PART) 

POLE 

( 1MA6  PART) 

RESIDUE 
(REAL  PART) 

RESIDUE  j 

(I MAG  PART)  l 
| 

825269 IE8 
-.636610PE8 
- . 247429CE9 
- . 1 196715E9 

A.  71 59227EB 
±.  1213223E9 
±.  21 23833E9 
± . 2833333E9 

. 47037 1 OEO 
.2435532E-1 
•3750488E0 
. 1486354E9 

±.2550093E0  ' 
T.1370926EO 
3F.2747979E0 
0.0000000 

15  cm  NORTi 

OF  CENTER 

N=12,  IBS-6 
e=0. 

, 1 1 1=5,  M=2 

287  -3 

POLE 

(REAL  PART) 

POLE 

(I MAG  PART) 

RESIDUE 
(REAL  PART) 

RESIDUE 
(I MAG  PART) 

- . 1 238699E9 
- . 1 455328E9 
- . 2976452E9 
- . 358251 7E9 
- . 3250828E9 

±. 5956967E8 
+.1044725E9 
+ . 2543575E9 
±.  3319750E9 
+. 425000QE9 

7568201 EO 
.4427186E0 
. 1 651 730E0 
. 3949756EO 
, 2408422E0 

^.4133825E0 
+.3773904E0 
±.  41 33825E0 
7346425E0 
0.0000000 

30  cm  EAST  OF  CENTER 

N=1 1 , IBS=6 

c=0. 

, 111=4,  M«2 
247E-2 

POLE 

(REAL  PART) 

POLE 

(I MAG  PART) 

RESIDUE 
(REAL  PART) 

RESIDUE 
(I MAG  PART) 

-.9395545E8 
1369431 E9 
- . 3 1H7709E9 
-2052859E9 
- . 2686828E9 

t.  71 31887E8 
t.  1220007E9 
±. 2553445E9 
±.  3777737E9 
±.  4250000E9 

- . 3951 31 8E0 
-.1 24341 8E1 
1835029E1 
. 1061 165E0 
2203826E0 

t,  1237333E1 

+.8057936E-1 

T.6562829E0 

+.9037383E-1 

0.0000000 

30  cm  SOUTH 

OF  CENTER 

N=13>  IBS-6 

E = 0. 

, 111=2,  M=2 
406E-3 

POLE 

(REAL  PART) 

POLE 

(I MAG  PART) 

RESIDUE 
(REAL  PART) 

RESIDUE 
(I MAG  PART) 

- . 7440484E8 
- . 2734229E9 
- . 2451866E9 
' - . 2480318E9 
, - . 2653G34E9 

+. 5768431 F8 
i . 1749879E9 
i . 2385293E9 
i.  2889975E9 
+..3814169E9 

. 2634359E0 
- . 5841 965EO 
. 1 394630E1 
. 681 1801 EO 
- . 6790937E0 

±. 3522781 EO 
1580445E1 
1721972E0 
+. 1 354053E1 
±. 1 1 B7006E1 

9 


Sfifc; 
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TABLE  17  (Cont.) 


30  cm  WEST 

OF  CENTER 

N=13,  IBS=6,  II 1=3,  M*2 
u=0.685E-2 

POLE 

(REAL  PART) 

POLE 

(IMAG  PART) 

RESIDUE 
(REAL  PART) 

RESIDUE 
(IMAG  PART) 

-.11 74326E9 
- . 27221 84E9 
-.2442735E9 

+ . 9037581 E8 

1..  1783297E9 

1.. 3781 778E9 

- .7656702E0 
. 3367912EO 
-.2437025E0 

+.87I9310EC 

T.1380838E0 

±.4363209E-1 

30  cm  NORTH 

OF  CENTER 

N=12,  IBS=6 

c=0. 

, II 1=2,  M=2 
171E-2 

POLE 

(REAL  PART) 

POLE 

(IMAG  PART) 

RESIDUE 
(REAL  PART) 

RESIDUE 
(IMAG  PART) 

- . 39291  71  EES 
-.3953903E8 
- . 2863D06E9 
- . 21 30829E9 
- . 3845860E9 
- . 3482679E9 

±.6274667E8 
+ . 1 162181E9 
+. 1847860E9 
+.2626060E9 
t.  2977974E9 
±. 3921 698E9 

.3064080E0 
-.8984933E-1 
- . 2726371 ET 
1131884E1 
.2613163E1 
. 1288635E1 

T.2378655E-2 
+.4503131E-1 
+. 3413519E0 
T.7116251E0 
t. 2442230E1 
+.1489745E1 

AVEBME 

POLE 

(REAL  PART) 

POLE 

(IMAG  PART) 

- . 74931 16E8 
- . 9981 995E8 
- . 241 6503  L9 
; -. 2809195E9 
j 2885261 E9 

+ . 6347621 E8 

1. 1 146405E9 

1.  .’535799E9 
+.3074991 E9 
+.4076659E9 

APPENDIX  D 


i appendix  tabulates  the  average  extracted  resonances  of 
target  waveforms  obtained  using  the  12m  long  antenna. 


Table  18  J 

WE RAGE  EXTRACTED  RESONANCES  OF  THE  MIME-LIKE  1 

TARGET  WAVEFORMS  OBTAINED  USING  | 

THE  12m  LONG  ANTENNA  1 


POLE 

POLE 

(REAL)* 

(IMAG) 

-242. 9760E6 

* 144. 3435E6 

-381 . 1 344E6 

1 233. 9871 E6 

-407. 1770E6 

' 31 9. 0055E6 

-373. 7482E6 

>413. 1869E6 

*Reai  and  imaginary  parts  of  the 
extracted  resonances  are  in  Nepers/s 
and  Hz,  Respectively. 
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APPENDIX  E 


This  appendix  lists  the  Fortran  program  for  performing  the  target 
identification  process  (SRT1D2) , It  also  tabulates  some  of  the  target 
identification  results  described  in  Chapter  IV.  Detail  correlation 
coefficient  values  are  given. 
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COEFFICIENT  FOR  THE 
MINE-LIKE  TARGET 


TABLE  21 

VALUES  OF  THE  CORRELATION  COEFFICIENTS  FOR  THE 
IDENTIFICATION  OF  THE  MINE-LIKE  TARGET 
AND  THE  BRASS  CYLINDER  IN  DRY  GROUND 


WAVEFORM 

ANTENNA 

LOCATION 

DESIRED  TARGET 
= MINE-LIKE 
TARGET 
p(8Tb) 

DESIRED  TARGET 
= BRASS 
CYLINDER 
p(5Tb) 

C 

.525 

.648 

15  cm,E 

.434 

.646 

15  cm,S 

.423 

.491 

15  cni,W 

.477 

.651 

MINE-LIKE 

15  cm.N 

.837 

.922 

TARGET 

30  cm,E 

.687 

.738 

30  cm,S 

.535 

.168 

30  cm,W 

30  cm,N 

.396 

.831 

C 

-.266 

.954 

15  cm.SN 

.228 

.983 

BRASS 

30  cm,SN 

.146 

.951 

CYLINDFR 

15  cm.EW 

.228 

.979 

30  ciii.EW 

-.049 

.991 

IK 


TABLE  23 

DETERMINING  THE  DETECTION  AND  IDENTIFICATION  THRESHOLDS  FOR  THE 
IDENTIFICATION  OF  THE  BRASS  CYLINDER.  RID  = 30  cin 


r 
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KAJt  is  units  of  200  ■*¥ 


APPENDIX  F 


This  appendix  tabulates  the  average  extracted  resonances  of 
the  different-size,  different-depth  cylinders  and  the  thin  wires 
(discussed  in  Chapter  V). 


APPENDIX  G 


This  appendix  tabulates  the  average  extracted  resonances  of 
the  waveforms  from  the  two  models  of  the  mine-like  target  and  the 
small-antenna  mine-like  target  waveforms. 


TABLE  27 

AVERAGE  EXTRACTED  RESONANCES  OF  THE  MINE-LIKE  TARGET 


POLE 

(REAL)* 

POLE 

(IMAG) 

MINE-LIKE  TARGET 
MODEL  NO,  1 

0.6m  ANTENNA 

-1 . 75363970EB 
-5. 73125538E7 
-2. 87494683E8 
-1.82908509E8 
-9. 76454733E7 

±fi .610191 33E7 
i-1 .321 22486E8 
+2.27030133E8 
+3.0741 1043E8 
+4. 28447900E8 

MINE-LIKE  TARGET 
MODEL  NO.  2 

0.6m  ANTENNA 

-1.49959862E8 
-9.64557560E7 
-2.09351880E8 
-1 . 91 461 498E8 
-1 . 91 729062E8 

*7. 56869867 E7 
.* 1 . 64502020E8 
*2.2771 1640E8 
±2.89778957  E8 
±4. T4181320E8 

MINE-LIKE  TARGET 
MODEL  NO.  1 

0.15m  ANTENNA 

-6.90468500E7 
-2 . 801 1 1 200E8 
-2.35432000E8 
-5.03972000E7 
-1 . 22024400E8 

n . 31 552700E8 
±2.22590300E8 
i 2.69826400E8 
--2.99420800E8 
-4.04685600E8 

*Real  Part  In  Nepers/s,  Imaginary  Part  In  Hz. 
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APPENDIX  H 


This  appendix  tabulates  the  Identification  results  for  Identi- 
fication of  the  mine-like  target  with  the  small-antenna  system. 


212 


TABLE  29 

cio  CUTER  COEFFICIENTS  USED  IN  THE  PREPROCESSOR 
OF  THE  MICROCOMPUTER  IDENTIFICATION  SYSTEM 
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APPENDIX  I 


This  appendix  gives  detail  descriptions  of  the  APU[68]  and  the 
micro-program  that  implements  the  various  system  control  and  target 
identification  processes. 


A . The  APU  and  Its  Interfa c e 
with  tine  SDK-80~ 

The  connection  diagram  of  the  APU  is  given  in  Figure  58. 


CONNECTION  DIAGRAM 


TOP  Vlf* 


Pi*  1 mp tor 

Figure  50.  Connection  diagram  for  the  APU. 
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Interfacing  the  APU  with  the  SDK-80  requires  the  generation  of 

the  various  control  signals.  In  this  study,  the  APU  is  interfaced 
to  the  SDK-80  as  a memory  location,  and  the  control  signals  are  gen- 
erated as  follows: 

1 .  CS:  chip  select 

The  chip-select  signal  is  generated  by  using  address 
lines  A13,  A14,  and  A15  nf  the  8080A  processor  in  the 
SDK-80.  The  signal  generation  circuit  is  shown  in 
Figure  5S.  The  chip-select  signal  CS  is  low  when  the 
address  lines  A13,  A14,  A15  are  all  high. 


Figure  59.  The  chip-select  (CS)  signal 
for  the  APU. 


2.  C/D:  command/data 

The  C/D  signal  is  tied  directly  to  the  address  line 
A<(>  of  the  8080A  processor. 

3.  TOR:  tied  directly  to  MFMET  of  the  8080A  processor. 

4.  IOW ; tied  directly  to  MFMW  line  of  the  808A  processor. 


5.  PAUSE:  tied  directly  to  the  ready  line  of  the  8080A 

processor. 

6.  CLK:  tied  directly  to  <\>2  of  the  8080A  processor. 

7.  Falk,  SVACK,  SVREq,  RESET,  END:  unused. 

8.  To  eliminate  possible  loading  problem,  two  additional 
siynals  are  generated  to  inhibit  the  ROM's  and  RAM's  of  the 
SDK-80  when  the  APU  is  being  addressed.  These  two  control 
signals  (El  and  E3)  are  given  in  Figure  CO. 


A suninary  of  the  APU  commands  is  given  in  Table  31 


TABLE  31 

APU  COMMAND  SUMMARY 


Command  Cod* 

Command 

□lL 

6 

s 

3 

2 

_n 

0 

Mntmoniw 

HXIO  POINT  SINGLE  PRECISION  ( 

m \ 

1 

1 1 

a 

n 

\ 

0 

0 

SAOO 

Adds  TOS  ro  NOS  Result  io  NOS  Pco  $»*ck 

R ■ 

1 . 

1 . 

0 

> 

1 

0 

1 

SSUB 

Subtracts  TOS  from  NOS  Result  to  NOS  Poo  Stack 

D l 

i i 

a 

i 

1 

’ 

0 

SMUL 

Mwit-pi.es  NOS  Jy  TOS  Result  to  NOS  Poo  Stack 

* 1 

1 . 

J_L 

a 

i 

1 

t 

1 

SOIV 

Divines  NOS  by  TOS  Result  to  NOS  Poo  Stae* 

[ FtXtO  POINT  DOUBLE  PRECISION 

; "T 

0T 

0 

1 

0 

0 

OAOO 

A. ».h  TOS  to  NOS  Result  to  NOS  Poo  Stack 

! A 1 

o t 

t 1 

0 

« 

1 

0 

1 

0SU9 

Subtract*  TOS  t*om  NOS  RrtuH  to  NOS  Pop  Stack 

' a 1 

C 1 

1 J 

0 

» 

1 

1 

0 

OMUL 

Muitipt.es  NOS  bv  TOS  Result  to  NOS  Pop  Stack 

a . 

- 

^ 1 

0 

1 

1 

1 

ODiV 

0'»nPt  NOS  bv  TOS  Result  to  NOS  Pop  Stac* 

floating  point 

a i 

t * 

° 

’ 

0 

0 

0 

0 

1*00 

AtHJ*  TOS  to  NOS  Result  to  NOS  Pop  Stack 

ft 

0 1 

5 1 

) 

0 

0 

0 

1 

FSUB 

Subtracts  TOS  from  NOS  Result  to  NOS  Poo  Stack 

ft  ' 

c , 

0 , 

1 

o 

0 

. 1 

0 

fmul 

Mwiiii»it«  NOS  hv  TOS  P*wM  to  NOS  Pop  Stack 

a 

0 , 

1 

0 

0 

1 

Lj 

few 

0«»>'iet  NOS  ov  TOS  Result  to  NOS  Poo  Stack 

OF  AIVIO  FLOATING  POINT  FUNCTIONS  <21 

, " 1 

0 

•5  1 

0 

0 

0 

i 0 

1 

S0;<  T 

S'iu«'f  Root  ot  TOS  P*wii  io  TOS 

ft 

0 

a 

0 

0 

0 

1 ’ 

0 

Sin 

S<n*  Ot  TOS  ftr*Vult  io  TOS 

» 

c 

0 

0 

o 

0 

1 > 

1 

COS 

Count  o*  TOS  ftaswlt  *n  TOS 

ft  ; 

c 

c . 

0 

0 

1 

: 0 

0 

TAN 

T^nfinnt  o*  TOS  Pa*u*t  i*  TOS 

ft 

0 

0 « 

c 

0 

* 

f 0 

AS»N 

• «¥¥*%»  ot  TOS  Rvtul*  *n  TOS 

ft 

a 

s 

0 

0 

1 

1 1 

0 

At  OS 

t'»w»W  Co*<ri»  ©♦  TOS  H»sutt  »n  TOS 

‘ « 

0 

. C ; 

a 

0 

» 

‘ \ 

1 

AT  AN 

ir***r**  Tini^m  o*  TOS  Result  •«  TOS 

ft  . 

a 

C 1 

0 

1 

0 

: a 

0 

lO<i 

Civn^on  LOt>dr«thrr>  <0»  O*  TOS  Retutt  m TOS. 

ft 

0 

0 

1 

0 

. 0 

1 

UN 

N«tu*ar  Lo^*a#^^^m  li>a\e  »>  o*  TOS  ResuM  m TOS 

ft 

0 

0 

1 

0 

i 

0 

t«» 

E »*K»‘»eAt*ai  t*11)  Ot  TOS  Result  TOS 

ft 

a 

: 

a 

1 

0 

: i 

\ 

PV.O 

NOS  »a*se<J  to  »be  oow#*  «n  TOS  Result  to  NOS  Pop  Stack 

DATA  MANIPULATION  COMMANDS  <31 

* i 

' » 

a 

1 0 

a 

0 

0 

NO* 

Nis  Otter^Tiom 

ft  * 

u 

i 

1 ! 

t i 

\ 

f.«5 

CunvtMS  T OS  from  **oatin»j  0O«mt  to  Srnple  OreC'SiO"  li*e*J  point  format 

• ft 

- 

d i 

i 

; \ 

: 1 

0 

F:  to 

C<>Avert*  TOS  f»om  float. po^t  to  double  precs'OA  ii«e<j  po.«t  format. 

ft 

j 

c • 

i 

i ’ 

\ 

0 

1 

f LTS 

Cnn*#r»|  TOS  from  s»«y»e  Qr»c>s.OA  t.n#d  oo»nt  to  tloet.«q  PO*nt  fotmjt. 

» ft 

* 

0 

• 

; ’ 

•*\ 

1 4 

0 

f l TO 

Converts  TOS  from  douWe  preesron  ttie<J  point  to  t<0attOQ  point 

\ « 

1 

1 

i 

0 

1 0 

0 

c»ss 

Cbanges  si^n  of  smgla  prectSion  fteeii  p ->■  n t opet»nd  on  TOS 

a 

1 

i 

; 0 

\ 1 

0 

0 

CNSO 

Changes  S»Q»s  Ot  double  precision  lived  po*nt  operand  On  TOS. 

ft 

0 

a 

? 

, 0 

i 0 

1 

CH5F 

CbjAi^s  Siqn  of  floatmq  pr>nt  operand  On  TOS 

ft 

! 

* 

i 

• 0 

i 

• 1 

1 

pros 

Push  Single  preesron  fi«e<s  point  Operjnd  on  TOS  to  NOS 

• ft 

0 

* 

* 

: o 

1 1 

1 

1 

PTOO 

Push  double  precision  ti»ed  OO.nt  ope»  snp  on  TOS  to  NOS 

ft 

c 

# 

i 

i o 

1 

1 

PTOF 

Push  boating  pent  operand  on  TQS  to  NOS 

* ft 

» 

t 

1 ’ 

a 

0 

0 

POPS 

Poo  s<ngle  precision  t«ved  pn«nt  oo^rand  born  TOS  NOS  bacomes  TOS 

o 

*J 

• 

t 

’ 0 

. 0 

0 

POPO 

Pop  dituttV  preciion  lived  pomt  ooarant)  »rom  TOS  NOS  become*  TOS, 

, M 

•• 

t 

1 1 

a 

0 

0 

POPF 

Pup  II, sating  ooint  operand  t»om  TOS  NOS  trecomes  TOS 

ft 

1 

• 

i 

! i 

0 

i 0 

1 

j xO<S 

E »cnang»  ungi*  preemun  lived  point  oueranos  TOS  and  NOS 

ft 

1 

* 

i 

0 

t 0 

1 

I XC*0 

CvChange  »toub*r  precdipn  lived  oorni  operands  TOS  and  NOS 

ft 

'* 

c 

t 

1 1 

c 

0 

1 

KCMF 

E»Change  floating  pomt  operands  TOS  and  NOS 

ft 

2 

* 

t 

1 

0 

1 

0 

l PuPl 

Push  lloatmg  pomt  constant  "•'*  nntp  TOS  Previous  TOS  becomes  NOS. 

A*  *ftf  ■»  Top  C 1 !tae*  V*>S  «s  Ne*t  0*  Stack 

) A<|  ,e»»  ****!.•• , »*•  -nr  lw>«ct>vA|  <W'o*  ih*  ot  the  stack  0**r  m*  result  can  be  counted  on  to  bt  valid  upon  eomnumd 

v* 

3 * V«“H  »t.m*«t*a*  •"  FIXO,  fs.TS,  *LTOt  ’*OVtr«  that  floating  pomt  rta»*  format  ba  specified  Itommand  pit*  6 »ndt 

»MI  he  Ol 


B. 


The  Microprogram  for  the 
Microcomputer  System 


The  various  functions  of  the  microcomputer  systems  are  imp! 
mented  as  commands  listed  in  Table  32.  The  microprogram  that 
implements  the  various  commands  is  given  following  Table  32. 


TAB'  E 32 

TABLE  OF  COMMANDS  IMPLEMENTED  IN 
THE  MICROCOMPUTER  SYSTEM 


COMMAND  DESCRIPTIONS 

CODE 


F Restart 

A Display  a memory  waveform  on  the 

Oscil Iosco pe 

8 Branch  to  another  ROM 

7 Displays  data  values  of  Interest 

6 Change  or  enter  the  number  of  waveform 

taken  to  form  an  average  waveform 

H Change  or  enter  number  of  samples  per 

waveform 

2 Change  or  enter  the  sequence  number 

for  the  next  waveform 

1 Initiate  Recording  sequence  via  the 

push  button 

H Record  a waveform,  and  perform  the 

identification  process 

B Dump  memory  onto  tape 

D Check  if  data  transmission  to 

recorder  is  error- free 
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Details  In  Obtaining  the  Second 
Set'  oflda'ntlYI' cation  Data  With 
tjjeTWT crbc'omp'u ter  System 

1-  Pi fference  equation  .coefficients 

Table  33  tabulates  the  difference  equation  coefficients  used 
In  the  microcomputer  system  for  identification  of  the  mine-like 
target. 


2.  Detection  thresholds  (Rj0«3O  cm) 
tm  RANGE  « 14.49  (TB) 

MAX  RANGE  « 12,82 
Em  RANGE  “ 0.639,  0.862 

The  same  detection  thresholds  were  used  for  Identification  of 
the  mine -like  tarqet  In  both  ground  conditions. 

3.  Identification  thresholds  (Rjp»30  cm) 

Table  34  tabulates  the  Identification  thresholds  for  identifi- 
cation of  the  mine-like  target  In  the  two  different  ground  conditions. 
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TABLE  34 

IDENTIFICATION  THRESHOLDS  FOR  THE  IDENTIFICATION  OF  THE 
MINE-LIKE  TARGET  WITH  THE  MICROCOMPUTER  SYSTEM 
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“T01 

T.4TB 

"T02 
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03 
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f T04 
T=7Tb 
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f,T 

T05 

T=8Tb 
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V 

CONDITION 

1 

(TUNED) 

.522 

.301 

.615 

.677 

.780 

.79177 

CONDITION 

2 

(UNTUNED) 

.500 

.500 

.250 

.125 

.500 

.40000 

» 


APPENDIX  J 


Table  3b  tabulates  the  extracted  resonances  and  their  corres 
ponding  residues  of  the  thin-wire  (30cm  long,  5cm  deep)  waveform. 


Table  35 

EXTRACTED  RESONANCES  FROM  THE  30cm  LONG,  5cm  DEEP  THIN  WIRE 
ANTENNA  LOCATION  = CENTER  OF  WIRE 


POLE  * 
(REAL) 

POLE 

(IMAG) 

RESIDUE 

(REAL) 

RESIDUE 

(IMAG) 

- . 1 308792E  9 

- . 2243403E  9 

- . 1 708435E  0 

. 6741 1 98E-1 

- . 4685945E  9 

28008 13E  9 

. 1377038E  1 

- . 2542809E  1 

-.21 5600 IE  9 

. '260306E  8 

. 2210962E  0 

. 1 263321 E 1 

4685945E  9 

280081 3E  9 

. 1377038E  1 

. 2542809E  1 

- . 2156001 E 9 

7260306E  8 

. 2210962E  0 

- . 1 263321 E 1 

- . 1 70807eE  9 

. 3030841 E 9 

. 545901 4E- 1 

.2335131E  0 

1209931 E 9 

. M26777E  9 

-.96311 55E  0 

-.2862594E  1 

- . 1 209931 E 9 

-.1  1 267 7 7 E 9 

- . 9631 167E  0 

. 2862594E  1 

- . 1 7 0807 3 E 9 

-. 3030841 E 9 

. 5459009E-1 

2335131 E 0 

- . 1 308792E  9 

. 2243403E  9 

-. 1708435E  0 

-.6741 193E-1 

- . 6272871 C 8 

. 1839096E-2 

2458968E-1 

-. 4475425E-7 

PARAMETERS:  N- 

11,  IBS-7,  M=2, 

1 1 1=4,  , =0. 391E-2 

Two  pairs  of  poles  appear  to  have  dominant  residues,  namely  the 
pole  pairs  at  280  MHz  and  112  MHz.  However,  the  pole  pair  at  280  MHz 
has  a much  more  negative  real  part,  thus,  in  the  late-time  region, 

* orl"y  the  pole  pair  at.  112  MHz  is  dominant.  .. 

Table  36  tahulai.es  the  extracted  resonances  and  their  correspond- 
ing residues  from  the  30cm  long  on-surface  wire. 

In  the  early-time  region,  the  pole  pair  at  302  MHz  is  dominant. 


* Real  Part  in  Nepers/s.  Imaginary  part  in  Hz. 
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TABLE  36 

EXTRACTED  RESONANCES  AND  FROM  THE  ON-SURFACE  THIN  WIRE 
(30cm  LONG).  ANTENNA  LOCATION  * CENTER  OF  WIRE 


POLE* 

(REAL) 

POLE 

(IMAG) 

RESIDUE 

(REAL) 

RESIDUE 

(IMAG) 

1462028EI0 

.OOOOOOOE  1 

. 24691 65E  1 

•6242044E-7 

-.3449277E  9 

-.3023622E  9 

, 1591370E  1 

- .1711432E  1 

- .3449277E  9 

.3C23622E  9 

. 1591370E  1 

.1711432E  1 

-.1003984E  9 

-.2404006E  9 

- . 1003876E  0 

. 1049173E-1 

-.2344375E  9 

.4529953E  9 

2735412E  0 

.1113815E-1 

-.4641263E  8 

.6872750E  8 

- . 6308170E-1 

. 7273762E-1 

- . 2344375E  9 

-.4529953E  9 

-2.035411E  0 

- . 1 113816E-1 

- . 1003984E  9 

. 2404006E  9 

1003877E  0 

1049174E-1 

- . 3901 977 E 8 

- . 1 401 51 3 E 9 

-.2077B24E-1 

-.2096944E-1 

-.3901977E  8 

.1401513E  9 

-.2077823E-1 

. 2096944E-1 

- .3232055E  9 

.3787546E  9 

-.3689420E  0 

-.6241565E  0 

- . 3232055E  9 

3787546E  9 

-.3680420E  0 

.6241 565E  0 

-.4641263E  8 

-.6872750E  8 

-.6308171E-1 

-.727376U-1 

PARAMETERS: 

N=13,  IBS=5,  M=2, 

111=1,  c=0.3839075E 
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♦Real  part  in  Nepers/s.  Imaginary  part  In  Hz. 


